var/obsproc/MAP_plot/Dir_map/ffdduv.F

   1       SUBROUTINE FFDDUV(F,D,U,V,YLON,XLONC,PHIC,IPROJ,II,ID)
   2 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   3 C                                                                 C
   4 C    WIND SPEED AND DIRECTION <-------> WIND COMPONENTS U AND V   C
   5 C                                                                 C
   6 C    ID = 1,      F, D ------> U, V                               C
   7 C    ID = 2,      U, V ------> F, D                               C
   8 C                                                                 C
   9 C    II = 0,      SAME UNIT IN BOTH INPUT AND OUTPUT              C
  10 C    II = 1,      INPUT UNIT = KNOT,     OUTPUT UNIT = M/S        C
  11 C    II = 2,      INPUT UNIT = M/S,      OUTPUT UNIT = KNOT       C
  12 C                                                                 C
  13 C  ** The definition of wind components at the South/North poles  C
  14 C     is based on the YLON provided by the station information as C
  15 C     usual:i.e. the D = 0 degree is the wind coming from NORTH   C
  16 C     along the longitude YLON.                                   C
  17 C                             Y.-R. Guo   02/14/2003              C
  18 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  19 C
  20       CHARACTER*6 NPROJ(3),IPROJ
  21       DATA NPROJ/'LAMCON','POLSTR','MERCAT'/
  22 C
  23 c      print '("F,D,U,V,YLON,XLONC,PHIC,IPROJ,II,ID:")'
  24 c      print '(7f10.2,1x,A,2I3)',F,D,U,V,YLON,XLONC,PHIC,IPROJ,II,ID
  25
  26       XN = 0.716
  27       IF(IPROJ .EQ.  NPROJ(2)) XN = 1.0
  28       DEGRAN = 3.1415927/180.
  29 C
  30       SKN = 0.5144
  31 C
  32       GO TO (10,20,40) ID
  33 C
  34 10    CONTINUE
  35 C
  36 C (1) TO CONVERT THE WIND FROM VELOCITY FORM TO COMPONENT FORM ON
  37 C     EARTH, AND THEN TO LAMBERT CONFORMAL PROJECTION, OR POLAR
  38 C     STEREOGRAPHIC PROJECTION, U AND V COMPONENTS.  CHANGEING
  39 C     PROJECTIONS REQUIRES ONLY A CHANGE OF THE CONE CONSTANT, XN
  40 C     -- THE EQUATIONS REMAIN THE SAME.
  41 C
  42       AEARTH = D*DEGRAN
  43       UEARTH = -F*SIN(AEARTH)
  44       VEARTH = -F*COS(AEARTH)
  45 C     FOR CONVERSION TO GRID COORDINATES,
  46 C     SEE PROGRAM DATAMAP, SUBR VECT, AND
  47 C     ANTHES' METEO. 597 NOTES, EQUA. 2.23, 2.25, 2.28.
  48 C     ALSO, CONVERT FROM KNOTS TO M/S.
  49       XLONRT = XLONC-YLON
  50       IF(XLONRT .GT. 180.) XLONRT=XLONRT-360.
  51       IF(XLONRT .LT.-180.) XLONRT=XLONRT+360.
  52       ANG=XLONRT*XN*DEGRAN
  53 C
  54 C   FOR MERCATOR PROJECTION, THE WINDS ARE AS IN EARTH COORDINATES
  55 C
  56       IF(IPROJ.EQ.NPROJ(3)) ANG=0.
  57       IF(PHIC.LT.0.0) ANG=-ANG
  58       U = VEARTH*SIN(ANG) + UEARTH*COS(ANG)
  59       V = VEARTH*COS(ANG) - UEARTH*SIN(ANG)
  60 C
  61       GO TO 30
  62 C
  63 20    CONTINUE
  64 C
  65       XLONRT = XLONC-YLON
  66       IF(XLONRT .GT. 180.) XLONRT=XLONRT-360.
  67       IF(XLONRT .LT.-180.) XLONRT=XLONRT+360.
  68       ANG=XLONRT*XN*DEGRAN
  69 C
  70 C   FOR MERCATOR PROJECTION, THE WINDS ARE AS IN EARTH COORDINATES
  71 C
  72       IF(IPROJ.EQ.NPROJ(3)) ANG=0.
  73       IF(PHIC.LT.0.0) ANG=-ANG
  74       UE = U*COS(ANG) - V*SIN(ANG)
  75       VE = U*SIN(ANG) + V*COS(ANG)
  76 C
  77       F = SQRT(UE*UE + VE*VE)
  78       IF (F.EQ.0.0) THEN
  79          D = 99999.0
  80          GO TO 30
  81       ENDIF
  82 C
  83       IF (VE.EQ.0.0) THEN
  84          IF (UE.GT.0.0) D = 270.
  85          IF (UE.LT.0.0) D =  90.
  86 C
  87       ELSE
  88       AEARTH = ATAN(UE/VE)/DEGRAN
  89       IF (UE.LE.0.0 .AND. VE.LE.0.0 ) D = AEARTH
  90       IF (UE.LE.0.0 .AND. VE.GE.0.0 ) D = AEARTH + 180.0
  91       IF (UE.GE.0.0 .AND. VE.GE.0.0 ) D = AEARTH + 180.0
  92       IF (UE.GE.0.0 .AND. VE.LE.0.0 ) D = AEARTH + 360.0
  93       ENDIF
  94 C
  95 30    CONTINUE
  96 C
  97       IF (II.EQ.0) GO TO 40
  98       IF (II.EQ.1) THEN
  99          IF (ID.EQ.1) THEN
 100             U = U*SKN
 101             V = V*SKN
 102          ELSE IF (ID.EQ.2) THEN
 103             F = F*SKN
 104          ENDIF
 105       ELSE IF (II.EQ.2) THEN
 106          IF (ID.EQ.1) THEN
 107             U = U/SKN
 108             V = V/SKN
 109          ELSE IF (ID.EQ.2) THEN
 110             F = F/SKN
 111          ENDIF
 112       ENDIF
 113 C
 114 40    RETURN
 115       END