74 |
ELSE |
ELSE |
75 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
76 |
CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid ) |
CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid ) |
77 |
|
DO bj = myByLo(myThid), myByHi(myThid) |
78 |
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
79 |
|
DO j=1,sNy |
80 |
|
DO i=1,sNx |
81 |
|
fu(i,j,bi,bj) = fu(i,j,bi,bj)/(delZ(1)*rhonil)*0.1 |
82 |
|
fu(i,j,bi,bj) = 0.1/(delZ(1)*rhonil) |
83 |
|
ENDDO |
84 |
|
ENDDO |
85 |
|
ENDDO |
86 |
|
ENDDO |
87 |
_END_MASTER(myThid) |
_END_MASTER(myThid) |
88 |
ENDIF |
ENDIF |
89 |
C |
C |
102 |
ELSE |
ELSE |
103 |
_BEGIN_MASTER(myThid) |
_BEGIN_MASTER(myThid) |
104 |
CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid ) |
CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid ) |
105 |
|
DO bj = myByLo(myThid), myByHi(myThid) |
106 |
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
107 |
|
DO j=1,sNy |
108 |
|
DO i=1,sNx |
109 |
|
fv(i,j,bi,bj) = fv(i,j,bi,bj)/(delZ(1)*rhonil)*0. |
110 |
|
ENDDO |
111 |
|
ENDDO |
112 |
|
ENDDO |
113 |
|
ENDDO |
114 |
_END_MASTER(myThid) |
_END_MASTER(myThid) |
115 |
ENDIF |
ENDIF |
116 |
C |
C |
117 |
_EXCH_XY_R4(fu , myThid ) |
_EXCH_XY_R4(fu , myThid ) |
118 |
_EXCH_XY_R4(fv , myThid ) |
_EXCH_XY_R4(fv , myThid ) |
119 |
|
|
120 |
CcnhDebugStarts |
CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid) |
121 |
Cdbg WRITE(0,*) ' distY = ', distY |
CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid) |
|
Cdbg WRITE(0,*) ' ly = ', lY |
|
|
Cdbg WRITE(0,*) ' tauMax= ', tauMax |
|
|
Cdbg CALL PLOT_FIELD_XYR8( fu, 'INI_FORCING FU',1,myThid) |
|
|
Cdbg STOP 'INI_FORCING' |
|
|
CcnhDebugEnds |
|
122 |
|
|
123 |
RETURN |
RETURN |
124 |
END |
END |