| File |
Line number |
Procedure |
Code |
|
./pkg/exch2/w2_eeboot.F |
40 |
W2_EEBOOT |
INTEGER ii, I, J |
|
./pkg/exch2/w2_eeboot.F |
62 |
W2_EEBOOT |
DO J = 1,W2_maxNeighbours |
|
./pkg/exch2/w2_eeboot.F |
63 |
W2_EEBOOT |
exch2_neighbourId(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
64 |
W2_EEBOOT |
exch2_opposingSend(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
66 |
W2_EEBOOT |
exch2_pij(ii,J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
68 |
W2_EEBOOT |
exch2_oi(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
69 |
W2_EEBOOT |
exch2_oj(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
70 |
W2_EEBOOT |
exch2_iLo(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
71 |
W2_EEBOOT |
exch2_iHi(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
72 |
W2_EEBOOT |
exch2_jLo(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
73 |
W2_EEBOOT |
exch2_jHi(J,I) = 0 |
|
./pkg/exch2/w2_map_procs.F |
102 |
W2_MAP_PROCS |
IF ( J .NE. nSx*nSy*nPx*nPy ) THEN |
|
./pkg/exch2/w2_map_procs.F |
149 |
W2_MAP_PROCS |
DO J=1,exch2_nNeighbours(myTileId) |
|
./pkg/exch2/w2_map_procs.F |
151 |
W2_MAP_PROCS |
jj = exch2_neighbourId(J,myTileId) |
|
./pkg/exch2/w2_map_procs.F |
155 |
W2_MAP_PROCS |
& ' NEIGHBOUR',J,' = TILE', exch2_neighbourId(J,myTileId), |
|
./pkg/exch2/w2_map_procs.F |
156 |
W2_MAP_PROCS |
& ' (n=', exch2_opposingSend(J,myTileId), ') Comm = MSG', |
|
./pkg/exch2/w2_map_procs.F |
157 |
W2_MAP_PROCS |
& ' (PROC=',W2_tileProc(exch2_neighbourId(J,myTileId)),')' |
|
./pkg/exch2/w2_map_procs.F |
162 |
W2_MAP_PROCS |
& ' NEIGHBOUR',J,' = TILE', exch2_neighbourId(J,myTileId), |
|
./pkg/exch2/w2_map_procs.F |
163 |
W2_MAP_PROCS |
& ' (n=', exch2_opposingSend(J,myTileId), ') Comm = PUT', |
|
./pkg/exch2/w2_map_procs.F |
164 |
W2_MAP_PROCS |
& ' (PROC=',W2_tileProc(exch2_neighbourId(J,myTileId)),')' |
|
./pkg/exch2/w2_map_procs.F |
167 |
W2_MAP_PROCS |
W2_myCommFlag(J,bi,bj) = commFlag |
|
./pkg/exch2/w2_map_procs.F |
45 |
W2_MAP_PROCS |
INTEGER I, J |
|
./pkg/exch2/w2_map_procs.F |
61 |
W2_MAP_PROCS |
DO J=1,W2_maxNeighbours |
|
./pkg/exch2/w2_map_procs.F |
62 |
W2_MAP_PROCS |
W2_myCommFlag(J,bi,bj) = ' ' |
|
./pkg/exch2/w2_map_procs.F |
76 |
W2_MAP_PROCS |
J = 0 |
|
./pkg/exch2/w2_map_procs.F |
85 |
W2_MAP_PROCS |
jj = J/(nSx*nPx) |
|
./pkg/exch2/w2_map_procs.F |
86 |
W2_MAP_PROCS |
ii = MOD(J,nSx*nPx) |
|
./pkg/exch2/w2_map_procs.F |
98 |
W2_MAP_PROCS |
J = J + 1 |
|
./pkg/exf/exf_mapfields.F |
342 |
EXF_MAPFIELDS |
saltFlux(I,J,bi,bj) = saltflx(I,J,bi,bj) |
|
./pkg/exf/exf_mapfields.F |
352 |
EXF_MAPFIELDS |
exf_iceFraction(I,J,bi,bj) = |
|
./pkg/exf/exf_mapfields.F |
353 |
EXF_MAPFIELDS |
& MIN( MAX(exf_iceFraction(I,J,bi,bj),zeroRS), oneRS ) |
|
./pkg/fizhi/fizhi_utils.F |
101 |
STRIPINT |
DO 100 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
103 |
STRIPINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
106 |
STRIPINT |
DO 200 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
108 |
STRIPINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
111 |
STRIPINT |
B(Lena+I,J) = A(Lena+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
129 |
PASTE |
DO 100 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
131 |
PASTE |
A(I+OFFSET-1,J) = B(I,J) |
|
./pkg/fizhi/fizhi_utils.F |
147 |
PSTBMP |
DO 100 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
149 |
PSTBMP |
A(I+OFFSET-1,J) = A(I+OFFSET-1,J) + B(I,J) |
|
./pkg/fizhi/fizhi_utils.F |
166 |
STRINT |
DO 100 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
168 |
STRINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
171 |
STRINT |
DO 200 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
173 |
STRINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
176 |
STRINT |
B(Lena+I,J) = A(Lena+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
72 |
STRIP |
DO 100 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
74 |
STRIP |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
77 |
STRIP |
DO 200 J=1,L |
|
./pkg/fizhi/fizhi_utils.F |
79 |
STRIP |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
82 |
STRIP |
B(Lena+I,J) = A(Lena+OFFSET-1,J) |
|
./pkg/fizhi/update_earth_exports.F |
888 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
889 |
GETALB |
ALBOC(I,J) = A0 + (A1 + (A2 + A3*cosz(I,J))*cosz(I,J))*cosz(I,J) |
|
./pkg/fizhi/update_earth_exports.F |
890 |
GETALB |
ALBVR(I,J,bi,bj) = ALBSI*FRACI(I,J) + ALBOC(I,J)*(ONE-FRACI(I,J)) |
|
./pkg/fizhi/update_earth_exports.F |
891 |
GETALB |
ALBNR(I,J,bi,bj) = ALBVR(I,J,bi,bj) |
|
./pkg/fizhi/update_earth_exports.F |
892 |
GETALB |
ALBVF(I,J,bi,bj) = ALBSI * FRACI(I,J) + OCNALB * (ONE-FRACI(I,J)) |
|
./pkg/fizhi/update_earth_exports.F |
893 |
GETALB |
ALBNF(I,J,bi,bj) = ALBVF(I,J,bi,bj) |
|
./pkg/fizhi/update_earth_exports.F |
909 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
917 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
922 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
929 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
934 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
941 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
946 |
GETALB |
DO J=1,JM |
|
./pkg/fizhi/update_earth_exports.F |
953 |
GETALB |
DO J=1,JM |
|
./pkg/flt/exch2_send_put_vec.F |
67 |
EXCH2_SEND_PUT_VEC_RL |
INTEGER I,J,sBi,sBj, nN, thisTile, tgT, tgN, dir, sDir |
|
./pkg/ggl90/ggl90_calc.F |
158 |
GGL90_CALC |
recip_hFacI(I,J,K)=0. _d 0 |
|
./pkg/ggl90/ggl90_calc.F |
160 |
GGL90_CALC |
& recip_hFacI(I,J,K)=1. _d 0/hFac |
|
./pkg/ggl90/ggl90_init_varia.F |
104 |
GGL90_INIT_VARIA |
DO J=1-OLy,sNy+OLy |
|
./pkg/ggl90/ggl90_init_varia.F |
123 |
GGL90_INIT_VARIA |
DO J=1-OLy,sNy+OLy |
|
./pkg/ggl90/ggl90_init_varia.F |
148 |
GGL90_INIT_VARIA |
DO J=1-OLy,sNy+OLy |
|
./pkg/ggl90/ggl90_init_varia.F |
150 |
GGL90_INIT_VARIA |
GGL90TKE(I,J,K,bi,bj) = MAX(GGL90TKE(I,J,K,bi,bj), |
|
./pkg/ggl90/ggl90_init_varia.F |
151 |
GGL90_INIT_VARIA |
& GGL90TKEmin)*maskC(I,J,K,bi,bj) |
|
./pkg/ggl90/ggl90_init_varia.F |
37 |
GGL90_INIT_VARIA |
INTEGER I, J, K |
|
./pkg/ggl90/ggl90_init_varia.F |
48 |
GGL90_INIT_VARIA |
DO J=1-OLy,sNy+OLy |
|
./pkg/ggl90/ggl90_init_varia.F |
50 |
GGL90_INIT_VARIA |
GGL90viscArU(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
51 |
GGL90_INIT_VARIA |
GGL90viscArV(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
52 |
GGL90_INIT_VARIA |
GGL90diffKr(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
54 |
GGL90_INIT_VARIA |
GGL90TKE(I,J,K,bi,bj)=GGL90eps*maskC(I,J,K,bi,bj) |
|
./pkg/ggl90/ggl90_init_varia.F |
57 |
GGL90_INIT_VARIA |
GGL90TKE(I,J,K,bi,bj)=GGL90TKEmin*maskC(I,J,K,bi,bj) |
|
./pkg/ggl90/ggl90_init_varia.F |
77 |
GGL90_INIT_VARIA |
DO J=1-OLy,sNy+OLy |
|
./pkg/ggl90/ggl90_init_varia.F |
79 |
GGL90_INIT_VARIA |
IDEMIX_E(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
80 |
GGL90_INIT_VARIA |
IDEMIX_V0(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
81 |
GGL90_INIT_VARIA |
IDEMIX_tau_d(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
86 |
GGL90_INIT_VARIA |
DO J=1-OLy,sNy+OLy |
|
./pkg/ggl90/ggl90_init_varia.F |
88 |
GGL90_INIT_VARIA |
IDEMIX_F_B(I,J,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
89 |
GGL90_INIT_VARIA |
IDEMIX_F_S(I,J,bi,bj) = 0. _d 0 |
|
./pkg/grdchk/grdchk_get_position.F |
187 |
GRDCHK_GET_POSITION |
& .and. i.eq. OB_Iw(J,bi,bj) ) then |
|
./pkg/grdchk/grdchk_get_position.F |
195 |
GRDCHK_GET_POSITION |
& .and. i.eq. OB_Ie(J,bi,bj) ) then |
|
./pkg/grdchk/grdchk_loc.F |
240 |
GRDCHK_LOC |
& .and. i.eq. OB_Iw(J,bi,bj) ) then |
|
./pkg/grdchk/grdchk_loc.F |
242 |
GRDCHK_LOC |
itmp = OB_Iw(J,bi,bj) |
|
./pkg/grdchk/grdchk_loc.F |
249 |
GRDCHK_LOC |
& .and. i.eq. OB_Ie(J,bi,bj) ) then |
|
./pkg/grdchk/grdchk_loc.F |
251 |
GRDCHK_LOC |
itmp = OB_Ie(J,bi,bj) |
|
./pkg/icefront/icefront_init_fixed.F |
100 |
ICEFRONT_INIT_FIXED |
& K_icefront(I,J,bi,bj) = K |
|
./pkg/icefront/icefront_init_fixed.F |
110 |
ICEFRONT_INIT_FIXED |
DO J = 1, sNy |
|
./pkg/icefront/icefront_init_fixed.F |
113 |
ICEFRONT_INIT_FIXED |
IF (ICEFRONTlength(I,J,bi,bj) .GT. 0. _d 0) THEN |
|
./pkg/icefront/icefront_init_fixed.F |
115 |
ICEFRONT_INIT_FIXED |
IsInterface=Isinterface + K_icefront(I,J,bi,bj) |
|
./pkg/icefront/icefront_init_fixed.F |
37 |
ICEFRONT_INIT_FIXED |
INTEGER I, J, K, bi, bj |
|
./pkg/icefront/icefront_init_fixed.F |
85 |
ICEFRONT_INIT_FIXED |
DO J = 1-OLy, sNy+OLy |
|
./pkg/icefront/icefront_init_fixed.F |
87 |
ICEFRONT_INIT_FIXED |
R_icefront(I,J,bi,bj) = ABS(R_icefront(I,J,bi,bj)) |
|
./pkg/icefront/icefront_init_fixed.F |
95 |
ICEFRONT_INIT_FIXED |
DO J = 1-OLy, sNy+OLy |
|
./pkg/icefront/icefront_init_fixed.F |
99 |
ICEFRONT_INIT_FIXED |
IF ( R_icefront(I,J,bi,bj) .GT. ABS(rF(K))) |
|
./pkg/icefront/icefront_init_varia.F |
31 |
ICEFRONT_INIT_VARIA |
INTEGER I, J, K, bi, bj |
|
./pkg/icefront/icefront_init_varia.F |
37 |
ICEFRONT_INIT_VARIA |
DO J = 1, sNy |
|
./pkg/icefront/icefront_thermodynamics.F |
100 |
ICEFRONT_THERMODYNAMICS |
& *ICEFRONTheatTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
103 |
ICEFRONT_THERMODYNAMICS |
eps1 = rUnit2mass*HeatCapacity_Cp*ICEFRONTheatTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
104 |
ICEFRONT_THERMODYNAMICS |
eps2 = rUnit2mass*ICEFRONTlatentHeat*ICEFRONTsaltTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
106 |
ICEFRONT_THERMODYNAMICS |
& *ICEFRONTsaltTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
113 |
ICEFRONT_THERMODYNAMICS |
tLoc = theta(I,J,K,bi,bj) |
|
./pkg/icefront/icefront_thermodynamics.F |
114 |
ICEFRONT_THERMODYNAMICS |
sLoc = MAX(salt(I,J,K,bi,bj), 0. _d 0) |
|
./pkg/icefront/icefront_thermodynamics.F |
121 |
ICEFRONT_THERMODYNAMICS |
IF ( K .EQ. K_icefront(I,J,bi,bj)) THEN |
|
./pkg/icefront/icefront_thermodynamics.F |
122 |
ICEFRONT_THERMODYNAMICS |
pLoc = 0.5*(ABS(R_icefront(I,J,bi,bj))+ABS(rF(K))) |
|
./pkg/icefront/icefront_thermodynamics.F |
125 |
ICEFRONT_THERMODYNAMICS |
& (R_icefront(I,J,bi,bj)-pLoc) |
|
./pkg/icefront/icefront_thermodynamics.F |
126 |
ICEFRONT_THERMODYNAMICS |
& / R_icefront(I,J,bi,bj) |
|
./pkg/icefront/icefront_thermodynamics.F |
146 |
ICEFRONT_THERMODYNAMICS |
FreshWaterFlux(I,J) = maskC(I,J,K,bi,bj) * |
|
./pkg/icefront/icefront_thermodynamics.F |
150 |
ICEFRONT_THERMODYNAMICS |
HeatFlux(I,J) = maskC(I,J,K,bi,bj) * HeatCapacity_Cp * |
|
./pkg/icefront/icefront_thermodynamics.F |
151 |
ICEFRONT_THERMODYNAMICS |
& ( -rUnit2mass*ICEFRONTheatTransCoeff(I,J) + |
|
./pkg/icefront/icefront_thermodynamics.F |
152 |
ICEFRONT_THERMODYNAMICS |
& FreshWaterFlux(I,J) ) * ( thetaFreeze - tLoc ) |
|
./pkg/icefront/icefront_thermodynamics.F |
155 |
ICEFRONT_THERMODYNAMICS |
icefront_TendT(i,j,K,bi,bj) = - HeatFlux(I,J)* eps5 |
|
./pkg/icefront/icefront_thermodynamics.F |
156 |
ICEFRONT_THERMODYNAMICS |
icefront_TendS(i,j,K,bi,bj) = FreshWaterFlux(I,J) * |
|
./pkg/icefront/icefront_thermodynamics.F |
171 |
ICEFRONT_THERMODYNAMICS |
& * (ABS(R_icefront(I,J,bi,bj))-ABS(rF(K))) |
|
./pkg/icefront/icefront_thermodynamics.F |
175 |
ICEFRONT_THERMODYNAMICS |
& * (ABS(R_icefront(I,J,bi,bj))-ABS(rF(K))) |
|
./pkg/icefront/icefront_thermodynamics.F |
181 |
ICEFRONT_THERMODYNAMICS |
HeatFlux (I,J) = 0. _d 0 |
|
./pkg/icefront/icefront_thermodynamics.F |
182 |
ICEFRONT_THERMODYNAMICS |
FreshWaterFlux(I,J) = 0. _d 0 |
|
./pkg/icefront/icefront_thermodynamics.F |
60 |
ICEFRONT_THERMODYNAMICS |
INTEGER I,J,K |
|
./pkg/icefront/icefront_thermodynamics.F |
88 |
ICEFRONT_THERMODYNAMICS |
DO J = 1, sNy |
|
./pkg/icefront/icefront_thermodynamics.F |
92 |
ICEFRONT_THERMODYNAMICS |
IF( ICEFRONTlength(I,J,bi,bj) .GT. 0. _d 0 |
|
./pkg/icefront/icefront_thermodynamics.F |
93 |
ICEFRONT_THERMODYNAMICS |
& .AND. K .LE. K_icefront(I,J,bi,bj) ) THEN |
|
./pkg/icefront/icefront_thermodynamics.F |
94 |
ICEFRONT_THERMODYNAMICS |
ICEFRONTheatTransCoeff(I,J) = 1.0 _d -02 |
|
./pkg/icefront/icefront_thermodynamics.F |
95 |
ICEFRONT_THERMODYNAMICS |
& *abs(wVEL(I,J,K,bi,bj)) |
|
./pkg/icefront/icefront_thermodynamics.F |
97 |
ICEFRONT_THERMODYNAMICS |
ICEFRONTheatTransCoeff(I,J) = max |
|
./pkg/icefront/icefront_thermodynamics.F |
98 |
ICEFRONT_THERMODYNAMICS |
& (ICEFRONTheatTransCoeff(I,J),1. _d -04) |
|
./pkg/icefront/icefront_thermodynamics.F |
99 |
ICEFRONT_THERMODYNAMICS |
ICEFRONTsaltTransCoeff(I,J) = 5.05 _d -3 |
|
./pkg/kl10/kl10_calc.F |
105 |
KL10_CALC |
IF ( (rhoS(K).LT.rhoS(K-1)).AND.(maskC(I,J,K,bi |
|
./pkg/kl10/kl10_calc.F |
122 |
KL10_CALC |
KLdiffKr(I,J,1,bi,bj) = MAX(KLviscAr(I,J,1,bi,bj), |
|
./pkg/kl10/kl10_calc.F |
124 |
KL10_CALC |
& diffKr(I,J,1,bi,bj) ) |
|
./pkg/kl10/kl10_calc.F |
165 |
KL10_CALC |
KLviscAr(I,J,K,bi,bj) = MAX(KTemp,viscArNr(K)) |
|
./pkg/kl10/kl10_calc.F |
166 |
KL10_CALC |
KLdiffKr(I,J,K,bi,bj) = MAX(KTemp, |
|
./pkg/kl10/kl10_calc.F |
168 |
KL10_CALC |
& diffKr(I,J,K,bi,bj) ) |
|
./pkg/kl10/kl10_calc.F |
187 |
KL10_CALC |
IF (hfacW(I+di,J-1,K,bi,bj).GT.0.000001) THEN |
|
./pkg/kl10/kl10_calc.F |
188 |
KL10_CALC |
dz = 0.5*(drF(K)*hfacW(I+di,J-1,K,bi,bj) |
|
./pkg/kl10/kl10_calc.F |
189 |
KL10_CALC |
& +drF(Km1)*hfacW(I+di,J-1,Km1,bi,bj)) |
|
./pkg/kl10/kl10_calc.F |
191 |
KL10_CALC |
tempu = (uVel(I+di,J-1,Km1,bi,bj)-uVel(I+di,J |
|
./pkg/kl10/kl10_calc.F |
193 |
KL10_CALC |
epsw(di)=tempu*tempu*KLviscAr(I+di,J-1,K,bi |
|
./pkg/kl10/kl10_calc.F |
201 |
KL10_CALC |
dz=0.5*(drF(Km1)*hfacW(I+di,J-1,Km1,bi ,bj)) |
|
./pkg/kl10/kl10_calc.F |
206 |
KL10_CALC |
IF (hfacS(I-1,J+di,K,bi,bj).GT.0.000001) THEN |
|
./pkg/kl10/kl10_calc.F |
207 |
KL10_CALC |
dz = 0.5*(drF(K)*hfacS(I-1,J+di,K,bi,bj) |
|
./pkg/kl10/kl10_calc.F |
208 |
KL10_CALC |
& +drF(Km1)*hfacS(I-1,J+di,Km1,bi,bj)) |
|
./pkg/kl10/kl10_calc.F |
210 |
KL10_CALC |
tempu = (vVel(I-1,J+di,Km1,bi,bj)-vVel(I-1,J |
|
./pkg/kl10/kl10_calc.F |
212 |
KL10_CALC |
epss(di)=tempu*tempu*KLviscAr(I-1,J+di,K,bi |
|
./pkg/kl10/kl10_calc.F |
220 |
KL10_CALC |
dz=+0.5*(drF(Km1)*hfacS(I-1,J+di,Km1 ,bi,bj)) |
|
./pkg/kl10/kl10_calc.F |
242 |
KL10_CALC |
KLeps(I-1,J-1,Km1,bi,bj) = 0.5*(ep+ec) |
|
./pkg/kl10/kl10_calc.F |
244 |
KL10_CALC |
KLeps(I-1,J-1,Km1,bi,bj) = ec |
|
./pkg/kl10/kl10_calc.F |
259 |
KL10_CALC |
KLeps(I-1,J-1,Nr,bi,bj) =ep |
|
./pkg/kl10/kl10_calc.F |
61 |
KL10_CALC |
INTEGER I, J, K, Km1, JJ |
|
./pkg/kl10/kl10_calc.F |
79 |
KL10_CALC |
DO J=jMin,jMax |
|
./pkg/kl10/kl10_calc.F |
82 |
KL10_CALC |
CALL FIND_RHO_SCALAR(theta(I,J,K,bi,bj), salt(I,J,K,bi,bj), |
|
./pkg/kl10/kl10_calc.F |
83 |
KL10_CALC |
& totPhiHyd(I,J,K,bi,bj),rhot,myThid ) |
|
./pkg/kl10/kl10_calc.F |
87 |
KL10_CALC |
KLeps(I-1,J-1,1,bi,bj)=0.0 |
|
./pkg/kl10/kl10_calc.F |
92 |
KL10_CALC |
KLviscAr(I,J,1,bi,bj) = viscArNr(1) |
|
./pkg/kl10/kl10_calc.F |
93 |
KL10_CALC |
KLviscold = KLviscAr(I,J,1,bi,bj) ! at previous cell center |
|
./pkg/kl10/kl10_calc.F |
97 |
KL10_CALC |
CALL FIND_RHO_SCALAR(theta(I,J,K,bi,bj), salt(I,J,K,bi |
|
./pkg/kl10/kl10_calc.F |
98 |
KL10_CALC |
& ,bj),totPhiHyd(I,J,K,bi,bj),rhot,myThid ) |
|
./pkg/kpp/kpp_forcing_surf.F |
168 |
KPP_FORCING_SURF |
bo (I,J) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
169 |
KPP_FORCING_SURF |
bosol(I,J) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
172 |
KPP_FORCING_SURF |
boplume(I,J,k) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
174 |
KPP_FORCING_SURF |
boplume(I,J,Nrp1) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
212 |
KPP_FORCING_SURF |
bo(I,J) = - gravity * |
|
./pkg/kpp/kpp_forcing_surf.F |
213 |
KPP_FORCING_SURF |
& ( TTALPHA(I,J,1) * (surfForcT(i,j,bi,bj)+ |
|
./pkg/kpp/kpp_forcing_surf.F |
215 |
KPP_FORCING_SURF |
& SSBETA(I,J,1) * surfForcS(i,j,bi,bj) ) |
|
./pkg/kpp/kpp_forcing_surf.F |
216 |
KPP_FORCING_SURF |
& / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
217 |
KPP_FORCING_SURF |
bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) * |
|
./pkg/kpp/kpp_forcing_surf.F |
219 |
KPP_FORCING_SURF |
& / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
232 |
KPP_FORCING_SURF |
temparray(I,J) = - gravity * |
|
./pkg/kpp/kpp_forcing_surf.F |
233 |
KPP_FORCING_SURF |
& ( SSBETA(I,J,k) * SPforcS(i,j,k) + |
|
./pkg/kpp/kpp_forcing_surf.F |
234 |
KPP_FORCING_SURF |
& TTALPHA(I,J,k)* SPforcT(i,j,k) / HeatCapacity_Cp ) |
|
./pkg/kpp/kpp_forcing_surf.F |
235 |
KPP_FORCING_SURF |
& * recip_rhoConst / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
236 |
KPP_FORCING_SURF |
boplume(I,J,kp1) = boplume(I,J,k)+temparray(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
244 |
KPP_FORCING_SURF |
boplume(I,J,k) = 0. _d 0 |
|
./pkg/kpp/kpp_forcing_surf.F |
246 |
KPP_FORCING_SURF |
boplume(I,J,1) = - gravity * SSBETA(I,J,1) |
|
./pkg/kpp/kpp_forcing_surf.F |
248 |
KPP_FORCING_SURF |
& * recip_rhoConst / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
307 |
KPP_FORCING_SURF |
& maskC(I,J,k,bi,bj) .GT. 0. .AND. |
|
./pkg/kpp/kpp_init_varia.F |
45 |
KPP_INIT_VARIA |
nzmax(i,j,bi,bj) = kLowC(I,J,bi,bj) |
|
./pkg/kpp/kpp_routines.F |
1698 |
STATEKPP |
INTEGER I, J, K |
|
./pkg/kpp/kpp_routines.F |
1731 |
STATEKPP |
DO J = 1-OLy, sNy+OLy |
|
./pkg/kpp/kpp_routines.F |
1733 |
STATEKPP |
RHO1(I,J) = WORK1(I,J) + rhoConst |
|
./pkg/kpp/kpp_routines.F |
1734 |
STATEKPP |
TTALPHA(I,J,1) = WORK2(I,J) |
|
./pkg/kpp/kpp_routines.F |
1735 |
STATEKPP |
SSBETA(I,J,1) = WORK3(I,J) |
|
./pkg/kpp/kpp_routines.F |
1736 |
STATEKPP |
DBSFC(I,J,1) = 0. |
|
./pkg/kpp/kpp_routines.F |
1800 |
STATEKPP |
DO J = 1-OLy, sNy+OLy |
|
./pkg/kpp/kpp_routines.F |
1802 |
STATEKPP |
TTALPHA(I,J,K) = WORK1 (I,J) |
|
./pkg/kpp/kpp_routines.F |
1803 |
STATEKPP |
SSBETA(I,J,K) = WORK2 (I,J) |
|
./pkg/kpp/kpp_routines.F |
1804 |
STATEKPP |
DBLOC(I,J,K-1) = gravity * (RHOK(I,J) - RHOKM1(I,J)) / |
|
./pkg/kpp/kpp_routines.F |
1805 |
STATEKPP |
& (RHOK(I,J) + rhoConst) |
|
./pkg/kpp/kpp_routines.F |
1806 |
STATEKPP |
DBSFC(I,J,K) = gravity * (RHOK(I,J) - RHO1K (I,J)) / |
|
./pkg/kpp/kpp_routines.F |
1807 |
STATEKPP |
& (RHOK(I,J) + rhoConst) |
|
./pkg/kpp/kpp_routines.F |
1814 |
STATEKPP |
DO J = 1-OLy, sNy+OLy |
|
./pkg/kpp/kpp_routines.F |
1816 |
STATEKPP |
TTALPHA(I,J,Nrp1) = TTALPHA(I,J,Nr) |
|
./pkg/kpp/kpp_routines.F |
1817 |
STATEKPP |
SSBETA(I,J,Nrp1) = SSBETA(I,J,Nr) |
|
./pkg/kpp/kpp_routines.F |
1818 |
STATEKPP |
DBLOC(I,J,Nr) = 0. |
|
./pkg/kpp/kpp_routines.F |
1892 |
KPP_DOUBLEDIFF |
INTEGER I, J, K |
|
./pkg/kpp/kpp_routines.F |
1924 |
KPP_DOUBLEDIFF |
DO J = 1-OLy, sNy+OLy |
|
./pkg/kpp/kpp_routines.F |
1926 |
KPP_DOUBLEDIFF |
alphaDT(I,J) = ( theta(I,J,Km1,bi,bj)-theta(I,J,K,bi,bj) ) |
|
./pkg/kpp/kpp_routines.F |
1927 |
KPP_DOUBLEDIFF |
& * 0.5 _d 0 * ABS( TTALPHA(I,J,Km1) + TTALPHA(I,J,K) ) |
|
./pkg/kpp/kpp_routines.F |
1928 |
KPP_DOUBLEDIFF |
betaDS(I,J) = ( salt(I,J,Km1,bi,bj)-salt(I,J,K,bi,bj) ) |
|
./pkg/kpp/kpp_routines.F |
1929 |
KPP_DOUBLEDIFF |
& * 0.5 _d 0 * ( SSBETA(I,J,Km1) + SSBETA(I,J,K) ) |
|
./pkg/kpp/kpp_routines.F |
1930 |
KPP_DOUBLEDIFF |
nuddt(I,J) = 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1931 |
KPP_DOUBLEDIFF |
nudds(I,J) = 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1935 |
KPP_DOUBLEDIFF |
DO J = jMin, jMax |
|
./pkg/kpp/kpp_routines.F |
1941 |
KPP_DOUBLEDIFF |
IF ( alphaDT(I,J) .GT. betaDS(I,J) |
|
./pkg/kpp/kpp_routines.F |
1942 |
KPP_DOUBLEDIFF |
& .AND. betaDS(I,J) .GT. 0. _d 0 ) THEN |
|
./pkg/kpp/kpp_routines.F |
1943 |
KPP_DOUBLEDIFF |
Rrho = MIN( alphaDT(I,J)/betaDS(I,J), Rrho0 ) |
|
./pkg/kpp/kpp_routines.F |
1947 |
KPP_DOUBLEDIFF |
nudds(I,J) = dsfmax * nutmp * nutmp * nutmp |
|
./pkg/kpp/kpp_routines.F |
1949 |
KPP_DOUBLEDIFF |
nuddt(I,J) = 0.7 _d 0 * nudds(I,J) |
|
./pkg/kpp/kpp_routines.F |
1950 |
KPP_DOUBLEDIFF |
ELSEIF ( alphaDT(I,J) .LT. 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1951 |
KPP_DOUBLEDIFF |
& .AND. betaDS(I,J) .LT. 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1952 |
KPP_DOUBLEDIFF |
& .AND.alphaDT(I,J) .GT. betaDS(I,J) ) THEN |
|
./pkg/kpp/kpp_routines.F |
1957 |
KPP_DOUBLEDIFF |
Rrho = alphaDT(I,J)/betaDS(I,J) |
|
./pkg/kpp/kpp_routines.F |
1959 |
KPP_DOUBLEDIFF |
nuddt(I,J) = numol * 0.909 _d 0 |
|
./pkg/kpp/kpp_routines.F |
1966 |
KPP_DOUBLEDIFF |
nudds(I,J) = nuddt(I,J) * MAX( 0.15 _d 0 * Rrho, |
|
./pkg/kpp/kpp_routines.F |
1978 |
KPP_DOUBLEDIFF |
DO J = 1-OLy, sNy+OLy |
|
./pkg/kpp/kpp_routines.F |
1980 |
KPP_DOUBLEDIFF |
kappaRT(I,J,K) = kappaRT(I,J,K) + nuddt(I,J) |
|
./pkg/kpp/kpp_routines.F |
1981 |
KPP_DOUBLEDIFF |
kappaRS(I,J,K) = kappaRS(I,J,K) + nudds(I,J) |
|
./pkg/kpp/kpp_transport_s.F |
63 |
KPP_TRANSPORT_S |
INTEGER I, J |
|
./pkg/land/land_ini_vars.F |
63 |
LAND_INI_VARS |
DO J=1-Oly,sNy+Oly |
|
./pkg/land/land_ini_vars.F |
73 |
LAND_INI_VARS |
DO J=1-Oly,sNy+Oly |
|
./pkg/mdsio/mdsio_gl.F |
1150 |
MDSREADFIELD_2D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
1152 |
MDSREADFIELD_2D_GL |
global(I,J) = xy_buffer_r4(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
1160 |
MDSREADFIELD_2D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
1162 |
MDSREADFIELD_2D_GL |
global(I,J) = xy_buffer_r8(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
1176 |
MDSREADFIELD_2D_GL |
DO J=1,sNy |
|
./pkg/mdsio/mdsio_gl.F |
1177 |
MDSREADFIELD_2D_GL |
JJ=((jp-1)*nSy+(bj-1))*sNy+J |
|
./pkg/mdsio/mdsio_gl.F |
1381 |
MDSWRITEFIELD_2D_GL |
DO J=1,sNy |
|
./pkg/mdsio/mdsio_gl.F |
1382 |
MDSWRITEFIELD_2D_GL |
JJ=((jp-1)*nSy+(bj-1))*sNy+J |
|
./pkg/mdsio/mdsio_gl.F |
1396 |
MDSWRITEFIELD_2D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
1398 |
MDSWRITEFIELD_2D_GL |
xy_buffer_r4(I,J) = global(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
1406 |
MDSWRITEFIELD_2D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
1408 |
MDSWRITEFIELD_2D_GL |
xy_buffer_r8(I,J) = global(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
363 |
MDSREADFIELD_3D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
365 |
MDSREADFIELD_3D_GL |
global(I,J) = xy_buffer_r4(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
373 |
MDSREADFIELD_3D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
375 |
MDSREADFIELD_3D_GL |
global(I,J) = xy_buffer_r8(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
389 |
MDSREADFIELD_3D_GL |
DO J=1,sNy |
|
./pkg/mdsio/mdsio_gl.F |
390 |
MDSREADFIELD_3D_GL |
JJ=((jp-1)*nSy+(bj-1))*sNy+J |
|
./pkg/mdsio/mdsio_gl.F |
592 |
MDSWRITEFIELD_3D_GL |
DO J=1,sNy |
|
./pkg/mdsio/mdsio_gl.F |
593 |
MDSWRITEFIELD_3D_GL |
JJ=((jp-1)*nSy+(bj-1))*sNy+J |
|
./pkg/mdsio/mdsio_gl.F |
607 |
MDSWRITEFIELD_3D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
609 |
MDSWRITEFIELD_3D_GL |
xy_buffer_r4(I,J) = global(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
617 |
MDSWRITEFIELD_3D_GL |
DO J=1,Ny |
|
./pkg/mdsio/mdsio_gl.F |
619 |
MDSWRITEFIELD_3D_GL |
xy_buffer_r8(I,J) = global(I,J) |
|
./pkg/mom_common/mom_calc_absvort3.F |
46 |
MOM_CALC_ABSVORT3 |
omega3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
114 |
MOM_CALC_RELVORT3 |
J=1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
117 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
118 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
119 |
MOM_CALC_RELVORT3 |
& ( vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
120 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
121 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
131 |
MOM_CALC_RELVORT3 |
vort3(I-1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
132 |
MOM_CALC_RELVORT3 |
& recip_rAz(I-1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
133 |
MOM_CALC_RELVORT3 |
& vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
134 |
MOM_CALC_RELVORT3 |
& -vFld(I-2,J)*dyC(I-2,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
135 |
MOM_CALC_RELVORT3 |
& -uFld(I-1,J)*dxC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
136 |
MOM_CALC_RELVORT3 |
& +vFld(I+0,J-1)*dyC(I+0,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
140 |
MOM_CALC_RELVORT3 |
vort3(I,J-1)=vort3(I-1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
155 |
MOM_CALC_RELVORT3 |
J=1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
159 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
160 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
161 |
MOM_CALC_RELVORT3 |
& (-uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
162 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
163 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
166 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
167 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
168 |
MOM_CALC_RELVORT3 |
& (-vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
169 |
MOM_CALC_RELVORT3 |
& +uFld(I,J-1)*dxC(I,J-1,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
170 |
MOM_CALC_RELVORT3 |
& - uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
173 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
174 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
175 |
MOM_CALC_RELVORT3 |
& (+uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
176 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
177 |
MOM_CALC_RELVORT3 |
& - vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
188 |
MOM_CALC_RELVORT3 |
vort3(I+1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
189 |
MOM_CALC_RELVORT3 |
& recip_rAz(I+1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
190 |
MOM_CALC_RELVORT3 |
& vFld(I+1,J)*dyC(I+1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
191 |
MOM_CALC_RELVORT3 |
& -vFld(I-0,J)*dyC(I-0,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
192 |
MOM_CALC_RELVORT3 |
& -uFld(I+1,J)*dxC(I+1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
193 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J-1)*dyC(I-1,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
197 |
MOM_CALC_RELVORT3 |
vort3(I,J-1)=vort3(I+1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
212 |
MOM_CALC_RELVORT3 |
J=sNy+1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
216 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
217 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
218 |
MOM_CALC_RELVORT3 |
& (+uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
219 |
MOM_CALC_RELVORT3 |
& +vFld(I,J)*dyC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
220 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
223 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
224 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
225 |
MOM_CALC_RELVORT3 |
& (-uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
226 |
MOM_CALC_RELVORT3 |
& +uFld(I,J-1)*dxC(I,J-1,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
227 |
MOM_CALC_RELVORT3 |
& + vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
230 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
231 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
232 |
MOM_CALC_RELVORT3 |
& (+vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
233 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
234 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
245 |
MOM_CALC_RELVORT3 |
vort3(I-1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
246 |
MOM_CALC_RELVORT3 |
& recip_rAz(I-1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
247 |
MOM_CALC_RELVORT3 |
& vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
248 |
MOM_CALC_RELVORT3 |
& -vFld(I-2,J)*dyC(I-2,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
249 |
MOM_CALC_RELVORT3 |
& +vFld(I-0,J+1)*dyC(I-0,J+1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
250 |
MOM_CALC_RELVORT3 |
& +uFld(I-1,J-1)*dxC(I-1,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
254 |
MOM_CALC_RELVORT3 |
vort3(I,J+1)=vort3(I-1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
269 |
MOM_CALC_RELVORT3 |
J=sNy+1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
273 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
274 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
275 |
MOM_CALC_RELVORT3 |
& (-uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
276 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
277 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
280 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
281 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
282 |
MOM_CALC_RELVORT3 |
& (+uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
283 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
284 |
MOM_CALC_RELVORT3 |
& - vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
295 |
MOM_CALC_RELVORT3 |
vort3(I+1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
296 |
MOM_CALC_RELVORT3 |
& recip_rAz(I+1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
297 |
MOM_CALC_RELVORT3 |
& vFld(I+1,J)*dyC(I+1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
298 |
MOM_CALC_RELVORT3 |
& -vFld(I-0,J)*dyC(I-0,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
299 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J+1)*dyC(I-1,J+1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
300 |
MOM_CALC_RELVORT3 |
& +uFld(I+1,J-1)*dxC(I+1,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
304 |
MOM_CALC_RELVORT3 |
vort3(I,J+1)=vort3(I+1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
46 |
MOM_CALC_RELVORT3 |
DO J=1-OLy,sNy+OLy |
|
./pkg/mom_common/mom_calc_relvort3.F |
48 |
MOM_CALC_RELVORT3 |
vort3(I,J) = 0. _d 0 |
|
./pkg/mom_common/mom_calc_relvort3.F |
53 |
MOM_CALC_RELVORT3 |
DO J=2-OLy,sNy+OLy |
|
./pkg/mom_common/mom_calc_relvort3.F |
57 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
58 |
MOM_CALC_RELVORT3 |
& recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
59 |
MOM_CALC_RELVORT3 |
& ( vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
60 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
61 |
MOM_CALC_RELVORT3 |
& -( uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
62 |
MOM_CALC_RELVORT3 |
& -uFld(I,J-1)*dxC(I,J-1,bi,bj) ) |
|
./pkg/mom_fluxform/mom_u_del2u.F |
46 |
MOM_U_DEL2U |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_u_del2u.F |
60 |
MOM_U_DEL2U |
& *sqCosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_del2u.F |
78 |
MOM_U_DEL2U |
& *sqCosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_xviscflux.F |
50 |
MOM_U_XVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_u_xviscflux.F |
60 |
MOM_U_XVISCFLUX |
& *cosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_xviscflux.F |
63 |
MOM_U_XVISCFLUX |
& *sqCosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_xviscflux.F |
65 |
MOM_U_XVISCFLUX |
& *cosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_yviscflux.F |
51 |
MOM_U_YVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_u_yviscflux.F |
62 |
MOM_U_YVISCFLUX |
& *cosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_yviscflux.F |
67 |
MOM_U_YVISCFLUX |
& *sqCosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_u_yviscflux.F |
69 |
MOM_U_YVISCFLUX |
& *cosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_del2v.F |
45 |
MOM_V_DEL2V |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_v_del2v.F |
59 |
MOM_V_DEL2V |
& *sqCosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_del2v.F |
74 |
MOM_V_DEL2V |
& *sqCosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_xviscflux.F |
51 |
MOM_V_XVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_v_xviscflux.F |
61 |
MOM_V_XVISCFLUX |
& *cosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_xviscflux.F |
64 |
MOM_V_XVISCFLUX |
& *sqCosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_xviscflux.F |
66 |
MOM_V_XVISCFLUX |
& *cosFacV(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_yviscflux.F |
50 |
MOM_V_YVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_v_yviscflux.F |
61 |
MOM_V_YVISCFLUX |
& *cosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_yviscflux.F |
66 |
MOM_V_YVISCFLUX |
& *sqCosFacU(J,bi,bj) |
|
./pkg/mom_fluxform/mom_v_yviscflux.F |
68 |
MOM_V_YVISCFLUX |
& *cosFacU(J,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
100 |
MOM_VI_CORIOLIS |
& )*uBarXY*recip_dyC(i,j,bi,bj)*_maskS(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
34 |
MOM_VI_CORIOLIS |
INTEGER I,J |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
41 |
MOM_VI_CORIOLIS |
DO J=1-Oly,sNy+Oly-1 |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
52 |
MOM_VI_CORIOLIS |
& )*vBarXY*recip_dxC(i,j,bi,bj)*_maskW(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
57 |
MOM_VI_CORIOLIS |
DO J=1-Oly,sNy+Oly-1 |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
67 |
MOM_VI_CORIOLIS |
& )*vBarXY*recip_dxC(i,j,bi,bj)*_maskW(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
74 |
MOM_VI_CORIOLIS |
DO J=2-Oly,sNy+Oly |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
85 |
MOM_VI_CORIOLIS |
& )*uBarXY*recip_dyC(i,j,bi,bj)*_maskS(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
90 |
MOM_VI_CORIOLIS |
DO J=2-Oly,sNy+Oly |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
103 |
MOM_VI_U_VERTSHEAR |
uZp=(mask_Kp1*uFld(I,J,Kp1,bi,bj)-uFld(I,J,K,bi,bj))*rkSign |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
123 |
MOM_VI_U_VERTSHEAR |
uShearTerm(I,J)=-0.5* |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
130 |
MOM_VI_U_VERTSHEAR |
uShearTerm(I,J)=-0.5*(wBarXp*uZp+wBarXm*uZm) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
31 |
MOM_VI_U_VERTSHEAR |
INTEGER I,J,Kp1,Km1 |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
56 |
MOM_VI_U_VERTSHEAR |
DO J=1-OLy,sNy+OLy |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
69 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(I,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
70 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,K,bi,bj)*rA(i-1,j,bi,bj)*maskC(I-1,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
76 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj)*rA(i,j,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
77 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,Kp1,bi,bj)*rA(i-1,j,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
83 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,K,bi,bj)*maskC(I,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
84 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,K,bi,bj)*maskC(I-1,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
89 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
90 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
96 |
MOM_VI_U_VERTSHEAR |
uZm=(uFld(I,J,K,bi,bj)-mask_Km1*uFld(I,J,Km1,bi,bj))*rkSign |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
103 |
MOM_VI_V_VERTSHEAR |
vZp=(mask_Kp1*vFld(I,J,Kp1,bi,bj)-vFld(I,J,K,bi,bj))*rkSign |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
123 |
MOM_VI_V_VERTSHEAR |
vShearTerm(I,J)=-0.5* |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
130 |
MOM_VI_V_VERTSHEAR |
vShearTerm(I,J)=-0.5*(wBarYp*vZp+wBarYm*vZm) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
31 |
MOM_VI_V_VERTSHEAR |
INTEGER I,J,Kp1,Km1 |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
56 |
MOM_VI_V_VERTSHEAR |
DO J=2-OLy,sNy+OLy |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
69 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(i,j,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
70 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,K,bi,bj)*rA(i,j-1,bi,bj)*maskC(i,j-1,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
76 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj)*rA(i,j,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
77 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,Kp1,bi,bj)*rA(i,j-1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
83 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,K,bi,bj)*maskC(i,j,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
84 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,K,bi,bj)*maskC(i,j-1,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
89 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
90 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
96 |
MOM_VI_V_VERTSHEAR |
vZm=(vFld(I,J,K,bi,bj)-mask_Km1*vFld(I,J,Km1,bi,bj))*rkSign |
|
./pkg/monitor/mon_advcfl.F |
33 |
MON_ADVCFL |
INTEGER bi,bj,I,J,K |
|
./pkg/monitor/mon_advcfl.F |
41 |
MON_ADVCFL |
DO J=1,sNy |
|
./pkg/monitor/mon_advcfl.F |
43 |
MON_ADVCFL |
tmpVal=abs(U(I,J,K,bi,bj))*rDx(I,J,BI,BJ)*dT |
|
./pkg/monitor/mon_advcflw.F |
33 |
MON_ADVCFLW |
INTEGER bi,bj,I,J,K |
|
./pkg/monitor/mon_advcflw.F |
41 |
MON_ADVCFLW |
DO J=1,sNy |
|
./pkg/monitor/mon_advcflw.F |
43 |
MON_ADVCFLW |
tmpVal=abs(W(I,J,K,bi,bj))*rDz(K)*dT |
|
./pkg/monitor/mon_advcflw2.F |
35 |
MON_ADVCFLW2 |
INTEGER bi,bj,I,J,K |
|
./pkg/monitor/mon_advcflw2.F |
43 |
MON_ADVCFLW2 |
DO J=1,sNy |
|
./pkg/monitor/mon_advcflw2.F |
45 |
MON_ADVCFLW2 |
tmpVal=abs(W(I,J,K,bi,bj))*dT |
|
./pkg/monitor/mon_advcflw2.F |
46 |
MON_ADVCFLW2 |
& *max( rDrF(K)*rHfac(I,J,K,bi,bj), |
|
./pkg/monitor/mon_advcflw2.F |
47 |
MON_ADVCFLW2 |
& rDrF(K-1)*rHfac(I,J,K-1,bi,bj) ) |
|
./pkg/monitor/mon_stats_rl.F |
134 |
MON_STATS_RL |
DO J=1,sNy |
|
./pkg/monitor/mon_stats_rl.F |
136 |
MON_STATS_RL |
tmpVal=arr(I,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
137 |
MON_STATS_RL |
IF (arrMask(I,J,K,bi,bj).NE.0.) THEN |
|
./pkg/monitor/mon_stats_rl.F |
138 |
MON_STATS_RL |
tmpVol=arrArea(I,J,bi,bj)*arrhFac(I,J,K,bi,bj)*arrDr(K) |
|
./pkg/monitor/mon_stats_rl.F |
139 |
MON_STATS_RL |
& *arrMask(I,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
35 |
MON_STATS_RL |
INTEGER bi,bj,I,J,K |
|
./pkg/monitor/mon_stats_rl.F |
68 |
MON_STATS_RL |
DO J=1,sNy |
|
./pkg/monitor/mon_stats_rl.F |
70 |
MON_STATS_RL |
tmpVal=arr(I,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
71 |
MON_STATS_RL |
IF (arrMask(I,J,K,bi,bj).NE.0. .AND. noPnts) THEN |
|
./pkg/monitor/mon_stats_rl.F |
76 |
MON_STATS_RL |
IF (arrMask(I,J,K,bi,bj).NE.0.) THEN |
|
./pkg/monitor/mon_stats_rl.F |
81 |
MON_STATS_RL |
& (arr(I+1,J,K,bi,bj)-tmpVal)*arrMask(I+1,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
82 |
MON_STATS_RL |
& +(arr(I-1,J,K,bi,bj)-tmpVal)*arrMask(I-1,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
83 |
MON_STATS_RL |
& +(arr(I,J+1,K,bi,bj)-tmpVal)*arrMask(I,J+1,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
84 |
MON_STATS_RL |
& +(arr(I,J-1,K,bi,bj)-tmpVal)*arrMask(I,J-1,K,bi,bj) |
|
./pkg/monitor/mon_stats_rl.F |
87 |
MON_STATS_RL |
tmpVol = arrArea(I,J,bi,bj)*arrhFac(I,J,K,bi,bj)*arrDr(K) |
|
./pkg/monitor/mon_stats_rl.F |
88 |
MON_STATS_RL |
& *arrMask(I,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rs.F |
101 |
MON_STATS_RS |
DO J=1,sNy |
|
./pkg/monitor/mon_stats_rs.F |
103 |
MON_STATS_RS |
tmpVal=arr(I,J,K,bi,bj) |
|
./pkg/monitor/mon_stats_rs.F |
32 |
MON_STATS_RS |
INTEGER bi,bj,I,J,K |
|
./pkg/monitor/mon_stats_rs.F |
54 |
MON_STATS_RS |
DO J=1,sNy |
|
./pkg/monitor/mon_stats_rs.F |
56 |
MON_STATS_RS |
tmpVal=arr(I,J,K,bi,bj) |
|
./pkg/my82/my82_calc.F |
101 |
MY82_CALC |
DO J=jMin,jMax |
|
./pkg/my82/my82_calc.F |
103 |
MY82_CALC |
RiTmp = MIN(RiNumber(I,J),RiMax) |
|
./pkg/my82/my82_calc.F |
111 |
MY82_CALC |
SH(I,J,K) = SHtmp |
|
./pkg/my82/my82_calc.F |
112 |
MY82_CALC |
SM(I,J,K) = SHtmp*(beta1-beta2*RiFlux)/(beta3-beta4*RiFlux) |
|
./pkg/my82/my82_calc.F |
115 |
MY82_CALC |
& b1*(SH(I,J,K)*GH(I,J) + SM(I,J,K)*GM(I,J)) ) |
|
./pkg/my82/my82_calc.F |
116 |
MY82_CALC |
tke(I,J,K) = sqrt(tkesquare) |
|
./pkg/my82/my82_calc.F |
128 |
MY82_CALC |
DO J=jMin,jMax |
|
./pkg/my82/my82_calc.F |
130 |
MY82_CALC |
GH(I,J) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
131 |
MY82_CALC |
GM(I,J) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
139 |
MY82_CALC |
DO J=jMin,jMax |
|
./pkg/my82/my82_calc.F |
141 |
MY82_CALC |
GM(I,J) = GM(I,J) + tke(I,J,K)*rF(K) |
|
./pkg/my82/my82_calc.F |
142 |
MY82_CALC |
GH(I,J) = GH(I,J) + tke(I,J,K) |
|
./pkg/my82/my82_calc.F |
148 |
MY82_CALC |
DO J=jMin,jMax |
|
./pkg/my82/my82_calc.F |
150 |
MY82_CALC |
IF ( GH(I,J) .EQ. 0. _d 0 ) THEN |
|
./pkg/my82/my82_calc.F |
151 |
MY82_CALC |
MYhbl(I,J,bi,bj) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
153 |
MY82_CALC |
MYhbl(I,J,bi,bj) = -GM(I,J)/GH(I,J)*MYhblScale |
|
./pkg/my82/my82_calc.F |
160 |
MY82_CALC |
DO J=jMin,jMax |
|
./pkg/my82/my82_calc.F |
162 |
MY82_CALC |
tkel = MYhbl(I,J,bi,bj)*tke(I,J,K) |
|
./pkg/my82/my82_calc.F |
164 |
MY82_CALC |
MYviscAr(I,J,K,bi,bj) = MYhbl(I,J,bi,bj)*tkel*SM(I,J,K) |
|
./pkg/my82/my82_calc.F |
165 |
MY82_CALC |
MYdiffKr(I,J,K,bi,bj) = MYhbl(I,J,bi,bj)*tkel*SH(I,J,K) |
|
./pkg/my82/my82_calc.F |
167 |
MY82_CALC |
MYviscAr(I,J,K,bi,bj) = MAX(MYviscAr(I,J,K,bi,bj), |
|
./pkg/my82/my82_calc.F |
169 |
MY82_CALC |
MYdiffKr(I,J,K,bi,bj) = MAX(MYdiffKr(I,J,K,bi,bj), |
|
./pkg/my82/my82_calc.F |
176 |
MY82_CALC |
MYviscAr(I,J,K,bi,bj) = MIN(MYviscAr(I,J,K,bi,bj),MYviscMax) |
|
./pkg/my82/my82_calc.F |
177 |
MY82_CALC |
& * maskC(I,J,K,bi,bj) |
|
./pkg/my82/my82_calc.F |
178 |
MY82_CALC |
MYdiffKr(I,J,K,bi,bj) = MIN(MYdiffKr(I,J,K,bi,bj),MYdiffMax) |
|
./pkg/my82/my82_calc.F |
179 |
MY82_CALC |
& * maskC(I,J,K,bi,bj) |
|
./pkg/my82/my82_calc.F |
62 |
MY82_CALC |
INTEGER I, J, K |
|
./pkg/my82/my82_calc.F |
79 |
MY82_CALC |
DO J=1-OLy,sNy+OLy |
|
./pkg/my82/my82_calc.F |
81 |
MY82_CALC |
GH(I,J) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
82 |
MY82_CALC |
GM(I,J) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
86 |
MY82_CALC |
DO J=1-OLy,sNy+OLy |
|
./pkg/my82/my82_calc.F |
88 |
MY82_CALC |
SH(I,J,K) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
89 |
MY82_CALC |
SM(I,J,K) = 0. _d 0 |
|
./pkg/my82/my82_calc.F |
90 |
MY82_CALC |
tke(I,J,K) = 0. _d 0 |
|
./pkg/my82/my82_ri_number.F |
102 |
MY82_RI_NUMBER |
DO J= jMin, jMax |
|
./pkg/my82/my82_ri_number.F |
104 |
MY82_RI_NUMBER |
tempu= .5 _d 0*( uVel(I,J,Km1,bi,bj)+uVel(I+1,J,Km1,bi,bj) |
|
./pkg/my82/my82_ri_number.F |
105 |
MY82_RI_NUMBER |
& - (uVel(I,J,K ,bi,bj)+uVel(I+1,J,K ,bi,bj)) ) |
|
./pkg/my82/my82_ri_number.F |
107 |
MY82_RI_NUMBER |
tempv= .5 _d 0*( vVel(I,J,Km1,bi,bj)+vVel(I,J+1,Km1,bi,bj) |
|
./pkg/my82/my82_ri_number.F |
108 |
MY82_RI_NUMBER |
& - (vVel(I,J,K ,bi,bj)+vVel(I,J+1,K ,bi,bj)) ) |
|
./pkg/my82/my82_ri_number.F |
110 |
MY82_RI_NUMBER |
vertShear(I,J) = tempu*tempu+tempv*tempv |
|
./pkg/my82/my82_ri_number.F |
115 |
MY82_RI_NUMBER |
buoyFreq(I,J) = gravity*mass2rUnit * |
|
./pkg/my82/my82_ri_number.F |
116 |
MY82_RI_NUMBER |
& (rhoKm1(I,J) - rhoK(I,J))*recip_drC(K) |
|
./pkg/my82/my82_ri_number.F |
121 |
MY82_RI_NUMBER |
RiNumber(I,J) = -buoyFreq(I,J)/max(vertShear(I,J),epsilon) |
|
./pkg/my82/my82_ri_number.F |
128 |
MY82_RI_NUMBER |
DO J=jMin,jMax |
|
./pkg/my82/my82_ri_number.F |
130 |
MY82_RI_NUMBER |
RiTmp(I,J) = RiNumber(I,J) |
|
./pkg/my82/my82_ri_number.F |
133 |
MY82_RI_NUMBER |
DO J=1-Oly+1,sNy+Oly-1 |
|
./pkg/my82/my82_ri_number.F |
135 |
MY82_RI_NUMBER |
RiNumber(I,J) = p5*RiTmp(I,J) |
|
./pkg/my82/my82_ri_number.F |
136 |
MY82_RI_NUMBER |
& + p125*RiTmp(I-1,J) + p125*RiTmp(I+1,J) |
|
./pkg/my82/my82_ri_number.F |
137 |
MY82_RI_NUMBER |
& + p125*RiTmp(I,J-1) + p125*RiTmp(I,J+1) |
|
./pkg/my82/my82_ri_number.F |
60 |
MY82_RI_NUMBER |
INTEGER I,J,Km1 |
|
./pkg/obcs/obcs_adjust_uvice.F |
146 |
OBCS_ADJUST_UVICE |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_adjust_uvice.F |
148 |
OBCS_ADJUST_UVICE |
IF (OB_Ie(J,bi,bj).NE.OB_indexNone) THEN |
|
./pkg/obcs/obcs_adjust_uvice.F |
150 |
OBCS_ADJUST_UVICE |
uFld(OB_Ie(J,bi,bj),J,bi,bj) = |
|
./pkg/obcs/obcs_adjust_uvice.F |
151 |
OBCS_ADJUST_UVICE |
& _maskW(OB_Ie(J,bi,bj),J,K,bi,bj) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
152 |
OBCS_ADJUST_UVICE |
& uFld(OB_Ie(J,bi,bj)-1,J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
153 |
OBCS_ADJUST_UVICE |
vFld(OB_Ie(J,bi,bj),J,bi,bj) = |
|
./pkg/obcs/obcs_adjust_uvice.F |
154 |
OBCS_ADJUST_UVICE |
& _maskS(OB_Ie(J,bi,bj),J,K,bi,bj) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
155 |
OBCS_ADJUST_UVICE |
& vFld(OB_Ie(J,bi,bj)-1,J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
158 |
OBCS_ADJUST_UVICE |
uFld(OB_Ie(J,bi,bj),J,bi,bj)= |
|
./pkg/obcs/obcs_adjust_uvice.F |
159 |
OBCS_ADJUST_UVICE |
& max(OBEuice(J,bi,bj),uFld(OB_Ie(J,bi,bj),J,bi,bj)) |
|
./pkg/obcs/obcs_adjust_uvice.F |
160 |
OBCS_ADJUST_UVICE |
& * _maskW(OB_Ie(J,bi,bj),J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
162 |
OBCS_ADJUST_UVICE |
uFld(OB_Ie(J,bi,bj),J,bi,bj)=OBEuice(J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
163 |
OBCS_ADJUST_UVICE |
& * _maskW(OB_Ie(J,bi,bj),J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
165 |
OBCS_ADJUST_UVICE |
vFld(OB_Ie(J,bi,bj),J,bi,bj)=OBEvice(J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
166 |
OBCS_ADJUST_UVICE |
& * _maskS(OB_Ie(J,bi,bj),J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
168 |
OBCS_ADJUST_UVICE |
if ( _maskW(OB_Ie(J,bi,bj) ,J,K,bi,bj) .NE. 0. .AND. |
|
./pkg/obcs/obcs_adjust_uvice.F |
169 |
OBCS_ADJUST_UVICE |
& _maskW(OB_Ie(J,bi,bj)-2,J,K,bi,bj) .NE. 0. ) |
|
./pkg/obcs/obcs_adjust_uvice.F |
170 |
OBCS_ADJUST_UVICE |
& uFld(OB_Ie(J,bi,bj)-1,J,bi,bj) = 0.5 _d 0 * |
|
./pkg/obcs/obcs_adjust_uvice.F |
171 |
OBCS_ADJUST_UVICE |
& ( uFld(OB_Ie(J,bi,bj) ,J,bi,bj) + |
|
./pkg/obcs/obcs_adjust_uvice.F |
172 |
OBCS_ADJUST_UVICE |
& uFld(OB_Ie(J,bi,bj)-2,J,bi,bj) ) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
173 |
OBCS_ADJUST_UVICE |
& _maskW(OB_Ie(J,bi,bj)-1,J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
176 |
OBCS_ADJUST_UVICE |
if ( _maskS(OB_Ie(J,bi,bj) ,J,K,bi,bj) .NE. 0. .AND. |
|
./pkg/obcs/obcs_adjust_uvice.F |
177 |
OBCS_ADJUST_UVICE |
& _maskS(OB_Ie(J,bi,bj)-2,J,K,bi,bj) .NE. 0. ) |
|
./pkg/obcs/obcs_adjust_uvice.F |
178 |
OBCS_ADJUST_UVICE |
& vFld(OB_Ie(J,bi,bj)-1,J,bi,bj) = 0.5 _d 0 * |
|
./pkg/obcs/obcs_adjust_uvice.F |
179 |
OBCS_ADJUST_UVICE |
& ( vFld(OB_Ie(J,bi,bj) ,J,bi,bj) + |
|
./pkg/obcs/obcs_adjust_uvice.F |
180 |
OBCS_ADJUST_UVICE |
& vFld(OB_Ie(J,bi,bj)-2,J,bi,bj) ) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
181 |
OBCS_ADJUST_UVICE |
& _maskS(OB_Ie(J,bi,bj)-1,J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
191 |
OBCS_ADJUST_UVICE |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_adjust_uvice.F |
193 |
OBCS_ADJUST_UVICE |
IF (OB_Iw(J,bi,bj).NE.OB_indexNone) THEN |
|
./pkg/obcs/obcs_adjust_uvice.F |
195 |
OBCS_ADJUST_UVICE |
uFld(OB_Iw(J,bi,bj)+1,J,bi,bj)= |
|
./pkg/obcs/obcs_adjust_uvice.F |
196 |
OBCS_ADJUST_UVICE |
& _maskW(OB_Iw(J,bi,bj)+1,J,K,bi,bj) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
197 |
OBCS_ADJUST_UVICE |
& uFld(OB_Iw(J,bi,bj)+2,J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
203 |
OBCS_ADJUST_UVICE |
vFld(OB_Iw(J,bi,bj)+1,J,bi,bj)= |
|
./pkg/obcs/obcs_adjust_uvice.F |
204 |
OBCS_ADJUST_UVICE |
& _maskS(OB_Iw(J,bi,bj)+1,J,K,bi,bj) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
205 |
OBCS_ADJUST_UVICE |
& vFld(OB_Iw(J,bi,bj)+2,J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
208 |
OBCS_ADJUST_UVICE |
uFld(OB_Iw(J,bi,bj)+1,J,bi,bj)= |
|
./pkg/obcs/obcs_adjust_uvice.F |
209 |
OBCS_ADJUST_UVICE |
& min(OBWuice(J,bi,bj),uFld(OB_Iw(J,bi,bj)+1,J,bi,bj)) |
|
./pkg/obcs/obcs_adjust_uvice.F |
210 |
OBCS_ADJUST_UVICE |
& * _maskW(OB_Iw(J,bi,bj)+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
212 |
OBCS_ADJUST_UVICE |
uFld(OB_Iw(J,bi,bj)+1,J,bi,bj)=OBWuice(J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
213 |
OBCS_ADJUST_UVICE |
& * _maskW(OB_Iw(J,bi,bj)+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
215 |
OBCS_ADJUST_UVICE |
vFld(OB_Iw(J,bi,bj),J,bi,bj)=OBWvice(J,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
216 |
OBCS_ADJUST_UVICE |
& * _maskS(OB_Iw(J,bi,bj),J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
218 |
OBCS_ADJUST_UVICE |
if ( _maskW(OB_Iw(J,bi,bj)+1,J,K,bi,bj) .NE. 0. .AND. |
|
./pkg/obcs/obcs_adjust_uvice.F |
219 |
OBCS_ADJUST_UVICE |
& _maskW(OB_Iw(J,bi,bj)+3,J,K,bi,bj) .NE. 0. ) |
|
./pkg/obcs/obcs_adjust_uvice.F |
220 |
OBCS_ADJUST_UVICE |
& uFld(OB_Iw(J,bi,bj)+2,J,bi,bj) = 0.5 _d 0 * |
|
./pkg/obcs/obcs_adjust_uvice.F |
221 |
OBCS_ADJUST_UVICE |
& ( uFld(OB_Iw(J,bi,bj)+1,J,bi,bj) + |
|
./pkg/obcs/obcs_adjust_uvice.F |
222 |
OBCS_ADJUST_UVICE |
& uFld(OB_Iw(J,bi,bj)+3,J,bi,bj) ) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
223 |
OBCS_ADJUST_UVICE |
& _maskW(OB_Iw(J,bi,bj)+2,J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
226 |
OBCS_ADJUST_UVICE |
if ( _maskS(OB_Iw(J,bi,bj) ,J,K,bi,bj) .NE. 0. .AND. |
|
./pkg/obcs/obcs_adjust_uvice.F |
227 |
OBCS_ADJUST_UVICE |
& _maskS(OB_Iw(J,bi,bj)+2,J,K,bi,bj) .NE. 0. ) |
|
./pkg/obcs/obcs_adjust_uvice.F |
228 |
OBCS_ADJUST_UVICE |
& vFld(OB_Iw(J,bi,bj)+1,J,bi,bj) = 0.5 _d 0 * |
|
./pkg/obcs/obcs_adjust_uvice.F |
229 |
OBCS_ADJUST_UVICE |
& ( vFld(OB_Iw(J,bi,bj) ,J,bi,bj) + |
|
./pkg/obcs/obcs_adjust_uvice.F |
230 |
OBCS_ADJUST_UVICE |
& vFld(OB_Iw(J,bi,bj)+2,J,bi,bj) ) * |
|
./pkg/obcs/obcs_adjust_uvice.F |
231 |
OBCS_ADJUST_UVICE |
& _maskS(OB_Iw(J,bi,bj)+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_adjust_uvice.F |
40 |
OBCS_ADJUST_UVICE |
INTEGER I,J,K,bi,bj |
|
./pkg/obcs/obcs_apply_seaice.F |
212 |
OBCS_APPLY_SEAICE |
Iobc = OB_Ie(J,bi,bj) |
|
./pkg/obcs/obcs_apply_seaice.F |
95 |
OBCS_APPLY_SEAICE |
Iobc = OB_Ie(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
110 |
OBCS_CALC_STEVENS |
INTEGER I, J, K |
|
./pkg/obcs/obcs_calc_stevens.F |
174 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
175 |
OBCS_CALC_STEVENS |
I = OB_Ie(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
185 |
OBCS_CALC_STEVENS |
uMerBar = uMerBar + uVel(I-1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
187 |
OBCS_CALC_STEVENS |
uMerBar = uMerBar + OBEuStevens(J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
189 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacW(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
190 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacW(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
195 |
OBCS_CALC_STEVENS |
uMerPri(K) = (uVel(I-1,J,K,bi,bj)-uMerBar) |
|
./pkg/obcs/obcs_calc_stevens.F |
197 |
OBCS_CALC_STEVENS |
uMerPri(K) = (OBEuStevens(J,K,bi,bj)-uMerBar) |
|
./pkg/obcs/obcs_calc_stevens.F |
199 |
OBCS_CALC_STEVENS |
& * _maskW(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
205 |
OBCS_CALC_STEVENS |
uMerBar = uMerBar + OBEu(J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
206 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacW(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
207 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacW(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
213 |
OBCS_CALC_STEVENS |
OBEu(J,K,bi,bj) = (uMerBar + uMerPri(K)) |
|
./pkg/obcs/obcs_calc_stevens.F |
214 |
OBCS_CALC_STEVENS |
& * _maskW(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
217 |
OBCS_CALC_STEVENS |
OBEw(J,K,bi,bj)=0. |
|
./pkg/obcs/obcs_calc_stevens.F |
225 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
226 |
OBCS_CALC_STEVENS |
I = OB_Ie(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
228 |
OBCS_CALC_STEVENS |
OBEeta(J,bi,bj) = etaN(I-1,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
236 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
237 |
OBCS_CALC_STEVENS |
I = OB_Ie(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
239 |
OBCS_CALC_STEVENS |
cflMer(J,K) = 0.5 _d 0 * _dxC(I-1,J,bi,bj)/dTtracerLev(K) |
|
./pkg/obcs/obcs_calc_stevens.F |
242 |
OBCS_CALC_STEVENS |
gFacM(J,K) = ABS(MIN(SIGN(1.D0,uVel(I,J,K,bi,bj)),0.D0)) |
|
./pkg/obcs/obcs_calc_stevens.F |
244 |
OBCS_CALC_STEVENS |
cflMer(J,K) = 0. _d 0 |
|
./pkg/obcs/obcs_calc_stevens.F |
245 |
OBCS_CALC_STEVENS |
gFacM (J,K) = 0. _d 0 |
|
./pkg/obcs/obcs_calc_stevens.F |
301 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
302 |
OBCS_CALC_STEVENS |
I = OB_Iw(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
312 |
OBCS_CALC_STEVENS |
uMerBar = uMerBar + uVel(I+2,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
314 |
OBCS_CALC_STEVENS |
uMerBar = uMerBar + OBWuStevens(J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
316 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacW(I+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
317 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacW(I+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
322 |
OBCS_CALC_STEVENS |
uMerPri(K) = (uVel(I+2,J,K,bi,bj)-uMerBar) |
|
./pkg/obcs/obcs_calc_stevens.F |
324 |
OBCS_CALC_STEVENS |
uMerPri(K) = (OBWuStevens(J,K,bi,bj)-uMerBar) |
|
./pkg/obcs/obcs_calc_stevens.F |
326 |
OBCS_CALC_STEVENS |
& * _maskW(I+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
332 |
OBCS_CALC_STEVENS |
uMerBar = uMerBar + OBWu(J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
333 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacW(I+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
334 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacW(I+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
340 |
OBCS_CALC_STEVENS |
OBWu(J,K,bi,bj) = (uMerBar + uMerPri(K)) |
|
./pkg/obcs/obcs_calc_stevens.F |
341 |
OBCS_CALC_STEVENS |
& * _maskW(I+1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
344 |
OBCS_CALC_STEVENS |
OBWw(J,K,bi,bj)=0. |
|
./pkg/obcs/obcs_calc_stevens.F |
352 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
353 |
OBCS_CALC_STEVENS |
I = OB_Iw(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
355 |
OBCS_CALC_STEVENS |
OBWeta(J,bi,bj) = etaN(I+1,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
363 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
364 |
OBCS_CALC_STEVENS |
I = OB_Iw(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
366 |
OBCS_CALC_STEVENS |
cflMer(J,K) = 0.5 _d 0 * _dxC(I+2,J,bi,bj)/dTtracerLev(K) |
|
./pkg/obcs/obcs_calc_stevens.F |
369 |
OBCS_CALC_STEVENS |
gFacM(J,K) = ABS(MAX(SIGN(1.D0,uVel(I+1,J,K,bi,bj)),0.D0)) |
|
./pkg/obcs/obcs_calc_stevens.F |
371 |
OBCS_CALC_STEVENS |
cflMer(J,K) = 0. _d 0 |
|
./pkg/obcs/obcs_calc_stevens.F |
372 |
OBCS_CALC_STEVENS |
gFacM (J,K) = 0. _d 0 |
|
./pkg/obcs/obcs_calc_stevens.F |
413 |
OBCS_CALC_STEVENS |
J = OB_Jn(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
414 |
OBCS_CALC_STEVENS |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
423 |
OBCS_CALC_STEVENS |
vZonBar = vZonBar + vVel(I,J-1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
427 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacS(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
428 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacS(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
433 |
OBCS_CALC_STEVENS |
vZonPri(K) = (vVel(I,J-1,K,bi,bj)-vZonBar) |
|
./pkg/obcs/obcs_calc_stevens.F |
437 |
OBCS_CALC_STEVENS |
& * _maskS(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
444 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacS(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
445 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacS(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
452 |
OBCS_CALC_STEVENS |
& * _maskS(I,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
464 |
OBCS_CALC_STEVENS |
J = OB_Jn(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
465 |
OBCS_CALC_STEVENS |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
466 |
OBCS_CALC_STEVENS |
OBNeta(I,bi,bj) = etaN(I,J-1,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
475 |
OBCS_CALC_STEVENS |
J = OB_Jn(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
476 |
OBCS_CALC_STEVENS |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
477 |
OBCS_CALC_STEVENS |
cflZon(I,K) = 0.5 _d 0 * _dyC(I,J-1,bi,bj)/dTtracerLev(K) |
|
./pkg/obcs/obcs_calc_stevens.F |
480 |
OBCS_CALC_STEVENS |
gFacZ(I,K) = ABS(MIN(SIGN(1.D0,vVel(I,J,K,bi,bj)),0.D0)) |
|
./pkg/obcs/obcs_calc_stevens.F |
525 |
OBCS_CALC_STEVENS |
J = OB_Js(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
526 |
OBCS_CALC_STEVENS |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
535 |
OBCS_CALC_STEVENS |
vZonBar = vZonBar + vVel(I,J+2,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
539 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacS(I,J+1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
540 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacS(I,J+1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
545 |
OBCS_CALC_STEVENS |
vZonPri(K) = (vVel(I,J+2,K,bi,bj)-vZonBar) |
|
./pkg/obcs/obcs_calc_stevens.F |
549 |
OBCS_CALC_STEVENS |
& * _maskS(I,J+1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
556 |
OBCS_CALC_STEVENS |
& *drF(K)* _hFacS(I,J+1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
557 |
OBCS_CALC_STEVENS |
drFBar = drFBar + drF(K)* _hFacS(I,J+1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
564 |
OBCS_CALC_STEVENS |
& * _maskS(I,J+1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
576 |
OBCS_CALC_STEVENS |
J = OB_Js(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
577 |
OBCS_CALC_STEVENS |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
578 |
OBCS_CALC_STEVENS |
OBSeta(I,bi,bj) = etaN(I,J+1,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
587 |
OBCS_CALC_STEVENS |
J = OB_Js(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
588 |
OBCS_CALC_STEVENS |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
589 |
OBCS_CALC_STEVENS |
cflZon(I,K) = 0.5 _d 0 * _dyC(I,J+2,bi,bj)/dTtracerLev(K) |
|
./pkg/obcs/obcs_calc_stevens.F |
592 |
OBCS_CALC_STEVENS |
gFacZ(I,K) = ABS(MAX(SIGN(1.D0,vVel(I,J+1,K,bi,bj)),0.D0)) |
|
./pkg/obcs/obcs_calc_stevens.F |
636 |
OBCS_CALC_STEVENS |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
639 |
OBCS_CALC_STEVENS |
I = OB_Ie(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
641 |
OBCS_CALC_STEVENS |
uPhase = uPhase + OBEu(J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
642 |
OBCS_CALC_STEVENS |
& *drF(k)* _hFacW(I,J,K,bi,bj)*dyG(I,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
643 |
OBCS_CALC_STEVENS |
I = OB_Iw(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
645 |
OBCS_CALC_STEVENS |
vPhase = vPhase + OBWu(J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
646 |
OBCS_CALC_STEVENS |
& *drF(k)* _hFacW(I+1,J,K,bi,bj)*dyG(I+1,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
726 |
OBCS_STEVENS_CALC_TRACER_EAST |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
727 |
OBCS_STEVENS_CALC_TRACER_EAST |
I = OB_I(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
729 |
OBCS_STEVENS_CALC_TRACER_EAST |
dTracSpace = (tracer(I-1,J,K,bi,bj)-tracer(I-2,J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
730 |
OBCS_STEVENS_CALC_TRACER_EAST |
& * _maskW(I-1,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
731 |
OBCS_STEVENS_CALC_TRACER_EAST |
dTracTime = (tracer(I-1,J,K,bi,bj)-OBE_Stevens(J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
732 |
OBCS_STEVENS_CALC_TRACER_EAST |
uPhase = cflMer(J,K) * pFac |
|
./pkg/obcs/obcs_calc_stevens.F |
734 |
OBCS_STEVENS_CALC_TRACER_EAST |
uPhase = MIN( cflMer(J,K), |
|
./pkg/obcs/obcs_calc_stevens.F |
735 |
OBCS_STEVENS_CALC_TRACER_EAST |
& MAX( 0.D0, -cflMer(J,K)*dTracTime/dTracSpace ) |
|
./pkg/obcs/obcs_calc_stevens.F |
740 |
OBCS_STEVENS_CALC_TRACER_EAST |
OBEf(J,K,bi,bj) = _maskW(I,J,K,bi,bj) * ( |
|
./pkg/obcs/obcs_calc_stevens.F |
741 |
OBCS_STEVENS_CALC_TRACER_EAST |
& - ( aFac*MAX(0.D0,uVel(I,J,K,bi,bj)) + uPhase ) |
|
./pkg/obcs/obcs_calc_stevens.F |
742 |
OBCS_STEVENS_CALC_TRACER_EAST |
& *(tracer(I,J,K,bi,bj)-tracer(I-1,J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
743 |
OBCS_STEVENS_CALC_TRACER_EAST |
& * _recip_dxC(I,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
744 |
OBCS_STEVENS_CALC_TRACER_EAST |
& - gFacM(J,K) * gammaf |
|
./pkg/obcs/obcs_calc_stevens.F |
745 |
OBCS_STEVENS_CALC_TRACER_EAST |
& * (tracer(I,J,K,bi,bj)-OBEf(J,K,bi,bj)) ) |
|
./pkg/obcs/obcs_calc_stevens.F |
808 |
OBCS_STEVENS_CALC_TRACER_WEST |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/obcs_calc_stevens.F |
809 |
OBCS_STEVENS_CALC_TRACER_WEST |
I = OB_I(J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
811 |
OBCS_STEVENS_CALC_TRACER_WEST |
dTracSpace = (tracer(I+2,J,K,bi,bj)-tracer(I+1,J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
812 |
OBCS_STEVENS_CALC_TRACER_WEST |
& * _maskW(I+2,J,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
813 |
OBCS_STEVENS_CALC_TRACER_WEST |
dTracTime = (tracer(I+1,J,K,bi,bj)-OBW_Stevens(J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
814 |
OBCS_STEVENS_CALC_TRACER_WEST |
uPhase = -cflMer(J,K) * pFac |
|
./pkg/obcs/obcs_calc_stevens.F |
816 |
OBCS_STEVENS_CALC_TRACER_WEST |
uPhase = MAX( -cflMer(J,K), |
|
./pkg/obcs/obcs_calc_stevens.F |
817 |
OBCS_STEVENS_CALC_TRACER_WEST |
& MIN( 0.D0, -cflMer(J,K)*dTracTime/dTracSpace ) |
|
./pkg/obcs/obcs_calc_stevens.F |
822 |
OBCS_STEVENS_CALC_TRACER_WEST |
OBWf(J,K,bi,bj) = _maskW(I+1,J,K,bi,bj) * ( |
|
./pkg/obcs/obcs_calc_stevens.F |
823 |
OBCS_STEVENS_CALC_TRACER_WEST |
& - ( aFac*MIN(0.D0,uVel(I+1,J,K,bi,bj)) + uPhase ) |
|
./pkg/obcs/obcs_calc_stevens.F |
824 |
OBCS_STEVENS_CALC_TRACER_WEST |
& *(tracer(I+1,J,K,bi,bj)-tracer(I,J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
825 |
OBCS_STEVENS_CALC_TRACER_WEST |
& * _recip_dxC(I+1,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
826 |
OBCS_STEVENS_CALC_TRACER_WEST |
& - gFacM(J,K) * gammaf |
|
./pkg/obcs/obcs_calc_stevens.F |
827 |
OBCS_STEVENS_CALC_TRACER_WEST |
& * (tracer(I,J,K,bi,bj)-OBWf(J,K,bi,bj)) ) |
|
./pkg/obcs/obcs_calc_stevens.F |
890 |
OBCS_STEVENS_CALC_TRACER_NORTH |
J = OB_J(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
891 |
OBCS_STEVENS_CALC_TRACER_NORTH |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
893 |
OBCS_STEVENS_CALC_TRACER_NORTH |
dTracSpace = (tracer(I,J-1,K,bi,bj)-tracer(I,J-2,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
894 |
OBCS_STEVENS_CALC_TRACER_NORTH |
& * _maskS(I,J-1,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
895 |
OBCS_STEVENS_CALC_TRACER_NORTH |
dTracTime = (tracer(I,J-1,K,bi,bj)-OBN_Stevens(I,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
904 |
OBCS_STEVENS_CALC_TRACER_NORTH |
OBNf(I,K,bi,bj) = _maskS(I,J,K,bi,bj) * ( |
|
./pkg/obcs/obcs_calc_stevens.F |
905 |
OBCS_STEVENS_CALC_TRACER_NORTH |
& - ( aFac*MAX(0.D0,vVel(I,J,K,bi,bj)) + vPhase ) |
|
./pkg/obcs/obcs_calc_stevens.F |
906 |
OBCS_STEVENS_CALC_TRACER_NORTH |
& *(tracer(I,J,K,bi,bj)-tracer(I,J-1,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
907 |
OBCS_STEVENS_CALC_TRACER_NORTH |
& * _recip_dyC(I,J,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
909 |
OBCS_STEVENS_CALC_TRACER_NORTH |
& * (tracer(I,J,K,bi,bj)-OBNf(I,K,bi,bj)) ) |
|
./pkg/obcs/obcs_calc_stevens.F |
972 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
J = OB_J(I,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
973 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
IF ( J.NE.OB_indexNone ) THEN |
|
./pkg/obcs/obcs_calc_stevens.F |
974 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
dTracSpace = (tracer(I,J+2,K,bi,bj)-tracer(I,J+1,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
975 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
& * _maskS(I,J+2,K,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
976 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
dTracTime = (tracer(I,J+1,K,bi,bj)-OBS_Stevens(I,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
985 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
OBSf(I,K,bi,bj) = _maskS(I,J+1,K,bi,bj) * ( |
|
./pkg/obcs/obcs_calc_stevens.F |
986 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
& - ( aFac*MIN(0.D0,vVel(I,J+1,K,bi,bj)) + vPhase ) |
|
./pkg/obcs/obcs_calc_stevens.F |
987 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
& *(tracer(I,J+1,K,bi,bj)-tracer(I,J,K,bi,bj)) |
|
./pkg/obcs/obcs_calc_stevens.F |
988 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
& * _recip_dyC(I,J+1,bi,bj) |
|
./pkg/obcs/obcs_calc_stevens.F |
990 |
OBCS_STEVENS_CALC_TRACER_SOUTH |
& * (tracer(I,J,K,bi,bj)-OBSf(I,K,bi,bj)) ) |
|
./pkg/obcs/obcs_readparms.F |
586 |
OBCS_READPARMS |
OB_Iw(J,bi,bj) = OB_indexNone |
|
./pkg/obcs/orlanski_east.F |
107 |
ORLANSKI_EAST |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/orlanski_east.F |
108 |
ORLANSKI_EAST |
I_obc=OB_Ie(J,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
111 |
ORLANSKI_EAST |
IF ((UE_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_east.F |
112 |
ORLANSKI_EAST |
& (UE_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_east.F |
115 |
ORLANSKI_EAST |
CL=-(uVel(I_obc-1,J,K,bi,bj)-UE_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_east.F |
116 |
ORLANSKI_EAST |
& (ab1*UE_STORE_2(J,K,bi,bj) + ab2*UE_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_east.F |
126 |
ORLANSKI_EAST |
CVEL_UE(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_east.F |
128 |
ORLANSKI_EAST |
CVEL_UE(J,K,bi,bj) = f1*(CL*dxF(I_obc-2,J,bi,bj)/deltaT |
|
./pkg/obcs/orlanski_east.F |
129 |
ORLANSKI_EAST |
& )+f2*CVEL_UE(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
132 |
ORLANSKI_EAST |
OBEu(J,K,bi,bj)=uVel(I_obc,J,K,bi,bj)- |
|
./pkg/obcs/orlanski_east.F |
133 |
ORLANSKI_EAST |
& CVEL_UE(J,K,bi,bj)*(deltaT*recip_dxF(I_obc-1,J,bi,bj))* |
|
./pkg/obcs/orlanski_east.F |
134 |
ORLANSKI_EAST |
& (ab1*(uVel(I_obc,J,K,bi,bj)-uVel(I_obc-1,J,K,bi,bj)) + |
|
./pkg/obcs/orlanski_east.F |
135 |
ORLANSKI_EAST |
& ab2*(UE_STORE_4(J,K,bi,bj)-UE_STORE_1(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_east.F |
137 |
ORLANSKI_EAST |
IF ((VE_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_east.F |
138 |
ORLANSKI_EAST |
& (VE_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_east.F |
141 |
ORLANSKI_EAST |
CL=-(vVel(I_obc-1,J,K,bi,bj)-VE_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_east.F |
142 |
ORLANSKI_EAST |
& (ab1*VE_STORE_2(J,K,bi,bj) + ab2*VE_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_east.F |
152 |
ORLANSKI_EAST |
CVEL_VE(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_east.F |
154 |
ORLANSKI_EAST |
CVEL_VE(J,K,bi,bj) = f1*(CL*dxV(I_obc-1,J,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
155 |
ORLANSKI_EAST |
$ /deltaT)+f2*CVEL_VE(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
158 |
ORLANSKI_EAST |
OBEv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)- |
|
./pkg/obcs/orlanski_east.F |
159 |
ORLANSKI_EAST |
& CVEL_VE(J,K,bi,bj)*(deltaT*recip_dxV(I_obc,J,bi,bj))* |
|
./pkg/obcs/orlanski_east.F |
160 |
ORLANSKI_EAST |
& (ab1*(vVel(I_obc,J,K,bi,bj)-vVel(I_obc-1,J,K,bi,bj)) + |
|
./pkg/obcs/orlanski_east.F |
161 |
ORLANSKI_EAST |
& ab2*(VE_STORE_4(J,K,bi,bj)-VE_STORE_1(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_east.F |
163 |
ORLANSKI_EAST |
IF ((TE_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_east.F |
164 |
ORLANSKI_EAST |
& (TE_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_east.F |
167 |
ORLANSKI_EAST |
CL=-(theta(I_obc-1,J,K,bi,bj)-TE_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_east.F |
168 |
ORLANSKI_EAST |
& (ab1*TE_STORE_2(J,K,bi,bj) + ab2*TE_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_east.F |
178 |
ORLANSKI_EAST |
CVEL_TE(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_east.F |
180 |
ORLANSKI_EAST |
CVEL_TE(J,K,bi,bj) = f1*(CL*dxC(I_obc-1,J,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
181 |
ORLANSKI_EAST |
$ /deltaT)+f2*CVEL_TE(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
184 |
ORLANSKI_EAST |
OBEt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)- |
|
./pkg/obcs/orlanski_east.F |
185 |
ORLANSKI_EAST |
& CVEL_TE(J,K,bi,bj)*(deltaT*recip_dxC(I_obc,J,bi,bj))* |
|
./pkg/obcs/orlanski_east.F |
186 |
ORLANSKI_EAST |
& (ab1*(theta(I_obc,J,K,bi,bj)-theta(I_obc-1,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_east.F |
187 |
ORLANSKI_EAST |
& ab2*(TE_STORE_4(J,K,bi,bj)-TE_STORE_1(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_east.F |
189 |
ORLANSKI_EAST |
IF ((SE_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_east.F |
190 |
ORLANSKI_EAST |
& (SE_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_east.F |
193 |
ORLANSKI_EAST |
CL=-(salt(I_obc-1,J,K,bi,bj)-SE_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_east.F |
194 |
ORLANSKI_EAST |
& (ab1*SE_STORE_2(J,K,bi,bj) + ab2*SE_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_east.F |
204 |
ORLANSKI_EAST |
CVEL_SE(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_east.F |
206 |
ORLANSKI_EAST |
CVEL_SE(J,K,bi,bj) = f1*(CL*dxC(I_obc-1,J,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
207 |
ORLANSKI_EAST |
$ /deltaT)+f2*CVEL_SE(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
210 |
ORLANSKI_EAST |
OBEs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)- |
|
./pkg/obcs/orlanski_east.F |
211 |
ORLANSKI_EAST |
& CVEL_SE(J,K,bi,bj)*(deltaT*recip_dxC(I_obc,J,bi,bj))* |
|
./pkg/obcs/orlanski_east.F |
212 |
ORLANSKI_EAST |
& (ab1*(salt(I_obc,J,K,bi,bj)-salt(I_obc-1,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_east.F |
213 |
ORLANSKI_EAST |
& ab2*(SE_STORE_4(J,K,bi,bj)-SE_STORE_1(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_east.F |
217 |
ORLANSKI_EAST |
IF ((WE_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_east.F |
218 |
ORLANSKI_EAST |
& (WE_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_east.F |
221 |
ORLANSKI_EAST |
CL=-(wVel(I_obc-1,J,K,bi,bj)-WE_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_east.F |
222 |
ORLANSKI_EAST |
& (ab1*WE_STORE_2(J,K,bi,bj)+ab2*WE_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_east.F |
232 |
ORLANSKI_EAST |
CVEL_WE(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_east.F |
234 |
ORLANSKI_EAST |
CVEL_WE(J,K,bi,bj)=f1*(CL*dxC(I_obc-1,J,bi,bj)/deltaT) |
|
./pkg/obcs/orlanski_east.F |
235 |
ORLANSKI_EAST |
& + f2*CVEL_WE(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
238 |
ORLANSKI_EAST |
OBEw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)- |
|
./pkg/obcs/orlanski_east.F |
239 |
ORLANSKI_EAST |
& CVEL_WE(J,K,bi,bj)*(deltaT*recip_dxC(I_obc,J,bi,bj))* |
|
./pkg/obcs/orlanski_east.F |
240 |
ORLANSKI_EAST |
& (ab1*(wVel(I_obc,J,K,bi,bj)-wVel(I_obc-1,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_east.F |
241 |
ORLANSKI_EAST |
& ab2*(WE_STORE_4(J,K,bi,bj)-WE_STORE_1(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_east.F |
247 |
ORLANSKI_EAST |
UE_STORE_3(J,K,bi,bj)=UE_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
249 |
ORLANSKI_EAST |
UE_STORE_2(J,K,bi,bj)=uVel(I_obc-1,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_east.F |
250 |
ORLANSKI_EAST |
& uVel(I_obc-2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
251 |
ORLANSKI_EAST |
UE_STORE_1(J,K,bi,bj)=uVel(I_obc-1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
252 |
ORLANSKI_EAST |
UE_STORE_4(J,K,bi,bj)=uVel(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
255 |
ORLANSKI_EAST |
VE_STORE_3(J,K,bi,bj)=VE_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
257 |
ORLANSKI_EAST |
VE_STORE_2(J,K,bi,bj)=vVel(I_obc-1,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_east.F |
258 |
ORLANSKI_EAST |
& vVel(I_obc-2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
259 |
ORLANSKI_EAST |
VE_STORE_1(J,K,bi,bj)=vVel(I_obc-1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
260 |
ORLANSKI_EAST |
VE_STORE_4(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
263 |
ORLANSKI_EAST |
TE_STORE_3(J,K,bi,bj)=TE_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
265 |
ORLANSKI_EAST |
TE_STORE_2(J,K,bi,bj)=theta(I_obc-1,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_east.F |
266 |
ORLANSKI_EAST |
& theta(I_obc-2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
267 |
ORLANSKI_EAST |
TE_STORE_1(J,K,bi,bj)=theta(I_obc-1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
268 |
ORLANSKI_EAST |
TE_STORE_4(J,K,bi,bj)=theta(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
271 |
ORLANSKI_EAST |
SE_STORE_3(J,K,bi,bj)=SE_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
273 |
ORLANSKI_EAST |
SE_STORE_2(J,K,bi,bj)=salt(I_obc-1,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_east.F |
274 |
ORLANSKI_EAST |
& salt(I_obc-2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
275 |
ORLANSKI_EAST |
SE_STORE_1(J,K,bi,bj)=salt(I_obc-1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
276 |
ORLANSKI_EAST |
SE_STORE_4(J,K,bi,bj)=salt(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
281 |
ORLANSKI_EAST |
WE_STORE_3(J,K,bi,bj)=WE_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
283 |
ORLANSKI_EAST |
WE_STORE_2(J,K,bi,bj)=wVel(I_obc-1,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_east.F |
284 |
ORLANSKI_EAST |
& wVel(I_obc-2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
285 |
ORLANSKI_EAST |
WE_STORE_1(J,K,bi,bj)=wVel(I_obc-1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
286 |
ORLANSKI_EAST |
WE_STORE_4(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_east.F |
93 |
ORLANSKI_EAST |
INTEGER J, K, I_obc |
|
./pkg/obcs/orlanski_init.F |
100 |
ORLANSKI_INIT |
UE_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
101 |
ORLANSKI_INIT |
VE_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
102 |
ORLANSKI_INIT |
TE_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
103 |
ORLANSKI_INIT |
SE_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
104 |
ORLANSKI_INIT |
UW_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
105 |
ORLANSKI_INIT |
VW_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
106 |
ORLANSKI_INIT |
TW_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
107 |
ORLANSKI_INIT |
SW_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
108 |
ORLANSKI_INIT |
UW_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
109 |
ORLANSKI_INIT |
VW_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
110 |
ORLANSKI_INIT |
TW_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
111 |
ORLANSKI_INIT |
SW_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
112 |
ORLANSKI_INIT |
UW_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
113 |
ORLANSKI_INIT |
VW_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
114 |
ORLANSKI_INIT |
TW_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
115 |
ORLANSKI_INIT |
SW_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
116 |
ORLANSKI_INIT |
UW_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
117 |
ORLANSKI_INIT |
VW_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
118 |
ORLANSKI_INIT |
TW_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
119 |
ORLANSKI_INIT |
SW_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
120 |
ORLANSKI_INIT |
CVEL_UE(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
121 |
ORLANSKI_INIT |
CVEL_VE(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
122 |
ORLANSKI_INIT |
CVEL_TE(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
123 |
ORLANSKI_INIT |
CVEL_SE(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
124 |
ORLANSKI_INIT |
CVEL_UW(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
125 |
ORLANSKI_INIT |
CVEL_VW(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
126 |
ORLANSKI_INIT |
CVEL_TW(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
127 |
ORLANSKI_INIT |
CVEL_SW(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
129 |
ORLANSKI_INIT |
WE_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
130 |
ORLANSKI_INIT |
WE_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
131 |
ORLANSKI_INIT |
WE_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
132 |
ORLANSKI_INIT |
WE_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
133 |
ORLANSKI_INIT |
WW_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
134 |
ORLANSKI_INIT |
WW_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
135 |
ORLANSKI_INIT |
WW_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
136 |
ORLANSKI_INIT |
WW_STORE_4(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
137 |
ORLANSKI_INIT |
CVEL_WE(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
138 |
ORLANSKI_INIT |
CVEL_WW(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
30 |
ORLANSKI_INIT |
INTEGER I, J, K |
|
./pkg/obcs/orlanski_init.F |
87 |
ORLANSKI_INIT |
DO J=1-Oly,sNy+Oly |
|
./pkg/obcs/orlanski_init.F |
88 |
ORLANSKI_INIT |
UE_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
89 |
ORLANSKI_INIT |
VE_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
90 |
ORLANSKI_INIT |
TE_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
91 |
ORLANSKI_INIT |
SE_STORE_1(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
92 |
ORLANSKI_INIT |
UE_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
93 |
ORLANSKI_INIT |
VE_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
94 |
ORLANSKI_INIT |
TE_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
95 |
ORLANSKI_INIT |
SE_STORE_2(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
96 |
ORLANSKI_INIT |
UE_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
97 |
ORLANSKI_INIT |
VE_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
98 |
ORLANSKI_INIT |
TE_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_init.F |
99 |
ORLANSKI_INIT |
SE_STORE_3(J,K,bi,bj)=0. |
|
./pkg/obcs/orlanski_west.F |
109 |
ORLANSKI_WEST |
DO J=1-OLy,sNy+OLy |
|
./pkg/obcs/orlanski_west.F |
110 |
ORLANSKI_WEST |
I_obc=OB_Iw(J,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
113 |
ORLANSKI_WEST |
IF ((UW_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_west.F |
114 |
ORLANSKI_WEST |
& (UW_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_west.F |
117 |
ORLANSKI_WEST |
CL=(uVel(I_obc+2,J,K,bi,bj)-UW_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_west.F |
118 |
ORLANSKI_WEST |
& (ab1*UW_STORE_2(J,K,bi,bj) + ab2*UW_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_west.F |
128 |
ORLANSKI_WEST |
CVEL_UW(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_west.F |
130 |
ORLANSKI_WEST |
CVEL_UW(J,K,bi,bj) = f1*(CL*dxF(I_obc+2,J,bi,bj)/deltaT |
|
./pkg/obcs/orlanski_west.F |
131 |
ORLANSKI_WEST |
& )+f2*CVEL_UW(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
134 |
ORLANSKI_WEST |
OBWu(J,K,bi,bj)=uVel(I_obc+1,J,K,bi,bj)+ |
|
./pkg/obcs/orlanski_west.F |
135 |
ORLANSKI_WEST |
& CVEL_UW(J,K,bi,bj)*(deltaT*recip_dxF(I_obc+1,J,bi,bj))* |
|
./pkg/obcs/orlanski_west.F |
136 |
ORLANSKI_WEST |
& (ab1*(uVel(I_obc+2,J,K,bi,bj)-uVel(I_obc+1,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_west.F |
137 |
ORLANSKI_WEST |
& ab2*(UW_STORE_1(J,K,bi,bj)-UW_STORE_4(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_west.F |
139 |
ORLANSKI_WEST |
IF ((VW_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_west.F |
140 |
ORLANSKI_WEST |
& (VW_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_west.F |
143 |
ORLANSKI_WEST |
CL=(vVel(I_obc+1,J,K,bi,bj)-VW_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_west.F |
144 |
ORLANSKI_WEST |
& (ab1*VW_STORE_2(J,K,bi,bj) + ab2*VW_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_west.F |
154 |
ORLANSKI_WEST |
CVEL_VW(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_west.F |
156 |
ORLANSKI_WEST |
CVEL_VW(J,K,bi,bj) = f1*(CL*dxV(I_obc+2,J,bi,bj)/deltaT |
|
./pkg/obcs/orlanski_west.F |
157 |
ORLANSKI_WEST |
& )+f2*CVEL_VW(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
160 |
ORLANSKI_WEST |
OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+ |
|
./pkg/obcs/orlanski_west.F |
161 |
ORLANSKI_WEST |
& CVEL_VW(J,K,bi,bj)*(deltaT*recip_dxV(I_obc+1,J,bi,bj))* |
|
./pkg/obcs/orlanski_west.F |
162 |
ORLANSKI_WEST |
& (ab1*(vVel(I_obc+1,J,K,bi,bj)-vVel(I_obc,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_west.F |
163 |
ORLANSKI_WEST |
& ab2*(VW_STORE_1(J,K,bi,bj)-VW_STORE_4(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_west.F |
165 |
ORLANSKI_WEST |
IF ((TW_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_west.F |
166 |
ORLANSKI_WEST |
& (TW_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_west.F |
169 |
ORLANSKI_WEST |
CL=(theta(I_obc+1,J,K,bi,bj)-TW_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_west.F |
170 |
ORLANSKI_WEST |
& (ab1*TW_STORE_2(J,K,bi,bj) + ab2*TW_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_west.F |
180 |
ORLANSKI_WEST |
CVEL_TW(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_west.F |
182 |
ORLANSKI_WEST |
CVEL_TW(J,K,bi,bj) = f1*(CL*dxC(I_obc+2,J,bi,bj)/deltaT |
|
./pkg/obcs/orlanski_west.F |
183 |
ORLANSKI_WEST |
& )+f2*CVEL_TW(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
186 |
ORLANSKI_WEST |
OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+ |
|
./pkg/obcs/orlanski_west.F |
187 |
ORLANSKI_WEST |
& CVEL_TW(J,K,bi,bj)*(deltaT*recip_dxC(I_obc+1,J,bi,bj))* |
|
./pkg/obcs/orlanski_west.F |
188 |
ORLANSKI_WEST |
& (ab1*(theta(I_obc+1,J,K,bi,bj)-theta(I_obc,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_west.F |
189 |
ORLANSKI_WEST |
& ab2*(TW_STORE_1(J,K,bi,bj)-TW_STORE_4(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_west.F |
191 |
ORLANSKI_WEST |
IF ((SW_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_west.F |
192 |
ORLANSKI_WEST |
& (SW_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_west.F |
195 |
ORLANSKI_WEST |
CL=(salt(I_obc+1,J,K,bi,bj)-SW_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_west.F |
196 |
ORLANSKI_WEST |
& (ab1*SW_STORE_2(J,K,bi,bj) + ab2*SW_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_west.F |
206 |
ORLANSKI_WEST |
CVEL_SW(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_west.F |
208 |
ORLANSKI_WEST |
CVEL_SW(J,K,bi,bj) = f1*(CL*dxC(I_obc+2,J,bi,bj)/deltaT |
|
./pkg/obcs/orlanski_west.F |
209 |
ORLANSKI_WEST |
& )+f2*CVEL_SW(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
212 |
ORLANSKI_WEST |
OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+ |
|
./pkg/obcs/orlanski_west.F |
213 |
ORLANSKI_WEST |
& CVEL_SW(J,K,bi,bj)*(deltaT*recip_dxC(I_obc+1,J,bi,bj))* |
|
./pkg/obcs/orlanski_west.F |
214 |
ORLANSKI_WEST |
& (ab1*(salt(I_obc+1,J,K,bi,bj)-salt(I_obc,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_west.F |
215 |
ORLANSKI_WEST |
& ab2*(SW_STORE_1(J,K,bi,bj)-SW_STORE_4(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_west.F |
219 |
ORLANSKI_WEST |
IF ((WW_STORE_2(J,K,bi,bj).eq.0.).and. |
|
./pkg/obcs/orlanski_west.F |
220 |
ORLANSKI_WEST |
& (WW_STORE_3(J,K,bi,bj).eq.0.)) THEN |
|
./pkg/obcs/orlanski_west.F |
223 |
ORLANSKI_WEST |
CL=(wVel(I_obc+1,J,K,bi,bj)-WW_STORE_1(J,K,bi,bj))/ |
|
./pkg/obcs/orlanski_west.F |
224 |
ORLANSKI_WEST |
& (ab1*WW_STORE_2(J,K,bi,bj)+ab2*WW_STORE_3(J,K,bi,bj)) |
|
./pkg/obcs/orlanski_west.F |
234 |
ORLANSKI_WEST |
CVEL_WW(J,K,bi,bj) = CFIX |
|
./pkg/obcs/orlanski_west.F |
236 |
ORLANSKI_WEST |
CVEL_WW(J,K,bi,bj)=f1*(CL*dxC(I_obc+2,J,bi,bj)/deltaT) |
|
./pkg/obcs/orlanski_west.F |
237 |
ORLANSKI_WEST |
& + f2*CVEL_WW(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
240 |
ORLANSKI_WEST |
OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+ |
|
./pkg/obcs/orlanski_west.F |
241 |
ORLANSKI_WEST |
& CVEL_WW(J,K,bi,bj)*(deltaT*recip_dxC(I_obc+1,J,bi,bj))* |
|
./pkg/obcs/orlanski_west.F |
242 |
ORLANSKI_WEST |
& (ab1*(wVel(I_obc+1,J,K,bi,bj)-wVel(I_obc,J,K,bi,bj))+ |
|
./pkg/obcs/orlanski_west.F |
243 |
ORLANSKI_WEST |
& ab2*(WW_STORE_1(J,K,bi,bj)-WW_STORE_4(J,K,bi,bj))) |
|
./pkg/obcs/orlanski_west.F |
249 |
ORLANSKI_WEST |
UW_STORE_3(J,K,bi,bj)=UW_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
251 |
ORLANSKI_WEST |
UW_STORE_2(J,K,bi,bj)=uVel(I_obc+3,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_west.F |
252 |
ORLANSKI_WEST |
& uVel(I_obc+2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
253 |
ORLANSKI_WEST |
UW_STORE_1(J,K,bi,bj)=uVel(I_obc+2,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
254 |
ORLANSKI_WEST |
UW_STORE_4(J,K,bi,bj)=uVel(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
257 |
ORLANSKI_WEST |
VW_STORE_3(J,K,bi,bj)=VW_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
259 |
ORLANSKI_WEST |
VW_STORE_2(J,K,bi,bj)=vVel(I_obc+2,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_west.F |
260 |
ORLANSKI_WEST |
& vVel(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
261 |
ORLANSKI_WEST |
VW_STORE_1(J,K,bi,bj)=vVel(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
262 |
ORLANSKI_WEST |
VW_STORE_4(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
265 |
ORLANSKI_WEST |
TW_STORE_3(J,K,bi,bj)=TW_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
267 |
ORLANSKI_WEST |
TW_STORE_2(J,K,bi,bj)=theta(I_obc+2,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_west.F |
268 |
ORLANSKI_WEST |
& theta(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
269 |
ORLANSKI_WEST |
TW_STORE_1(J,K,bi,bj)=theta(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
270 |
ORLANSKI_WEST |
TW_STORE_4(J,K,bi,bj)=theta(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
273 |
ORLANSKI_WEST |
SW_STORE_3(J,K,bi,bj)=SW_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
275 |
ORLANSKI_WEST |
SW_STORE_2(J,K,bi,bj)=salt(I_obc+2,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_west.F |
276 |
ORLANSKI_WEST |
& salt(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
277 |
ORLANSKI_WEST |
SW_STORE_1(J,K,bi,bj)=salt(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
278 |
ORLANSKI_WEST |
SW_STORE_4(J,K,bi,bj)=salt(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
283 |
ORLANSKI_WEST |
WW_STORE_3(J,K,bi,bj)=WW_STORE_2(J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
285 |
ORLANSKI_WEST |
WW_STORE_2(J,K,bi,bj)=wVel(I_obc+2,J,K,bi,bj) - |
|
./pkg/obcs/orlanski_west.F |
286 |
ORLANSKI_WEST |
& wVel(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
287 |
ORLANSKI_WEST |
WW_STORE_1(J,K,bi,bj)=wVel(I_obc+1,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
288 |
ORLANSKI_WEST |
WW_STORE_4(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj) |
|
./pkg/obcs/orlanski_west.F |
95 |
ORLANSKI_WEST |
INTEGER J, K, I_obc |
|
./pkg/ocn_compon_interf/cpl_ini_vars.F |
49 |
CPL_INI_VARS |
DO J=1-OLy,sNy+OLy |
|
./pkg/opps/opps_calc.F |
22 |
OPPS_CALC |
I I, J, bi, bj, myTime, myIter, myThid ) |
|
./pkg/opps/opps_calc.F |
420 |
OPPS_CALC |
OPPSconvectCount(I,J,K,bi,bj) = |
|
./pkg/opps/opps_calc.F |
421 |
OPPS_CALC |
& OPPSconvectCount(I,J,K,bi,bj) + 1. _d 0 |
|
./pkg/opps/opps_calc.F |
81 |
OPPS_CALC |
INTEGER I, J, bi, bj, KMax, nTracer, nTracerInUse |
|
./pkg/opps/opps_init.F |
42 |
OPPS_INIT |
OPPSconvectCount(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/opps/opps_interface.F |
103 |
OPPS_INTERFACE |
DO J=jMin,jMax |
|
./pkg/opps/opps_interface.F |
105 |
OPPS_INTERFACE |
IF ( kSurfC(I,J,bi,bj) .LE. Nr ) THEN |
|
./pkg/opps/opps_interface.F |
108 |
OPPS_INTERFACE |
tracerLoc(K,1) = theta(I,J,K,bi,bj) |
|
./pkg/opps/opps_interface.F |
109 |
OPPS_INTERFACE |
tracerLoc(K,2) = salt(I,J,K,bi,bj) |
|
./pkg/opps/opps_interface.F |
110 |
OPPS_INTERFACE |
wVelLoc(K) = wVel(I,J,K,bi,bj) |
|
./pkg/opps/opps_interface.F |
114 |
OPPS_INTERFACE |
tracerLoc(K,1) = theta(I,J,K,bi,bj) |
|
./pkg/opps/opps_interface.F |
115 |
OPPS_INTERFACE |
tracerLoc(K,2) = salt(I,J,K,bi,bj) |
|
./pkg/opps/opps_interface.F |
122 |
OPPS_INTERFACE |
tracerLoc(K,itr) = ptracer(I,J,K,bi,bj,itr-2) |
|
./pkg/opps/opps_interface.F |
141 |
OPPS_INTERFACE |
kMax = kLowC(I,J,bi,bj) |
|
./pkg/opps/opps_interface.F |
145 |
OPPS_INTERFACE |
I I,J,bi,bj,myTime,myIter,myThid) |
|
./pkg/opps/opps_interface.F |
165 |
OPPS_INTERFACE |
WRITE(msgBuf,'(A,2I5)') ' for (i,j) = ', I,J |
|
./pkg/opps/opps_interface.F |
172 |
OPPS_INTERFACE |
theta(I,J,K,bi,bj) = tracerLoc(K,1) |
|
./pkg/opps/opps_interface.F |
173 |
OPPS_INTERFACE |
salt(I,J,K,bi,bj) = tracerLoc(K,2) |
|
./pkg/opps/opps_interface.F |
178 |
OPPS_INTERFACE |
ptracer(I,J,K,bi,bj,itr-2) = tracerLoc(K,itr) |
|
./pkg/opps/opps_interface.F |
69 |
OPPS_INTERFACE |
INTEGER I, J, K, kMax |
|
./pkg/opps/opps_interface.F |
95 |
OPPS_INTERFACE |
DO J=1-Oly,sNy+Oly |
|
./pkg/opps/opps_interface.F |
97 |
OPPS_INTERFACE |
OPPSconvectCount(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/pp81/pp81_calc.F |
103 |
PP81_CALC |
DO J=jMin,jMax |
|
./pkg/pp81/pp81_calc.F |
105 |
PP81_CALC |
PPviscAr(I,J,K,bi,bj) = MAX(PPviscMin,PPviscAr(I,J,K,bi,bj)) |
|
./pkg/pp81/pp81_calc.F |
106 |
PP81_CALC |
PPdiffKr(I,J,K,bi,bj) = MAX(PPdiffMin,PPdiffKr(I,J,K,bi,bj)) |
|
./pkg/pp81/pp81_calc.F |
112 |
PP81_CALC |
DO J=jMin,jMax |
|
./pkg/pp81/pp81_calc.F |
114 |
PP81_CALC |
PPviscAr(I,J,K,bi,bj) = PPviscAr(I,J,K,bi,bj) |
|
./pkg/pp81/pp81_calc.F |
115 |
PP81_CALC |
& * maskC(I,J,K,bi,bj) |
|
./pkg/pp81/pp81_calc.F |
116 |
PP81_CALC |
PPdiffKr(I,J,K,bi,bj) = PPdiffKr(I,J,K,bi,bj) |
|
./pkg/pp81/pp81_calc.F |
117 |
PP81_CALC |
& * maskC(I,J,K,bi,bj) |
|
./pkg/pp81/pp81_calc.F |
60 |
PP81_CALC |
INTEGER I, J, K |
|
./pkg/pp81/pp81_calc.F |
75 |
PP81_CALC |
DO J=jMin,jMax |
|
./pkg/pp81/pp81_calc.F |
77 |
PP81_CALC |
IF ( RiNumber(I,J) .LT. RiLimit ) THEN |
|
./pkg/pp81/pp81_calc.F |
81 |
PP81_CALC |
denom = 1.0 + PPalpha*RiNumber(I,J) |
|
./pkg/pp81/pp81_calc.F |
85 |
PP81_CALC |
PPviscAr(I,J,K,bi,bj) = MAX(PPviscTmp,viscArNr(k)) |
|
./pkg/pp81/pp81_calc.F |
86 |
PP81_CALC |
PPdiffKr(I,J,K,bi,bj) = MAX(PPviscAr(I,J,K,bi,bj)/denom, |
|
./pkg/pp81/pp81_ri_number.F |
106 |
PP81_RI_NUMBER |
DO J= jMin, jMax |
|
./pkg/pp81/pp81_ri_number.F |
108 |
PP81_RI_NUMBER |
tempu= .5*( uVel(I,J,Km1,bi,bj)+uVel(I+1,J,Km1,bi,bj) |
|
./pkg/pp81/pp81_ri_number.F |
109 |
PP81_RI_NUMBER |
& - (uVel(I,J,K ,bi,bj)+uVel(I+1,J,K ,bi,bj)) ) |
|
./pkg/pp81/pp81_ri_number.F |
111 |
PP81_RI_NUMBER |
tempv= .5*( vVel(I,J,Km1,bi,bj)+vVel(I,J+1,Km1,bi,bj) |
|
./pkg/pp81/pp81_ri_number.F |
112 |
PP81_RI_NUMBER |
& - (vVel(I,J,K ,bi,bj)+vVel(I,J+1,K ,bi,bj)) ) |
|
./pkg/pp81/pp81_ri_number.F |
120 |
PP81_RI_NUMBER |
& (rhoKm1(I,J) - rhoK(I,J))*recip_drC(K) |
|
./pkg/pp81/pp81_ri_number.F |
125 |
PP81_RI_NUMBER |
RiNumber(I,J) = buoyFreq/max(RiFlux,epsilon) |
|
./pkg/pp81/pp81_ri_number.F |
131 |
PP81_RI_NUMBER |
DO J=jMin,jMax |
|
./pkg/pp81/pp81_ri_number.F |
133 |
PP81_RI_NUMBER |
RiTmp(I,J) = RiNumber(I,J) |
|
./pkg/pp81/pp81_ri_number.F |
136 |
PP81_RI_NUMBER |
DO J=1-Oly+1,sNy+Oly-1 |
|
./pkg/pp81/pp81_ri_number.F |
138 |
PP81_RI_NUMBER |
RiNumber(I,J) = p5*RiTmp(I,J) |
|
./pkg/pp81/pp81_ri_number.F |
139 |
PP81_RI_NUMBER |
& + p125*RiTmp(I-1,J) + p125*RiTmp(I+1,J) |
|
./pkg/pp81/pp81_ri_number.F |
140 |
PP81_RI_NUMBER |
& + p125*RiTmp(I,J-1) + p125*RiTmp(I,J+1) |
|
./pkg/pp81/pp81_ri_number.F |
62 |
PP81_RI_NUMBER |
INTEGER I,J,Km1 |
|
./pkg/salt_plume/salt_plume_calc_depth.F |
143 |
SALT_PLUME_CALC_DEPTH |
SaltPlumeDepth(i,j,bi,bj) = rF(1)-R_low(I,J,bi,bj) |
|
./pkg/salt_plume/salt_plume_calc_depth.F |
82 |
SALT_PLUME_CALC_DEPTH |
SaltPlumeDepth(i,j,bi,bj) = rF(1)-R_low(I,J,bi,bj) |
|
./pkg/salt_plume/salt_plume_calc_depth.F |
95 |
SALT_PLUME_CALC_DEPTH |
SaltPlumeDepth(i,j,bi,bj) = rF(1)-R_low(I,J,bi,bj) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
107 |
SALT_PLUME_TENDENCY_APPLY_S |
plumetend(I,J) = saltPlumeFlux(i,j,bi,bj)*plumefrac(I,J) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
109 |
SALT_PLUME_TENDENCY_APPLY_S |
plumetend(I,J) = SPforcingS(i,j,k,bi,bj) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
111 |
SALT_PLUME_TENDENCY_APPLY_S |
gS_arr(i,j) = gS_arr(i,j) + plumetend(I,J) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
115 |
SALT_PLUME_TENDENCY_APPLY_S |
plumetend(I,J) = 0. _d 0 |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
125 |
SALT_PLUME_TENDENCY_APPLY_S |
plumefrac(I,J) = 0. _d 0 |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
144 |
SALT_PLUME_TENDENCY_APPLY_S |
plumefrac(I,J) = (plumekb(2)-plumekb(1))*maskC(i,j,k,bi,bj) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
145 |
SALT_PLUME_TENDENCY_APPLY_S |
plumetend(I,J) = saltPlumeFlux(i,j,bi,bj)*plumefrac(I,J) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
149 |
SALT_PLUME_TENDENCY_APPLY_S |
gS_arr(i,j) = gS_arr(i,j) + plumetend(I,J) |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
152 |
SALT_PLUME_TENDENCY_APPLY_S |
plumefrac(I,J) = 0. _d 0 |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
153 |
SALT_PLUME_TENDENCY_APPLY_S |
plumetend(I,J) = 0. _d 0 |
|
./pkg/salt_plume/salt_plume_tendency_apply_s.F |
88 |
SALT_PLUME_TENDENCY_APPLY_S |
plumefrac(I,J) = plumekb2D(i,j) |
|
./pkg/salt_plume/salt_plume_tendency_apply_t.F |
60 |
SALT_PLUME_TENDENCY_APPLY_T |
gT_arr(i,j) = gT_arr(i,j) + plumeEtend(I,J) |
|
./pkg/seaice/advect.F |
109 |
ADVECT |
DO J=1-Oly,sNy+Oly |
|
./pkg/seaice/dynsolver.F |
107 |
DYNSOLVER |
AMASS(I,J,bi,bj)=RHOICE*QUART*( |
|
./pkg/seaice/dynsolver.F |
110 |
DYNSOLVER |
COR_ICE(I,J,bi,bj)=AMASS(I,J,bi,bj) * _fCoriG(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
124 |
DYNSOLVER |
U1=QUART*(UWIND(I-1,J-1,bi,bj)+UWIND(I-1,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
125 |
DYNSOLVER |
& +UWIND(I ,J-1,bi,bj)+UWIND(I ,J,bi,bj)) |
|
./pkg/seaice/dynsolver.F |
126 |
DYNSOLVER |
V1=QUART*(VWIND(I-1,J-1,bi,bj)+VWIND(I-1,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
127 |
DYNSOLVER |
& +VWIND(I ,J-1,bi,bj)+VWIND(I ,J,bi,bj)) |
|
./pkg/seaice/dynsolver.F |
135 |
DYNSOLVER |
DAIRN(I,J,bi,bj)=RHOAIR*OCEAN_drag |
|
./pkg/seaice/dynsolver.F |
137 |
DYNSOLVER |
WINDX(I,J,bi,bj)=DAIRN(I,J,bi,bj)* |
|
./pkg/seaice/dynsolver.F |
138 |
DYNSOLVER |
& (COSWIN*U1-SIGN(SINWIN, _fCori(I,J,bi,bj))*V1) |
|
./pkg/seaice/dynsolver.F |
139 |
DYNSOLVER |
WINDY(I,J,bi,bj)=DAIRN(I,J,bi,bj)* |
|
./pkg/seaice/dynsolver.F |
140 |
DYNSOLVER |
& (SIGN(SINWIN, _fCori(I,J,bi,bj))*U1+COSWIN*V1) |
|
./pkg/seaice/dynsolver.F |
143 |
DYNSOLVER |
IF ( YC(I,J,bi,bj) .LT. ZERO ) THEN |
|
./pkg/seaice/dynsolver.F |
144 |
DYNSOLVER |
DAIRN(I,J,bi,bj) = |
|
./pkg/seaice/dynsolver.F |
145 |
DYNSOLVER |
& RHOAIR*(SEAICE_drag_south*AAA*AREA(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
147 |
DYNSOLVER |
& +0.0764 _d 0*AAA*AAA)*(ONE-AREA(I,J,bi,bj))) |
|
./pkg/seaice/dynsolver.F |
149 |
DYNSOLVER |
DAIRN(I,J,bi,bj) = |
|
./pkg/seaice/dynsolver.F |
150 |
DYNSOLVER |
& RHOAIR*(SEAICE_drag*AAA*AREA(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
152 |
DYNSOLVER |
& +0.0764 _d 0*AAA*AAA)*(ONE-AREA(I,J,bi,bj))) |
|
./pkg/seaice/dynsolver.F |
154 |
DYNSOLVER |
FORCEX(I,J,bi,bj)=DAIRN(I,J,bi,bj)* |
|
./pkg/seaice/dynsolver.F |
155 |
DYNSOLVER |
& (COSWIN*U1-SIGN(SINWIN, _fCori(I,J,bi,bj))*V1) |
|
./pkg/seaice/dynsolver.F |
156 |
DYNSOLVER |
FORCEY(I,J,bi,bj)=DAIRN(I,J,bi,bj)* |
|
./pkg/seaice/dynsolver.F |
157 |
DYNSOLVER |
& (SIGN(SINWIN, _fCori(I,J,bi,bj))*U1+COSWIN*V1) |
|
./pkg/seaice/dynsolver.F |
196 |
DYNSOLVER |
FORCEX(I,J,bi,bj)=FORCEX(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
197 |
DYNSOLVER |
& -AMASS(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
200 |
DYNSOLVER |
& )*HALF*_recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
201 |
DYNSOLVER |
FORCEY(I,J,bi,bj)=FORCEY(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
202 |
DYNSOLVER |
& -AMASS(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
205 |
DYNSOLVER |
& )*HALF*_recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
207 |
DYNSOLVER |
FORCEX0(I,J,bi,bj)=FORCEX(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
208 |
DYNSOLVER |
FORCEY0(I,J,bi,bj)=FORCEY(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
216 |
DYNSOLVER |
FORCEX(I,J,bi,bj)=FORCEX(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
217 |
DYNSOLVER |
& -COR_ICE(I,J,bi,bj)*GWATY(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
218 |
DYNSOLVER |
FORCEY(I,J,bi,bj)=FORCEY(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
219 |
DYNSOLVER |
& +COR_ICE(I,J,bi,bj)*GWATX(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
221 |
DYNSOLVER |
FORCEX0(I,J,bi,bj)=FORCEX(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
222 |
DYNSOLVER |
FORCEY0(I,J,bi,bj)=FORCEY(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
225 |
DYNSOLVER |
PRESS0(I,J,bi,bj)=PSTAR*HEFF(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
228 |
DYNSOLVER |
ZMAX(I,J,bi,bj)=SEAICE_zetaMaxFac*PRESS0(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
230 |
DYNSOLVER |
ZMIN(I,J,bi,bj)=SEAICE_zetaMin |
|
./pkg/seaice/dynsolver.F |
231 |
DYNSOLVER |
PRESS0(I,J,bi,bj)=PRESS0(I,J,bi,bj)*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
257 |
DYNSOLVER |
UICENM1(I,J,bi,bj)=UICE(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
258 |
DYNSOLVER |
VICENM1(I,J,bi,bj)=VICE(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
259 |
DYNSOLVER |
UICEC(I,J,bi,bj)=UICE(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
260 |
DYNSOLVER |
VICEC(I,J,bi,bj)=VICE(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
278 |
DYNSOLVER |
UICE(I,J,bi,bj)=HALF*(UICE(I,J,bi,bj)+UICENM1(I,J,bi,bj)) |
|
./pkg/seaice/dynsolver.F |
279 |
DYNSOLVER |
VICE(I,J,bi,bj)=HALF*(VICE(I,J,bi,bj)+VICENM1(I,J,bi,bj)) |
|
./pkg/seaice/dynsolver.F |
280 |
DYNSOLVER |
UICEC(I,J,bi,bj)=UICE(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
281 |
DYNSOLVER |
VICEC(I,J,bi,bj)=VICE(I,J,bi,bj) |
|
./pkg/seaice/dynsolver.F |
81 |
DYNSOLVER |
GWATX(I,J,bi,bj)=HALF*(uVel(i,j,KGEO(I,J,bi,bj),bi,bj) |
|
./pkg/seaice/dynsolver.F |
82 |
DYNSOLVER |
& +uVel(i,j-1,KGEO(I,J,bi,bj),bi,bj)) |
|
./pkg/seaice/dynsolver.F |
83 |
DYNSOLVER |
GWATY(I,J,bi,bj)=HALF*(vVel(i,j,KGEO(I,J,bi,bj),bi,bj) |
|
./pkg/seaice/dynsolver.F |
84 |
DYNSOLVER |
& +vVel(i-1,j,KGEO(I,J,bi,bj),bi,bj)) |
|
./pkg/seaice/lsr.F |
1000 |
LSR |
DO J=1-OLy,sNy+OLy-1 |
|
./pkg/seaice/lsr.F |
1002 |
LSR |
dUdx(I,J) = 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
1003 |
LSR |
& ( UICEC(I+1,J+1,bi,bj) - UICEC(I ,J+1,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1004 |
LSR |
& * _recip_dxG(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
1005 |
LSR |
& +(UICEC(I+1,J ,bi,bj) - UICEC(I ,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1006 |
LSR |
& * _recip_dxG(I,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1007 |
LSR |
dUdy(I,J) = 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
1008 |
LSR |
& ( UICEC(I+1,J+1,bi,bj) - UICEC(I+1,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1009 |
LSR |
& * _recip_dyG(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1010 |
LSR |
& +(UICEC(I ,J+1,bi,bj) - UICEC(I ,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1011 |
LSR |
& * _recip_dyG(I, J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1024 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
1026 |
LSR |
rhsV(I,J,bi,bj)=-DRAGA(I,J,bi,bj)*UICEC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1027 |
LSR |
& +FORCEY(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1038 |
LSR |
& + ( zetaMinusEta(I ,J ,bi,bj) * dUdx(I ,J ) |
|
./pkg/seaice/lsr.F |
1039 |
LSR |
& + zetaMinusEta(I-1,J ,bi,bj) * dUdx(I-1,J ) |
|
./pkg/seaice/lsr.F |
1040 |
LSR |
& - zetaMinusEta(I ,J-1,bi,bj) * dUdx(I ,J-1) |
|
./pkg/seaice/lsr.F |
1041 |
LSR |
& - zetaMinusEta(I-1,J-1,bi,bj) * dUdx(I-1,J-1) |
|
./pkg/seaice/lsr.F |
1042 |
LSR |
& )* 0.5 _d 0 * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1045 |
LSR |
& + ( ETA (I ,J ,bi,bj) * dUdy(I ,J ) |
|
./pkg/seaice/lsr.F |
1046 |
LSR |
& + ETA (I ,J-1,bi,bj) * dUdy(I ,J-1) |
|
./pkg/seaice/lsr.F |
1047 |
LSR |
& - ETA (I-1,J ,bi,bj) * dUdy(I-1,J ) |
|
./pkg/seaice/lsr.F |
1048 |
LSR |
& - ETA (I-1,J-1,bi,bj) * dUdy(I-1,J-1) |
|
./pkg/seaice/lsr.F |
1049 |
LSR |
& )*0.5 _d 0* _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1051 |
LSR |
& +(ETA(I ,J ,bi,bj) + ETA(I ,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
1052 |
LSR |
& -ETA(I-1,J-1,bi,bj) - ETA(I-1,J ,bi,bj)) |
|
./pkg/seaice/lsr.F |
1053 |
LSR |
& * UICEC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1054 |
LSR |
& * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1055 |
LSR |
& * 0.5 _d 0 * _recip_dxV(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
1056 |
LSR |
& +ETAMEAN(I,J,bi,bj) * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1057 |
LSR |
& *(UICEC(I+1,J,bi,bj)-UICEC(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
1058 |
LSR |
& * 0.5 _d 0 * _recip_dxV(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
1060 |
LSR |
& +ETAMEAN(I,J,bi,bj)*TWO * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1061 |
LSR |
& *(UICEC(I+1,J,bi,bj)-UICEC(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
1062 |
LSR |
& * 1. _d 0 /( _dxG(I,J,bi,bj) + _dxG(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
1064 |
LSR |
VRT1(I,J,bi,bj)= 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
1065 |
LSR |
& ETA(I-1,J-1,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1066 |
LSR |
& +ETA(I-1,J ,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1067 |
LSR |
& ) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1068 |
LSR |
VRT2(I,J,bi,bj)= 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
1069 |
LSR |
& ETA(I ,J-1,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1070 |
LSR |
& +ETA(I ,J ,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1071 |
LSR |
& ) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1096 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
1098 |
LSR |
VTMP(I,J,bi,bj)=VICE(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1103 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
1104 |
LSR |
IF(J.EQ.1) THEN |
|
./pkg/seaice/lsr.F |
1105 |
LSR |
AA2= _recip_dyU(I,J,bi,bj) * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1107 |
LSR |
& etaPlusZeta(I-1,J-1,bi,bj) + etaPlusZeta(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
1110 |
LSR |
& +( ZETAMEAN(I,J,bi,bj)-ETAMEAN(I,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
1111 |
LSR |
& * _tanPhiAtV(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
1112 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
1113 |
LSR |
& + ETAMEAN(I,J,bi,bj)*TWO* _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1114 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere ) |
|
./pkg/seaice/lsr.F |
1115 |
LSR |
& *VICE(I,J-1,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1116 |
LSR |
ELSE IF(J.EQ.sNy) THEN |
|
./pkg/seaice/lsr.F |
1117 |
LSR |
AA1= _recip_dyU(I,J,bi,bj) * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1119 |
LSR |
& etaPlusZeta(I-1,J,bi,bj) + etaPlusZeta(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1122 |
LSR |
& -( ZETAMEAN(I,J,bi,bj)-ETAMEAN(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
1123 |
LSR |
& * _tanPhiAtV(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
1124 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
1125 |
LSR |
& - ETAMEAN(I,J,bi,bj)*TWO* _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1126 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere ) |
|
./pkg/seaice/lsr.F |
1127 |
LSR |
& *VICE(I,J+1,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1132 |
LSR |
VRT(J)=rhsV(I,J,bi,bj)+AA3+VRT1(I,J,bi,bj)*VICE(I-1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1133 |
LSR |
& +VRT2(I,J,bi,bj)*VICE(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1134 |
LSR |
VRT(J)=VRT(J)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1137 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
1138 |
LSR |
CVV(J)=CV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1141 |
LSR |
DO J=2,sNy |
|
./pkg/seaice/lsr.F |
1142 |
LSR |
JM=J-1 |
|
./pkg/seaice/lsr.F |
1143 |
LSR |
CVV(J)=CVV(J)/(BV(I,J,bi,bj)-AV(I,J,bi,bj)*CVV(JM)) |
|
./pkg/seaice/lsr.F |
1144 |
LSR |
VRT(J)=(VRT(J)-AV(I,J,bi,bj)*VRT(JM)) |
|
./pkg/seaice/lsr.F |
1145 |
LSR |
& /(BV(I,J,bi,bj)-AV(I,J,bi,bj)*CVV(JM)) |
|
./pkg/seaice/lsr.F |
1147 |
LSR |
DO J=1,sNy-1 |
|
./pkg/seaice/lsr.F |
1148 |
LSR |
J1=sNy-J |
|
./pkg/seaice/lsr.F |
1152 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
1153 |
LSR |
VICE(I,J,bi,bj)=VTMP(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1154 |
LSR |
& +WFAV*(VRT(J)-VTMP(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
1165 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
1167 |
LSR |
UERR=(VICE(I,J,bi,bj)-VTMP(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
1168 |
LSR |
& *UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1202 |
LSR |
DO J=1-OLy,sNy+OLy |
|
./pkg/seaice/lsr.F |
1204 |
LSR |
UICE(I,J,bi,bj)=UICE(I,J,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
1205 |
LSR |
VICE(I,J,bi,bj)=VICE(I,J,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
174 |
LSR |
ucLoc(I,J,bi,bj) = 0.25 _d 0 * ( |
|
./pkg/seaice/lsr.F |
175 |
LSR |
& + uIceC(I+1,J,bi,bj) + uIceC(I+1,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
176 |
LSR |
& + uIceC(I ,J,bi,bj) + uIceC(I, J+1,bi,bj) ) |
|
./pkg/seaice/lsr.F |
177 |
LSR |
vcLoc(I,J,bi,bj) = 0.25 _d 0 * ( |
|
./pkg/seaice/lsr.F |
178 |
LSR |
& + vIceC(I+1,J,bi,bj) + vIceC(I+1,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
179 |
LSR |
& + vIceC(I ,J,bi,bj) + vIceC(I, J+1,bi,bj) ) |
|
./pkg/seaice/lsr.F |
182 |
LSR |
DO J=1-OLy,sNy+OLy-1 |
|
./pkg/seaice/lsr.F |
184 |
LSR |
e11loc = 0.5 _d 0 * _recip_dxF(I,J,bi,bj) * |
|
./pkg/seaice/lsr.F |
185 |
LSR |
& (uIceC(I+1,J+1,bi,bj)+uIceC(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
186 |
LSR |
& -uIceC(I, J+1,bi,bj)-uIceC(I, J,bi,bj)) |
|
./pkg/seaice/lsr.F |
187 |
LSR |
& + vcLoc(I,J,bi,bj) * k2AtC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
188 |
LSR |
e22loc = 0.5 _d 0 * _recip_dyF(I,J,bi,bj) * |
|
./pkg/seaice/lsr.F |
189 |
LSR |
& (vIceC(I+1,J+1,bi,bj)+vIceC(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
190 |
LSR |
& -vIceC(I+1,J, bi,bj)-vIceC(I,J, bi,bj)) |
|
./pkg/seaice/lsr.F |
191 |
LSR |
& + ucLoc(I,J,bi,bj) * k1AtC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
193 |
LSR |
& 0.5 _d 0 * _recip_dyF(I,J,bi,bj) * |
|
./pkg/seaice/lsr.F |
194 |
LSR |
& (uIceC(I+1,J+1,bi,bj)+uIceC(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
195 |
LSR |
& -uIceC(I+1,J, bi,bj)-uIceC(I,J, bi,bj)) |
|
./pkg/seaice/lsr.F |
196 |
LSR |
& + 0.5 _d 0 * _recip_dxF(I,J,bi,bj) * |
|
./pkg/seaice/lsr.F |
197 |
LSR |
& (vIceC(I+1,J+1,bi,bj)+vIceC(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
198 |
LSR |
& -vIceC(I, J+1,bi,bj)-vIceC(I, J,bi,bj)) |
|
./pkg/seaice/lsr.F |
199 |
LSR |
& - vcLoc(I,J,bi,bj)*k1AtC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
200 |
LSR |
& - ucLoc(I,J,bi,bj)*k2AtC(I,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
210 |
LSR |
ZETA(I,J,bi,bj) = 0.5 _d 0*PRESS0(I,J,bi,bj)/DELT2 |
|
./pkg/seaice/lsr.F |
212 |
LSR |
ZETA(I,J,bi,bj) = MIN(ZMAX(I,J,bi,bj),ZETA(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
213 |
LSR |
ZETA(I,J,bi,bj) = MAX(ZMIN(I,J,bi,bj),ZETA(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
215 |
LSR |
ZETA(I,J,bi,bj) = ZETA(I,J,bi,bj)*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
216 |
LSR |
ETA(I,J,bi,bj) = ECM2*ZETA(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
217 |
LSR |
PRESS(I,J,bi,bj) = 2.0 _d 0*ZETA(I,J,bi,bj)*DELT2 |
|
./pkg/seaice/lsr.F |
223 |
LSR |
TEMPVAR=(uIce(I,J,bi,bj)-GWATX(I,J,bi,bj))**2 |
|
./pkg/seaice/lsr.F |
224 |
LSR |
& +(vIce(I,J,bi,bj)-GWATY(I,J,bi,bj))**2 |
|
./pkg/seaice/lsr.F |
225 |
LSR |
IF ( YC(I,J,bi,bj) .LT. ZERO ) THEN |
|
./pkg/seaice/lsr.F |
227 |
LSR |
DWATN(I,J,bi,bj)=QUART |
|
./pkg/seaice/lsr.F |
229 |
LSR |
DWATN(I,J,bi,bj)=SEAICE_waterDrag_south*SQRT(TEMPVAR) |
|
./pkg/seaice/lsr.F |
233 |
LSR |
DWATN(I,J,bi,bj)=QUART |
|
./pkg/seaice/lsr.F |
235 |
LSR |
DWATN(I,J,bi,bj)=SEAICE_waterDrag*SQRT(TEMPVAR) |
|
./pkg/seaice/lsr.F |
239 |
LSR |
DRAGS(I,J,bi,bj)=DWATN(I,J,bi,bj)*COSWAT |
|
./pkg/seaice/lsr.F |
241 |
LSR |
DRAGA(I,J,bi,bj)=DWATN(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
242 |
LSR |
& *SIGN(SINWAT, _fCori(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
243 |
LSR |
& + AMASS(I,J,bi,bj) * _fCoriG(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
245 |
LSR |
FORCEX(I,J,bi,bj)=FORCEX0(I,J,bi,bj)+DWATN(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
246 |
LSR |
& *(COSWAT*GWATX(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
247 |
LSR |
& -SIGN(SINWAT, _fCori(I,J,bi,bj))*GWATY(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
248 |
LSR |
FORCEY(I,J,bi,bj)=FORCEY0(I,J,bi,bj)+DWATN(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
249 |
LSR |
& *(SIGN(SINWAT, _fCori(I,J,bi,bj))*GWATX(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
250 |
LSR |
& +COSWAT*GWATY(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
255 |
LSR |
FORCEX(I,J,bi,bj)=FORCEX(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
256 |
LSR |
& -QUART * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
257 |
LSR |
& *(PRESS(I, J,bi,bj) + PRESS(I, J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
258 |
LSR |
& - PRESS(I-1,J,bi,bj) - PRESS(I-1,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
259 |
LSR |
FORCEY(I,J,bi,bj)=FORCEY(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
260 |
LSR |
& -QUART * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
261 |
LSR |
& *(PRESS(I,J, bi,bj) + PRESS(I-1,J, bi,bj) |
|
./pkg/seaice/lsr.F |
262 |
LSR |
& - PRESS(I,J-1,bi,bj) - PRESS(I-1,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
264 |
LSR |
FORCEX(I,J,bi,bj)=FORCEX(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
265 |
LSR |
& +AMASS(I,J,bi,bj)/SEAICE_deltaTdyn*uIceNm1(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
266 |
LSR |
FORCEY(I,J,bi,bj)=FORCEY(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
267 |
LSR |
& +AMASS(I,J,bi,bj)/SEAICE_deltaTdyn*vIceNm1(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
268 |
LSR |
FORCEX(I,J,bi,bj)=FORCEX(I,J,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
269 |
LSR |
FORCEY(I,J,bi,bj)=FORCEY(I,J,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
281 |
LSR |
etaU(I,J) = 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
282 |
LSR |
& + eta(I,J,bi,bj) + eta(I-1,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
283 |
LSR |
zetaU(I,J)= 0.5 _d 0 *( |
|
./pkg/seaice/lsr.F |
284 |
LSR |
& + zeta(I,J,bi,bj) + zeta(I-1,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
289 |
LSR |
etaV(I,J) = 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
290 |
LSR |
& + eta(I,J,bi,bj) + eta(I,J-1,bi,bj) ) |
|
./pkg/seaice/lsr.F |
291 |
LSR |
zetaV(I,J)= 0.5 _d 0 *( |
|
./pkg/seaice/lsr.F |
292 |
LSR |
& + zeta(I,J,bi,bj) + zeta(I,J-1,bi,bj) ) |
|
./pkg/seaice/lsr.F |
296 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
299 |
LSR |
UXX(I,J) = _dyC(I,J,bi,bj) * (zetaV(I,J)+etaV(I,J)) |
|
./pkg/seaice/lsr.F |
300 |
LSR |
& * _recip_dxG(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
302 |
LSR |
UXM(I,J) = _dyC(I,J,bi,bj) * (zetaV(I,J)-etaV(I,J)) |
|
./pkg/seaice/lsr.F |
303 |
LSR |
& * k1AtV(I,J,bi,bj) * 0.5 _d 0 |
|
./pkg/seaice/lsr.F |
306 |
LSR |
DO J=0,sNy |
|
./pkg/seaice/lsr.F |
309 |
LSR |
UYY(I,J) = _dxC(I,J,bi,bj) * etaU(I,J) |
|
./pkg/seaice/lsr.F |
310 |
LSR |
& * _recip_dyG(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
312 |
LSR |
UYM(I,J) = _dxC(I,J,bi,bj) * etaU(I,J) |
|
./pkg/seaice/lsr.F |
313 |
LSR |
& * k2AtU(I,J,bi,bj) * 0.5 _d 0 |
|
./pkg/seaice/lsr.F |
316 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
319 |
LSR |
VXX(I,J) = _dyC(I,J,bi,bj) * etaV(I,J) |
|
./pkg/seaice/lsr.F |
320 |
LSR |
& * _recip_dxG(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
322 |
LSR |
VXM(I,J) = _dyC(I,J,bi,bj) * etaV(I,J) |
|
./pkg/seaice/lsr.F |
323 |
LSR |
& * k1AtV(I,J,bi,bj) * 0.5 _d 0 |
|
./pkg/seaice/lsr.F |
326 |
LSR |
DO J=0,sNy |
|
./pkg/seaice/lsr.F |
329 |
LSR |
VYY(I,J) = _dxC(I,J,bi,bj) * (zetaU(I,J)+etaU(I,J)) |
|
./pkg/seaice/lsr.F |
330 |
LSR |
& * _recip_dyG(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
332 |
LSR |
VYM(I,J) = _dxC(I,J,bi,bj) * (zetaU(I,J)-etaU(I,J)) |
|
./pkg/seaice/lsr.F |
333 |
LSR |
& * k2AtU(I,J,bi,bj) * 0.5 _d 0 |
|
./pkg/seaice/lsr.F |
341 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
344 |
LSR |
AU(I,J,bi,bj)= ( - UXX(I-1,J) + UXM(I-1,J) ) |
|
./pkg/seaice/lsr.F |
345 |
LSR |
& * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
347 |
LSR |
CU(I,J,bi,bj)= ( - UXX(I ,J) - UXM(I ,J) ) |
|
./pkg/seaice/lsr.F |
348 |
LSR |
& * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
350 |
LSR |
BU(I,J,bi,bj)=(ONE - UVM(I,J,bi,bj)) + |
|
./pkg/seaice/lsr.F |
351 |
LSR |
& ( UXX(I-1,J) + UXX(I,J) + UYY(I,J-1) + UYY(I,J) |
|
./pkg/seaice/lsr.F |
352 |
LSR |
& + UXM(I-1,J) - UXM(I,J) - UYM(I,J-1) + UYM(I,J) |
|
./pkg/seaice/lsr.F |
353 |
LSR |
& ) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
355 |
LSR |
uRt1(I,J,bi,bj)= UYY(I,J-1) + UYM(I,J-1) |
|
./pkg/seaice/lsr.F |
357 |
LSR |
uRt2(I,J,bi,bj)= UYY(I,J ) - UYM(I,J ) |
|
./pkg/seaice/lsr.F |
362 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
364 |
LSR |
AU(I,J,bi,bj) = AU(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
365 |
LSR |
CU(I,J,bi,bj) = CU(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
368 |
LSR |
BU(I,J,bi,bj) = BU(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
369 |
LSR |
& + UVM(I,J,bi,bj) * |
|
./pkg/seaice/lsr.F |
370 |
LSR |
& ( AMASS(I,J,bi,bj)/SEAICE_deltaTdyn |
|
./pkg/seaice/lsr.F |
371 |
LSR |
& + DRAGS(I,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
372 |
LSR |
uRt1(I,J,bi,bj) = uRt1(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
373 |
LSR |
uRt2(I,J,bi,bj) = uRt2(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
381 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
384 |
LSR |
AV(I,J,bi,bj)=( - VYY(I,J-1) + VYM(I,J-1) |
|
./pkg/seaice/lsr.F |
385 |
LSR |
& ) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
387 |
LSR |
CV(I,J,bi,bj)=( - VYY(I,J ) - VYM(I,J ) |
|
./pkg/seaice/lsr.F |
388 |
LSR |
& ) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
390 |
LSR |
BV(I,J,bi,bj)= (ONE - UVM(I,J,bi,bj)) + |
|
./pkg/seaice/lsr.F |
391 |
LSR |
& ( VXX(I,J) + VXX(I-1,J) + VYY(I,J) + VYY(I,J-1) |
|
./pkg/seaice/lsr.F |
392 |
LSR |
& + VXM(I,J) - VXM(I-1,J) - VYM(I,J) + VYM(I,J-1) |
|
./pkg/seaice/lsr.F |
393 |
LSR |
& ) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
395 |
LSR |
vRt1(I,J,bi,bj) = VXX(I-1,J) + VXM(I-1,J) |
|
./pkg/seaice/lsr.F |
397 |
LSR |
vRt2(I,J,bi,bj) = VXX(I ,J) - VXM(I ,J) |
|
./pkg/seaice/lsr.F |
402 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
404 |
LSR |
AV(I,J,bi,bj) = AV(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
405 |
LSR |
CV(I,J,bi,bj) = CV(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
408 |
LSR |
BV(I,J,bi,bj) = BV(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
409 |
LSR |
& + UVM(I,J,bi,bj) * |
|
./pkg/seaice/lsr.F |
410 |
LSR |
& ( AMASS(I,J,bi,bj)/SEAICE_deltaTdyn |
|
./pkg/seaice/lsr.F |
411 |
LSR |
& + DRAGS(I,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
412 |
LSR |
vRt1(I,J,bi,bj) = vRt1(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
413 |
LSR |
vRt2(I,J,bi,bj) = vRt2(I,J,bi,bj) * recip_rAz(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
421 |
LSR |
sig11(I,J) = |
|
./pkg/seaice/lsr.F |
422 |
LSR |
& (zetaV(I,J)-etaV(I,J)) |
|
./pkg/seaice/lsr.F |
423 |
LSR |
& * (vcLoc(I,J,bi,bj)-vcLoc(I,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
424 |
LSR |
& * _recip_dyC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
425 |
LSR |
& + (zetaV(I,J)+etaV(I,J))*k2atV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
426 |
LSR |
& * 0.5 _d 0 * (vIceC(I,J,bi,bj)+vIceC(I+1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
427 |
LSR |
& - 0.5 _d 0 * (PRESS(I,J,bi,bj)+PRESS(I,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
428 |
LSR |
sig12(I,J) = etaV(I,J) |
|
./pkg/seaice/lsr.F |
429 |
LSR |
& * (ucLoc(I,J,bi,bj)-ucLoc(I,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
430 |
LSR |
& * _recip_dyC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
431 |
LSR |
& - etaV(I,J) * k2AtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
432 |
LSR |
& * 0.5 _d 0 * (uIceC(I,J,bi,bj)+uIceC(I+1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
438 |
LSR |
sig22(I,J) = |
|
./pkg/seaice/lsr.F |
439 |
LSR |
& (zetaU(I,J)-etaU(I,J)) |
|
./pkg/seaice/lsr.F |
440 |
LSR |
& * (ucLoc(I,J,bi,bj)-ucLoc(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
441 |
LSR |
& * _recip_dxC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
442 |
LSR |
& + (zetaU(I,J)+etaU(I,J))*k2atU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
443 |
LSR |
& * 0.5 _d 0 * (uIceC(I,J,bi,bj)+uIceC(I,J+1,bi,bj)) |
|
./pkg/seaice/lsr.F |
444 |
LSR |
& - 0.5 _d 0 * (PRESS(I,J,bi,bj)+PRESS(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
445 |
LSR |
sig21(I,J) = etaU(I,J) |
|
./pkg/seaice/lsr.F |
446 |
LSR |
& * (vcLoc(I,J,bi,bj)-vcLoc(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
447 |
LSR |
& * _recip_dxC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
448 |
LSR |
& - etaU(I,J) * k1AtU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
449 |
LSR |
& * 0.5 _d 0 * (vIceC(I,J,bi,bj)+vIceC(I,J+1,bi,bj)) |
|
./pkg/seaice/lsr.F |
454 |
LSR |
rhsU(I,J,bi,bj) = DRAGA(I,J,bi,bj)*vIceC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
455 |
LSR |
& +FORCEX(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
456 |
LSR |
& + ( _dyC(I, J, bi,bj) * sig11(I, J ) |
|
./pkg/seaice/lsr.F |
457 |
LSR |
& - _dyC(I-1,J, bi,bj) * sig11(I-1,J ) |
|
./pkg/seaice/lsr.F |
458 |
LSR |
& + _dxC(I, J, bi,bj) * sig21(I, J ) |
|
./pkg/seaice/lsr.F |
459 |
LSR |
& - _dxC(I, J-1,bi,bj) * sig21(I, J-1) ) |
|
./pkg/seaice/lsr.F |
460 |
LSR |
& * recip_rAz(I,J,bi,bj) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
461 |
LSR |
rhsV(I,J,bi,bj) = - DRAGA(I,J,bi,bj)*uIceC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
462 |
LSR |
& +FORCEY(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
463 |
LSR |
& + ( _dyC(I, J, bi,bj) * sig12(I, J ) |
|
./pkg/seaice/lsr.F |
464 |
LSR |
& - _dyC(I-1,J, bi,bj) * sig12(I-1,J ) |
|
./pkg/seaice/lsr.F |
465 |
LSR |
& + _dxC(I, J, bi,bj) * sig22(I, J ) |
|
./pkg/seaice/lsr.F |
466 |
LSR |
& - _dxC(I, J-1,bi,bj) * sig22(I, J-1) ) |
|
./pkg/seaice/lsr.F |
467 |
LSR |
& * recip_rAz(I,J,bi,bj) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
496 |
LSR |
uTmp(I,J,bi,bj)=uIce(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
502 |
LSR |
vTmp(I,J,bi,bj)=vIce(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
508 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
511 |
LSR |
IF (I.EQ.1) AA3 = AA3 - AU(I,J,bi,bj)*uIce(I-1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
512 |
LSR |
IF (I.EQ.sNx) AA3 = AA3 - CU(I,J,bi,bj)*uIce(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
514 |
LSR |
URT(I,J)=rhsU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
517 |
LSR |
& + uRt1(I,J,bi,bj)*uTmp(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
518 |
LSR |
& + uRt2(I,J,bi,bj)*uTmp(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
520 |
LSR |
& + uRt1(I,J,bi,bj)*uIce(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
521 |
LSR |
& + uRt2(I,J,bi,bj)*uIce(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
523 |
LSR |
URT(I,J)=URT(I,J)* UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
527 |
LSR |
CUU(I,J)=CU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
529 |
LSR |
CUU(1,J)=CUU(1,J)/BU(1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
530 |
LSR |
URT(1,J)=URT(1,J)/BU(1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
536 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
541 |
LSR |
CUU(I,J)=CUU(I,J)/(BU(I,J,bi,bj)-AU(I,J,bi,bj)*CUU(IM,J)) |
|
./pkg/seaice/lsr.F |
542 |
LSR |
URT(I,J)=(URT(I,J)-AU(I,J,bi,bj)*URT(IM,J)) |
|
./pkg/seaice/lsr.F |
543 |
LSR |
& /(BU(I,J,bi,bj)-AU(I,J,bi,bj)*CUU(IM,J)) |
|
./pkg/seaice/lsr.F |
548 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
553 |
LSR |
URT(J1,J)=URT(J1,J)-CUU(J1,J)*URT(J2,J) |
|
./pkg/seaice/lsr.F |
556 |
LSR |
uIce(I,J,bi,bj)=uTmp(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
557 |
LSR |
& +WFAU*(URT(I,J)-uTmp(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
566 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
568 |
LSR |
IF (J.EQ.1) AA3 = AA3 - AV(I,J,bi,bj)*vIce(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
569 |
LSR |
IF (J.EQ.sNy) AA3 = AA3 - CV(I,J,bi,bj)*vIce(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
571 |
LSR |
VRT(I,J)=rhsV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
574 |
LSR |
& + vRt1(I,J,bi,bj)*vTmp(I-1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
575 |
LSR |
& + vRt2(I,J,bi,bj)*vTmp(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
577 |
LSR |
& + vRt1(I,J,bi,bj)*vIce(I-1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
578 |
LSR |
& + vRt2(I,J,bi,bj)*vIce(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
580 |
LSR |
VRT(I,J)=VRT(I,J)* UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
583 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
584 |
LSR |
CVV(I,J)=CV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
588 |
LSR |
DO J=2,sNy |
|
./pkg/seaice/lsr.F |
589 |
LSR |
JM=J-1 |
|
./pkg/seaice/lsr.F |
590 |
LSR |
CVV(I,J)=CVV(I,J)/(BV(I,J,bi,bj)-AV(I,J,bi,bj)*CVV(I,JM)) |
|
./pkg/seaice/lsr.F |
591 |
LSR |
VRT(I,J)=(VRT(I,J)-AV(I,J,bi,bj)*VRT(I,JM)) |
|
./pkg/seaice/lsr.F |
592 |
LSR |
& /(BV(I,J,bi,bj)-AV(I,J,bi,bj)*CVV(I,JM)) |
|
./pkg/seaice/lsr.F |
594 |
LSR |
DO J=1,sNy-1 |
|
./pkg/seaice/lsr.F |
595 |
LSR |
J1=sNy-J |
|
./pkg/seaice/lsr.F |
599 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
600 |
LSR |
vIce(I,J,bi,bj)=vTmp(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
601 |
LSR |
& +WFAV*(VRT(I,J)-vTmp(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
615 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
617 |
LSR |
UERR=(uIce(I,J,bi,bj)-uTmp(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
618 |
LSR |
& * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
640 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
642 |
LSR |
UERR=(vIce(I,J,bi,bj)-vTmp(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
643 |
LSR |
& * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
704 |
LSR |
etaPlusZeta(I,J,bi,bj) = ETA(I,J,bi,bj)+ZETA(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
705 |
LSR |
zetaMinusEta(I,J,bi,bj) = ZETA(I,J,bi,bj)-ETA(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
710 |
LSR |
ETAMEAN(I,J,bi,bj) =QUART*( |
|
./pkg/seaice/lsr.F |
711 |
LSR |
& ETA(I,J-1,bi,bj) + ETA(I-1,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
712 |
LSR |
& +ETA(I,J ,bi,bj) + ETA(I-1,J ,bi,bj)) |
|
./pkg/seaice/lsr.F |
713 |
LSR |
ZETAMEAN(I,J,bi,bj)=QUART*( |
|
./pkg/seaice/lsr.F |
714 |
LSR |
& ZETA(I,J-1,bi,bj) + ZETA(I-1,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
715 |
LSR |
& +ZETA(I,J ,bi,bj) + ZETA(I-1,J ,bi,bj)) |
|
./pkg/seaice/lsr.F |
724 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
726 |
LSR |
AA1=( etaPlusZeta(I ,J-1,bi,bj) * _recip_dxF(I ,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
727 |
LSR |
& +etaPlusZeta(I ,J ,bi,bj) * _recip_dxF(I ,J ,bi,bj) |
|
./pkg/seaice/lsr.F |
728 |
LSR |
& )*0.5 _d 0 * _recip_dxV(I,J,bi,bj) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
729 |
LSR |
AA2=( etaPlusZeta(I-1,J-1,bi,bj) * _recip_dxF(I-1,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
730 |
LSR |
& +etaPlusZeta(I-1,J ,bi,bj) * _recip_dxF(I-1,J ,bi,bj) |
|
./pkg/seaice/lsr.F |
731 |
LSR |
& )*0.5 _d 0 * _recip_dxV(I,J,bi,bj) * UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
732 |
LSR |
AA3= 0.5 _d 0 *(ETA(I-1,J ,bi,bj)+ETA(I,J ,bi,bj)) |
|
./pkg/seaice/lsr.F |
733 |
LSR |
AA4= 0.5 _d 0 *(ETA(I-1,J-1,bi,bj)+ETA(I,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
734 |
LSR |
AA5= -(AA3-AA4) * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
735 |
LSR |
& * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
736 |
LSR |
AA6=TWO*ETAMEAN(I,J,bi,bj) *recip_rSphere*recip_rSphere |
|
./pkg/seaice/lsr.F |
737 |
LSR |
& * _tanPhiAtV(I,J,bi,bj) * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
738 |
LSR |
AU(I,J,bi,bj)=-AA2 |
|
./pkg/seaice/lsr.F |
739 |
LSR |
CU(I,J,bi,bj)=-AA1 |
|
./pkg/seaice/lsr.F |
740 |
LSR |
BU(I,J,bi,bj)=(ONE-UVM(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
741 |
LSR |
& - AU(I,J,bi,bj) - CU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
742 |
LSR |
& + ((AA3+AA4)*_recip_dyU(I,J,bi,bj)*_recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
744 |
LSR |
& + AMASS(I,J,bi,bj)/SEAICE_deltaTdyn |
|
./pkg/seaice/lsr.F |
745 |
LSR |
& + DRAGS(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
746 |
LSR |
& )*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
750 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
751 |
LSR |
AU(1,J,bi,bj)=ZERO |
|
./pkg/seaice/lsr.F |
752 |
LSR |
CU(sNx,J,bi,bj)=ZERO |
|
./pkg/seaice/lsr.F |
753 |
LSR |
CU(1,J,bi,bj)=CU(1,J,bi,bj)/BU(1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
757 |
LSR |
DO J=1-OLy,sNy+OLy-1 |
|
./pkg/seaice/lsr.F |
759 |
LSR |
dVdy(I,J) = 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
760 |
LSR |
& ( VICEC(I+1,J+1,bi,bj) - VICEC(I+1,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
761 |
LSR |
& * _recip_dyG(I+1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
762 |
LSR |
& +(VICEC(I ,J+1,bi,bj) - VICEC(I ,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
763 |
LSR |
& * _recip_dyG(I, J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
764 |
LSR |
dVdx(I,J) = 0.5 _d 0 * ( |
|
./pkg/seaice/lsr.F |
765 |
LSR |
& ( VICEC(I+1,J+1,bi,bj) - VICEC(I ,J+1,bi,bj) ) |
|
./pkg/seaice/lsr.F |
766 |
LSR |
& * _recip_dxG(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
767 |
LSR |
& +(VICEC(I+1,J ,bi,bj) - VICEC(I ,J ,bi,bj) ) |
|
./pkg/seaice/lsr.F |
768 |
LSR |
& * _recip_dxG(I,J, bi,bj) ) |
|
./pkg/seaice/lsr.F |
781 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
783 |
LSR |
rhsU(I,J,bi,bj)=DRAGA(I,J,bi,bj)*VICEC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
784 |
LSR |
& +FORCEX(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
794 |
LSR |
& + ( zetaMinusEta(I ,J ,bi,bj) * dVdy(I ,J ) |
|
./pkg/seaice/lsr.F |
795 |
LSR |
& + zetaMinusEta(I ,J-1,bi,bj) * dVdy(I ,J-1) |
|
./pkg/seaice/lsr.F |
796 |
LSR |
& - zetaMinusEta(I-1,J ,bi,bj) * dVdy(I-1,J ) |
|
./pkg/seaice/lsr.F |
797 |
LSR |
& - zetaMinusEta(I-1,J-1,bi,bj) * dVdy(I-1,J-1) |
|
./pkg/seaice/lsr.F |
798 |
LSR |
& )* 0.5 _d 0 * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
801 |
LSR |
& + ( ETA (I ,J ,bi,bj) * dVdx(I ,J ) |
|
./pkg/seaice/lsr.F |
802 |
LSR |
& + ETA (I-1,J ,bi,bj) * dVdx(I-1,J ) |
|
./pkg/seaice/lsr.F |
803 |
LSR |
& - ETA (I ,J-1,bi,bj) * dVdx(I ,J-1) |
|
./pkg/seaice/lsr.F |
804 |
LSR |
& - ETA (I-1,J-1,bi,bj) * dVdx(I-1,J-1) |
|
./pkg/seaice/lsr.F |
805 |
LSR |
& ) * 0.5 _d 0 * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
807 |
LSR |
& -(etaPlusZeta(I ,J ,bi,bj)+etaPlusZeta(I ,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
808 |
LSR |
& -etaPlusZeta(I-1,J-1,bi,bj)-etaPlusZeta(I-1,J ,bi,bj)) |
|
./pkg/seaice/lsr.F |
809 |
LSR |
& * VICEC(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
810 |
LSR |
& * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
811 |
LSR |
& * 0.5 _d 0 * _recip_dxV(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
813 |
LSR |
& -(ETAMEAN(I,J,bi,bj)+ZETAMEAN(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
814 |
LSR |
& *(VICEC(I+1,J,bi,bj) - VICEC(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
815 |
LSR |
& * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
816 |
LSR |
& * 1.0 _d 0 /( _dxG(I,J,bi,bj) + _dxG(I-1,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
819 |
LSR |
& -ETAMEAN(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
820 |
LSR |
& *(VICEC(I+1,J,bi,bj) - VICEC(I-1,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
821 |
LSR |
& *TWO* _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
822 |
LSR |
& * 1.0 _d 0 /( _dxG(I,J,bi,bj) + _dxG(I-1,J,bi,bj) ) |
|
./pkg/seaice/lsr.F |
825 |
LSR |
URT1(I,J,bi,bj)= |
|
./pkg/seaice/lsr.F |
826 |
LSR |
& 0.5 _d 0 * (ETA(I-1,J-1,bi,bj)+ETA(I,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
827 |
LSR |
& * _recip_dyU(I,J,bi,bj) * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
828 |
LSR |
& - ETAMEAN(I,J,bi,bj) * _tanPhiAtV(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
829 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
830 |
LSR |
& + TWO*ETAMEAN(I,J,bi,bj) * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
831 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
832 |
LSR |
URT2(I,J,bi,bj)= |
|
./pkg/seaice/lsr.F |
833 |
LSR |
& 0.5 _d 0 * (ETA(I-1,J,bi,bj)+ETA(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
834 |
LSR |
& * _recip_dyU(I,J,bi,bj) * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
835 |
LSR |
& + ETAMEAN(I,J,bi,bj) * _tanPhiAtV(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
836 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
837 |
LSR |
& - TWO*ETAMEAN(I,J,bi,bj) * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
838 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
861 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
863 |
LSR |
UTMP(I,J,bi,bj)=UICE(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
867 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
870 |
LSR |
AA2=(etaPlusZeta(I-1,J-1,bi,bj) * _recip_dxF(I-1,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
871 |
LSR |
& +etaPlusZeta(I-1,J ,bi,bj) * _recip_dxF(I-1,J ,bi,bj) |
|
./pkg/seaice/lsr.F |
872 |
LSR |
& )*0.5 _d 0 * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
873 |
LSR |
AA3=AA2*UICE(I-1,J,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
875 |
LSR |
AA1=(etaPlusZeta(I ,J-1,bi,bj) * _recip_dxF(I ,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
876 |
LSR |
& +etaPlusZeta(I ,J ,bi,bj) * _recip_dxF(I ,J ,bi,bj) |
|
./pkg/seaice/lsr.F |
877 |
LSR |
& )*0.5 _d 0 * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
878 |
LSR |
AA3=AA1*UICE(I+1,J,bi,bj)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
882 |
LSR |
URT(I)=rhsU(I,J,bi,bj)+AA3 |
|
./pkg/seaice/lsr.F |
883 |
LSR |
& +URT1(I,J,bi,bj)*UICE(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
884 |
LSR |
& +URT2(I,J,bi,bj)*UICE(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
885 |
LSR |
URT(I)=URT(I)*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
889 |
LSR |
CUU(I)=CU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
891 |
LSR |
URT(1)=URT(1)/BU(1,J,bi,bj) |
|
./pkg/seaice/lsr.F |
894 |
LSR |
CUU(I)=CUU(I)/(BU(I,J,bi,bj)-AU(I,J,bi,bj)*CUU(IM)) |
|
./pkg/seaice/lsr.F |
895 |
LSR |
URT(I)=(URT(I)-AU(I,J,bi,bj)*URT(IM)) |
|
./pkg/seaice/lsr.F |
896 |
LSR |
& /(BU(I,J,bi,bj)-AU(I,J,bi,bj)*CUU(IM)) |
|
./pkg/seaice/lsr.F |
904 |
LSR |
UICE(I,J,bi,bj)=UTMP(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
905 |
LSR |
& +WFAU*(URT(I)-UTMP(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
917 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
919 |
LSR |
UERR=(UICE(I,J,bi,bj)-UTMP(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
920 |
LSR |
& *UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
951 |
LSR |
DO J=1,sNy |
|
./pkg/seaice/lsr.F |
953 |
LSR |
AA1=0.5 _d 0 * _recip_dyU(I,J,bi,bj) * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
954 |
LSR |
& * (etaPlusZeta(I-1,J ,bi,bj) + etaPlusZeta(I,J ,bi,bj)) |
|
./pkg/seaice/lsr.F |
955 |
LSR |
AA2=0.5 _d 0 * _recip_dyU(I,J,bi,bj) * _recip_dyU(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
956 |
LSR |
& * (etaPlusZeta(I-1,J-1,bi,bj) + etaPlusZeta(I,J-1,bi,bj)) |
|
./pkg/seaice/lsr.F |
957 |
LSR |
AA3= (ETA(I ,J-1,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
958 |
LSR |
& +ETA(I ,J ,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
959 |
LSR |
& )* 0.5 _d 0 * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
960 |
LSR |
AA4= (ETA(I-1,J-1,bi,bj)+ETA(I-1,J,bi,bj))*0.5 _d 0 |
|
./pkg/seaice/lsr.F |
961 |
LSR |
& *_recip_dxV(I,J,bi,bj) * _recip_dxV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
962 |
LSR |
AA5=(zetaMinusEta(I-1,J ,bi,bj) + zetaMinusEta(I,J ,bi,bj) |
|
./pkg/seaice/lsr.F |
963 |
LSR |
& -zetaMinusEta(I-1,J-1,bi,bj) - zetaMinusEta(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
964 |
LSR |
& )* _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
965 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
967 |
LSR |
AA6=TWO*ETAMEAN(I,J,bi,bj) * recip_rSphere*recip_rSphere |
|
./pkg/seaice/lsr.F |
968 |
LSR |
& * _tanPhiAtV(I,J,bi,bj) * _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
970 |
LSR |
AV(I,J,bi,bj)=( |
|
./pkg/seaice/lsr.F |
972 |
LSR |
& - (ZETAMEAN(I,J,bi,bj)-ETAMEAN(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
973 |
LSR |
& * _tanPhiAtV(I,J-1,bi,bj) |
|
./pkg/seaice/lsr.F |
974 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
975 |
LSR |
& -ETAMEAN(I,J,bi,bj)*TWO* _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
976 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
977 |
LSR |
& )*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
978 |
LSR |
CV(I,J,bi,bj)=( |
|
./pkg/seaice/lsr.F |
980 |
LSR |
& +(ZETAMEAN(I,J,bi,bj)-ETAMEAN(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
981 |
LSR |
& * _tanPhiAtV(I,J+1,bi,bj) |
|
./pkg/seaice/lsr.F |
982 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
983 |
LSR |
& +ETAMEAN(I,J,bi,bj)*TWO* _tanPhiAtV(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
984 |
LSR |
& * 0.5 _d 0 * _recip_dyU(I,J,bi,bj)*recip_rSphere |
|
./pkg/seaice/lsr.F |
985 |
LSR |
& )*UVM(I,J,bi,bj) |
|
./pkg/seaice/lsr.F |
986 |
LSR |
BV(I,J,bi,bj)= (ONE-UVM(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
988 |
LSR |
& +AMASS(I,J,bi,bj)/SEAICE_deltaTdyn+DRAGS(I,J,bi,bj)) |
|
./pkg/seaice/lsr.F |
989 |
LSR |
& *UVM(I,J,bi,bj) |
|
./pkg/seaice/ostres.F |
55 |
OSTRES |
fu(I,J,bi,bj)=HALF |
|
./pkg/seaice/ostres.F |
56 |
OSTRES |
& *(WINDX(I,J+1,bi,bj)+WINDX(I,J,bi,bj)) |
|
./pkg/seaice/ostres.F |
57 |
OSTRES |
fv(I,J,bi,bj)=HALF |
|
./pkg/seaice/ostres.F |
58 |
OSTRES |
& *(WINDY(I+1,J,bi,bj)+WINDY(I,J,bi,bj)) |
|
./pkg/seaice/ostres.F |
72 |
OSTRES |
fuIce=QUART*( DWATN(I,J,bi,bj)+DWATN(I,J+1,bi,bj) )*( |
|
./pkg/seaice/ostres.F |
74 |
OSTRES |
& ( UICE(I,J, bi,bj)-GWATX(I,J, bi,bj) |
|
./pkg/seaice/ostres.F |
75 |
OSTRES |
& + UICE(I,J+1,bi,bj)-GWATX(I,J+1,bi,bj) ) |
|
./pkg/seaice/ostres.F |
76 |
OSTRES |
& -SIGN(SINWAT,COR_ICE(I,J,bi,bj)) * |
|
./pkg/seaice/ostres.F |
77 |
OSTRES |
& ( VICE(I,J, bi,bj)-GWATY(I,J, bi,bj) |
|
./pkg/seaice/ostres.F |
78 |
OSTRES |
& + VICE(I,J+1,bi,bj)-GWATY(I,J+1,bi,bj) ) |
|
./pkg/seaice/ostres.F |
80 |
OSTRES |
fvIce=QUART*( DWATN(I,J,bi,bj)+DWATN(I+1,J,bi,bj) )*( |
|
./pkg/seaice/ostres.F |
81 |
OSTRES |
& SIGN(SINWAT,COR_ICE(I,J,bi,bj)) * |
|
./pkg/seaice/ostres.F |
82 |
OSTRES |
& ( UICE(I, J,bi,bj)-GWATX(I, J,bi,bj) |
|
./pkg/seaice/ostres.F |
83 |
OSTRES |
& + UICE(I+1,J,bi,bj)-GWATX(I+1,J,bi,bj) ) |
|
./pkg/seaice/ostres.F |
85 |
OSTRES |
& ( VICE(I, J,bi,bj)-GWATY(I, J,bi,bj) |
|
./pkg/seaice/ostres.F |
86 |
OSTRES |
& + VICE(I+1,J,bi,bj)-GWATY(I+1,J,bi,bj) ) |
|
./pkg/seaice/ostres.F |
88 |
OSTRES |
fu(I,J,bi,bj)=(ONE-AREA(I,J,bi,bj))*fu(I,J,bi,bj)+ |
|
./pkg/seaice/ostres.F |
89 |
OSTRES |
& AREA(I,J,bi,bj)*fuIce |
|
./pkg/seaice/ostres.F |
90 |
OSTRES |
fv(I,J,bi,bj)=(ONE-AREA(I,J,bi,bj))*fv(I,J,bi,bj)+ |
|
./pkg/seaice/ostres.F |
91 |
OSTRES |
& AREA(I,J,bi,bj)*fvIce |
|
./pkg/seaice/seaice_advdiff.F |
435 |
SEAICE_ADVDIFF |
diagArray(I,J,2+(iTr-1)*5) = SItrExt(i,j,bi,bj) |
|
./pkg/seaice/seaice_advdiff.F |
467 |
SEAICE_ADVDIFF |
if (HEFF(I,J,bi,bj).GE.siEps) then |
|
./pkg/seaice/seaice_advdiff.F |
468 |
SEAICE_ADVDIFF |
SItracer(i,j,bi,bj,iTr)=SItrExt(i,j,bi,bj)/HEFF(I,J,bi,bj) |
|
./pkg/seaice/seaice_advdiff.F |
483 |
SEAICE_ADVDIFF |
& +tmpscal2*HEFF(I,J,bi,bj) |
|
./pkg/seaice/seaice_advdiff.F |
486 |
SEAICE_ADVDIFF |
if (HEFF(I,J,bi,bj).GE.siEps) then |
|
./pkg/seaice/seaice_advdiff.F |
490 |
SEAICE_ADVDIFF |
& +HEFF(I,J,bi,bj)*tmpscal1 |
|
./pkg/seaice/seaice_advdiff.F |
493 |
SEAICE_ADVDIFF |
diagArray(I,J,1+(iTr-1)*5)= - SItrBucket(i,j,bi,bj,iTr) |
|
./pkg/seaice/seaice_advdiff.F |
494 |
SEAICE_ADVDIFF |
& *HEFFM(I,J,bi,bj)/SEAICE_deltaTtherm*SEAICE_rhoIce |
|
./pkg/seaice/seaice_advdiff.F |
495 |
SEAICE_ADVDIFF |
tmpscal1= ( HEFF(I,J,bi,bj)*SItracer(i,j,bi,bj,iTr) |
|
./pkg/seaice/seaice_advdiff.F |
496 |
SEAICE_ADVDIFF |
& + SItrBucket(i,j,bi,bj,iTr) )*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/seaice_advdiff.F |
497 |
SEAICE_ADVDIFF |
diagArray(I,J,2+(iTr-1)*5)= tmpscal1-diagArray(I,J,2+(iTr-1)*5) |
|
./pkg/seaice/seaice_advdiff.F |
498 |
SEAICE_ADVDIFF |
diagArray(I,J,3+(iTr-1)*5)=HEFFM(i,j,bi,bj) * |
|
./pkg/seaice/seaice_advdiff.F |
507 |
SEAICE_ADVDIFF |
if (AREA(I,J,bi,bj).GE.SEAICE_area_floor) then |
|
./pkg/seaice/seaice_advdiff.F |
508 |
SEAICE_ADVDIFF |
SItracer(i,j,bi,bj,iTr)=SItrExt(i,j,bi,bj)/AREA(I,J,bi,bj) |
|
./pkg/seaice/seaice_advdiff.F |
521 |
SEAICE_ADVDIFF |
if (AREA(I,J,bi,bj).GE.SEAICE_area_floor) then |
|
./pkg/seaice/seaice_advdiff.F |
526 |
SEAICE_ADVDIFF |
diagArray(I,J,1+(iTr-1)*5)= 0. _d 0 |
|
./pkg/seaice/seaice_advdiff.F |
527 |
SEAICE_ADVDIFF |
diagArray(I,J,2+(iTr-1)*5)= - diagArray(I,J,2+(iTr-1)*5) |
|
./pkg/seaice/seaice_advdiff.F |
528 |
SEAICE_ADVDIFF |
& + AREA(I,J,bi,bj)*SItracer(i,j,bi,bj,iTr)*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/seaice_advdiff.F |
529 |
SEAICE_ADVDIFF |
diagArray(I,J,3+(iTr-1)*5)=HEFFM(i,j,bi,bj) * |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
104 |
SEAICE_BOTTOMDRAG_COEFFS |
IF ( AREA(I,J,bi,bj) .GT. 0.01 _d 0 ) THEN |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
114 |
SEAICE_BOTTOMDRAG_COEFFS |
& + uIce(i+1,J,bi,bj)*maskInW(i+1,j,bi,bj) )**2 |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
118 |
SEAICE_BOTTOMDRAG_COEFFS |
tmpFld(I,J) = SEAICEbasalDragK2 / SQRT(tmp + u0sq) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
127 |
SEAICE_BOTTOMDRAG_COEFFS |
IF ( AREA(I,J,bi,bj) .GT. 0.01 _d 0 ) THEN |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
133 |
SEAICE_BOTTOMDRAG_COEFFS |
hActual = HEFFITD(I,J,k,bi,bj) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
135 |
SEAICE_BOTTOMDRAG_COEFFS |
hCrit = ABS(R_low(I,J,bi,bj))*AREAITD(I,J,k,bi,bj) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
137 |
SEAICE_BOTTOMDRAG_COEFFS |
hActual = HEFF(I,J,bi,bj) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
138 |
SEAICE_BOTTOMDRAG_COEFFS |
hCrit = ABS(R_low(I,J,bi,bj))*AREA(I,J,bi,bj)*recip_k1 |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
148 |
SEAICE_BOTTOMDRAG_COEFFS |
CbotC(I,J,bi,bj) = CbotC(I,J,bi,bj) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
149 |
SEAICE_BOTTOMDRAG_COEFFS |
& + tmpFld(I,J) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
152 |
SEAICE_BOTTOMDRAG_COEFFS |
& *(SEAICE_area_max - AREA(I,J,bi,bj)) ) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
153 |
SEAICE_BOTTOMDRAG_COEFFS |
& * maskC(I,J,kSrf,bi,bj) |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
98 |
SEAICE_BOTTOMDRAG_COEFFS |
CbotC(I,J,bi,bj) = 0. _d 0 |
|
./pkg/seaice/seaice_bottomdrag_coeffs.F |
99 |
SEAICE_BOTTOMDRAG_COEFFS |
tmpFld(I,J) = 0. _d 0 |
|
./pkg/seaice/seaice_budget_ocean.F |
106 |
SEAICE_BUDGET_OCEAN |
DO J=1,sNy |
|
./pkg/seaice/seaice_budget_ocean.F |
109 |
SEAICE_BUDGET_OCEAN |
netHeatFlux(I,J) = Qnet(I,J,bi,bj) |
|
./pkg/seaice/seaice_budget_ocean.F |
110 |
SEAICE_BUDGET_OCEAN |
SWHeatFlux (I,J) = Qsw(I,J,bi,bj) |
|
./pkg/seaice/seaice_budget_ocean.F |
120 |
SEAICE_BUDGET_OCEAN |
SWHeatFlux(I,J) = (SEAICE_waterAlbedo - 1. _d 0) |
|
./pkg/seaice/seaice_budget_ocean.F |
121 |
SEAICE_BUDGET_OCEAN |
& *swdown(I,J,bi,bj) |
|
./pkg/seaice/seaice_budget_ocean.F |
125 |
SEAICE_BUDGET_OCEAN |
lwflux(I,J,bi,bj) = - lwdownLoc(I,J)*SEAICE_emissivity |
|
./pkg/seaice/seaice_budget_ocean.F |
126 |
SEAICE_BUDGET_OCEAN |
& + D3*tsurfLoc(I,J)**4 |
|
./pkg/seaice/seaice_budget_ocean.F |
127 |
SEAICE_BUDGET_OCEAN |
sstdegC = tsurfLoc(I,J) - TMELT |
|
./pkg/seaice/seaice_budget_ocean.F |
129 |
SEAICE_BUDGET_OCEAN |
hs(I,J,bi,bj) = D1*UG(I,J)*(atempLoc(I,J)-tsurfLoc(I,J)) |
|
./pkg/seaice/seaice_budget_ocean.F |
134 |
SEAICE_BUDGET_OCEAN |
hl(I,J,bi,bj) = D1W*UG(I,J)*(AQH(I,J,bi,bj)-ssq) |
|
./pkg/seaice/seaice_budget_ocean.F |
136 |
SEAICE_BUDGET_OCEAN |
netHeatFlux(I,J)=SWHeatFlux(I,J) |
|
./pkg/seaice/seaice_budget_ocean.F |
137 |
SEAICE_BUDGET_OCEAN |
& + lwflux(I,J,bi,bj) |
|
./pkg/seaice/seaice_budget_ocean.F |
138 |
SEAICE_BUDGET_OCEAN |
& - hs(I,J,bi,bj) - hl(I,J,bi,bj) |
|
./pkg/seaice/seaice_budget_ocean.F |
141 |
SEAICE_BUDGET_OCEAN |
evap(i,j,bi,bj) = -hl(I,J,bi,bj) |
|
./pkg/seaice/seaice_budget_ocean.F |
85 |
SEAICE_BUDGET_OCEAN |
DO J=1,sNy |
|
./pkg/seaice/seaice_budget_ocean.F |
87 |
SEAICE_BUDGET_OCEAN |
netHeatFlux(I,J) = 0. _d 0 |
|
./pkg/seaice/seaice_budget_ocean.F |
88 |
SEAICE_BUDGET_OCEAN |
SWHeatFlux (I,J) = 0. _d 0 |
|
./pkg/seaice/seaice_budget_ocean.F |
92 |
SEAICE_BUDGET_OCEAN |
tsurfLoc (I,J) = TSURF(I,J) |
|
./pkg/seaice/seaice_budget_ocean.F |
95 |
SEAICE_BUDGET_OCEAN |
atempLoc (I,J) = MAX(celsius2K+MIN_ATEMP,ATEMP(I,J,bi,bj)) |
|
./pkg/seaice/seaice_budget_ocean.F |
98 |
SEAICE_BUDGET_OCEAN |
lwdownLoc(I,J) = MAX(MIN_LWDOWN,LWDOWN(I,J,bi,bj)) |
|
./pkg/seaice/seaice_calc_ice_strength.F |
117 |
SEAICE_CALC_ICE_STRENGTH |
PRESS0(I,J,bi,bj)=SEAICE_strength*tmpscal2 |
|
./pkg/seaice/seaice_calc_ice_strength.F |
119 |
SEAICE_CALC_ICE_STRENGTH |
ZMAX(I,J,bi,bj) = SEAICE_zetaMaxFac*PRESS0(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_ice_strength.F |
120 |
SEAICE_CALC_ICE_STRENGTH |
ZMIN(I,J,bi,bj) = SEAICE_zetaMin |
|
./pkg/seaice/seaice_calc_ice_strength.F |
121 |
SEAICE_CALC_ICE_STRENGTH |
PRESS0(I,J,bi,bj)= PRESS0(I,J,bi,bj)*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_ice_strength.F |
174 |
SEAICE_CALC_ICE_STRENGTH |
ZMAX(I,J,bi,bj) = SEAICE_zetaMaxFac*PRESS0(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_ice_strength.F |
175 |
SEAICE_CALC_ICE_STRENGTH |
ZMIN(I,J,bi,bj) = SEAICE_zetaMin |
|
./pkg/seaice/seaice_calc_ice_strength.F |
176 |
SEAICE_CALC_ICE_STRENGTH |
PRESS0(I,J,bi,bj)= PRESS0(I,J,bi,bj)*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
109 |
SEAICE_CALC_LHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_lhs.F |
111 |
SEAICE_CALC_LHS |
areaW(I,J) = 1. _d 0 |
|
./pkg/seaice/seaice_calc_lhs.F |
112 |
SEAICE_CALC_LHS |
areaS(I,J) = 1. _d 0 |
|
./pkg/seaice/seaice_calc_lhs.F |
121 |
SEAICE_CALC_LHS |
dragSym(I,J) = DWATN(I,J,bi,bj)*COSWAT |
|
./pkg/seaice/seaice_calc_lhs.F |
123 |
SEAICE_CALC_LHS |
& +CbotC(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
137 |
SEAICE_CALC_LHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_lhs.F |
139 |
SEAICE_CALC_LHS |
areaW(I,J) = 0.5 _d 0*(AREA(I,J,bi,bj)+AREA(I-1,J,bi,bj)) |
|
./pkg/seaice/seaice_calc_lhs.F |
140 |
SEAICE_CALC_LHS |
areaS(I,J) = 0.5 _d 0*(AREA(I,J,bi,bj)+AREA(I,J-1,bi,bj)) |
|
./pkg/seaice/seaice_calc_lhs.F |
144 |
SEAICE_CALC_LHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_lhs.F |
147 |
SEAICE_CALC_LHS |
uIceLHS(I,J,bi,bj) = |
|
./pkg/seaice/seaice_calc_lhs.F |
148 |
SEAICE_CALC_LHS |
& bdfAlpha*seaiceMassU(I,J,bi,bj)*recip_deltaT |
|
./pkg/seaice/seaice_calc_lhs.F |
149 |
SEAICE_CALC_LHS |
& *uIceLoc(I,J,bi,bj) - stressDivergenceX(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
151 |
SEAICE_CALC_LHS |
vIceLHS(I,J,bi,bj) = |
|
./pkg/seaice/seaice_calc_lhs.F |
152 |
SEAICE_CALC_LHS |
& bdfAlpha*seaiceMassV(I,J,bi,bj)*recip_deltaT |
|
./pkg/seaice/seaice_calc_lhs.F |
153 |
SEAICE_CALC_LHS |
& *vIceLoc(I,J,bi,bj) - stressDivergenceY(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
155 |
SEAICE_CALC_LHS |
uIceLHS(I,J,bi,bj) = uIceLHS(I,J,bi,bj) - 0.5 _d 0*( |
|
./pkg/seaice/seaice_calc_lhs.F |
156 |
SEAICE_CALC_LHS |
& seaiceMassC(I ,J,bi,bj) * _fCori(I ,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
157 |
SEAICE_CALC_LHS |
& * 0.5 _d 0*( vIceLoc(I ,J,bi,bj)+vIceLoc(I ,J+1,bi,bj) ) |
|
./pkg/seaice/seaice_calc_lhs.F |
158 |
SEAICE_CALC_LHS |
& + seaiceMassC(I-1,J,bi,bj) * _fCori(I-1,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
159 |
SEAICE_CALC_LHS |
& * 0.5 _d 0*( vIceLoc(I-1,J,bi,bj)+vIceLoc(I-1,J+1,bi,bj) ) |
|
./pkg/seaice/seaice_calc_lhs.F |
162 |
SEAICE_CALC_LHS |
vIceLHS(I,J,bi,bj) = vIceLHS(I,J,bi,bj) + 0.5 _d 0*( |
|
./pkg/seaice/seaice_calc_lhs.F |
163 |
SEAICE_CALC_LHS |
& seaiceMassC(I,J ,bi,bj) * _fCori(I,J ,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
164 |
SEAICE_CALC_LHS |
& * 0.5 _d 0*( uIceLoc(I,J ,bi,bj)+uIceLoc(I+1, J,bi,bj) ) |
|
./pkg/seaice/seaice_calc_lhs.F |
165 |
SEAICE_CALC_LHS |
& + seaiceMassC(I,J-1,bi,bj) * _fCori(I,J-1,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
166 |
SEAICE_CALC_LHS |
& * 0.5 _d 0*( uIceLoc(I,J-1,bi,bj)+uIceLoc(I+1,J-1,bi,bj) ) |
|
./pkg/seaice/seaice_calc_lhs.F |
169 |
SEAICE_CALC_LHS |
uIceLHS(I,J,bi,bj) = uIceLHS(I,J,bi,bj) + ( |
|
./pkg/seaice/seaice_calc_lhs.F |
170 |
SEAICE_CALC_LHS |
& 0.5 _d 0 * ( dragSym(I,J)+dragSym(I-1,J) ) |
|
./pkg/seaice/seaice_calc_lhs.F |
171 |
SEAICE_CALC_LHS |
& * uIceLoc(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
172 |
SEAICE_CALC_LHS |
& - SIGN(SINWAT, _fCori(I,J,bi,bj))* 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_lhs.F |
173 |
SEAICE_CALC_LHS |
& ( DWATN(I ,J,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_lhs.F |
174 |
SEAICE_CALC_LHS |
& (vIceLoc(I ,J,bi,bj)+vIceLoc(I ,J+1,bi,bj)) |
|
./pkg/seaice/seaice_calc_lhs.F |
175 |
SEAICE_CALC_LHS |
& + DWATN(I-1,J,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_lhs.F |
176 |
SEAICE_CALC_LHS |
& (vIceLoc(I-1,J,bi,bj)+vIceLoc(I-1,J+1,bi,bj)) |
|
./pkg/seaice/seaice_calc_lhs.F |
177 |
SEAICE_CALC_LHS |
& ) ) * areaW(I,J) |
|
./pkg/seaice/seaice_calc_lhs.F |
179 |
SEAICE_CALC_LHS |
vIceLHS(I,J,bi,bj) = vIceLHS(I,J,bi,bj) + ( |
|
./pkg/seaice/seaice_calc_lhs.F |
180 |
SEAICE_CALC_LHS |
& 0.5 _d 0 * ( dragSym(I,J)+dragSym(I,J-1) ) |
|
./pkg/seaice/seaice_calc_lhs.F |
181 |
SEAICE_CALC_LHS |
& * vIceLoc(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
182 |
SEAICE_CALC_LHS |
& + SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_lhs.F |
183 |
SEAICE_CALC_LHS |
& ( DWATN(I,J ,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_lhs.F |
184 |
SEAICE_CALC_LHS |
& (uIceLoc(I,J ,bi,bj)+uIceLoc(I+1,J ,bi,bj)) |
|
./pkg/seaice/seaice_calc_lhs.F |
185 |
SEAICE_CALC_LHS |
& + DWATN(I,J-1,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_lhs.F |
186 |
SEAICE_CALC_LHS |
& (uIceLoc(I,J-1,bi,bj)+uIceLoc(I+1,J-1,bi,bj)) |
|
./pkg/seaice/seaice_calc_lhs.F |
187 |
SEAICE_CALC_LHS |
& ) ) * areaS(I,J) |
|
./pkg/seaice/seaice_calc_lhs.F |
189 |
SEAICE_CALC_LHS |
uIceLHS(I,J,bi,bj) = uIceLHS(I,J,bi,bj)*maskinW(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
190 |
SEAICE_CALC_LHS |
vIceLHS(I,J,bi,bj) = vIceLHS(I,J,bi,bj)*maskinS(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_lhs.F |
195 |
SEAICE_CALC_LHS |
DO J=1-Oly,sNy+Oly |
|
./pkg/seaice/seaice_calc_lhs.F |
197 |
SEAICE_CALC_LHS |
gUmom(I,J) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_lhs.F |
198 |
SEAICE_CALC_LHS |
gVmom(I,J) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_lhs.F |
208 |
SEAICE_CALC_LHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_lhs.F |
210 |
SEAICE_CALC_LHS |
uIceLHS(I,J,bi,bj) = uIceLHS(I,J,bi,bj) - gUmom(I,J) |
|
./pkg/seaice/seaice_calc_lhs.F |
211 |
SEAICE_CALC_LHS |
vIceLHS(I,J,bi,bj) = vIceLHS(I,J,bi,bj) - gVmom(I,J) |
|
./pkg/seaice/seaice_calc_residual.F |
137 |
SEAICE_CALC_RESIDUAL |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_residual.F |
139 |
SEAICE_CALC_RESIDUAL |
uIceRes(I,J,bi,bj) = uIceLHS(I,J,bi,bj) - uIceRHS(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_residual.F |
140 |
SEAICE_CALC_RESIDUAL |
vIceRes(I,J,bi,bj) = vIceLHS(I,J,bi,bj) - vIceRHS(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_residual.F |
70 |
SEAICE_CALC_RESIDUAL |
DO J=1-Oly,sNy+Oly |
|
./pkg/seaice/seaice_calc_residual.F |
72 |
SEAICE_CALC_RESIDUAL |
uIceLHS(I,J,bi,bj) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_residual.F |
73 |
SEAICE_CALC_RESIDUAL |
vIceLHS(I,J,bi,bj) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_residual.F |
74 |
SEAICE_CALC_RESIDUAL |
uIceRHS(I,J,bi,bj) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_residual.F |
75 |
SEAICE_CALC_RESIDUAL |
vIceRHS(I,J,bi,bj) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_rhs.F |
100 |
SEAICE_CALC_RHS |
& ) )*areaW(I,J) |
|
./pkg/seaice/seaice_calc_rhs.F |
101 |
SEAICE_CALC_RHS |
vIceRHS(I,J,bi,bj) = FORCEY(I,J,bi,bj) + ( |
|
./pkg/seaice/seaice_calc_rhs.F |
102 |
SEAICE_CALC_RHS |
& 0.5 _d 0 * ( DWATN(I,J,bi,bj)+DWATN(I,J-1,bi,bj) ) * |
|
./pkg/seaice/seaice_calc_rhs.F |
103 |
SEAICE_CALC_RHS |
& COSWAT * vVel(I,J,kSrf,bi,bj) |
|
./pkg/seaice/seaice_calc_rhs.F |
104 |
SEAICE_CALC_RHS |
& + SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_rhs.F |
105 |
SEAICE_CALC_RHS |
& ( DWATN(I,J ,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_rhs.F |
106 |
SEAICE_CALC_RHS |
& ( uVel(I,J ,kSrf,bi,bj)+uVel(I+1,J ,kSrf,bi,bj)) |
|
./pkg/seaice/seaice_calc_rhs.F |
107 |
SEAICE_CALC_RHS |
& + DWATN(I,J-1,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_rhs.F |
108 |
SEAICE_CALC_RHS |
& ( uVel(I,J-1,kSrf,bi,bj)+uVel(I+1,J-1,kSrf,bi,bj)) |
|
./pkg/seaice/seaice_calc_rhs.F |
109 |
SEAICE_CALC_RHS |
& ) )*areaS(I,J) |
|
./pkg/seaice/seaice_calc_rhs.F |
111 |
SEAICE_CALC_RHS |
uIceRHS(I,J,bi,bj) = uIceRHS(I,J,bi,bj)*maskinW(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_rhs.F |
112 |
SEAICE_CALC_RHS |
vIceRHS(I,J,bi,bj) = vIceRHS(I,J,bi,bj)*maskinS(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_rhs.F |
70 |
SEAICE_CALC_RHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_rhs.F |
72 |
SEAICE_CALC_RHS |
areaW(I,J) = 1. _d 0 |
|
./pkg/seaice/seaice_calc_rhs.F |
73 |
SEAICE_CALC_RHS |
areaS(I,J) = 1. _d 0 |
|
./pkg/seaice/seaice_calc_rhs.F |
80 |
SEAICE_CALC_RHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_rhs.F |
82 |
SEAICE_CALC_RHS |
areaW(I,J) = 0.5 _d 0*(AREA(I,J,bi,bj)+AREA(I-1,J,bi,bj)) |
|
./pkg/seaice/seaice_calc_rhs.F |
83 |
SEAICE_CALC_RHS |
areaS(I,J) = 0.5 _d 0*(AREA(I,J,bi,bj)+AREA(I,J-1,bi,bj)) |
|
./pkg/seaice/seaice_calc_rhs.F |
87 |
SEAICE_CALC_RHS |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_rhs.F |
92 |
SEAICE_CALC_RHS |
uIceRHS(I,J,bi,bj) = FORCEX(I,J,bi,bj) + ( |
|
./pkg/seaice/seaice_calc_rhs.F |
93 |
SEAICE_CALC_RHS |
& 0.5 _d 0 * ( DWATN(I,J,bi,bj)+DWATN(I-1,J,bi,bj) ) * |
|
./pkg/seaice/seaice_calc_rhs.F |
94 |
SEAICE_CALC_RHS |
& COSWAT * uVel(I,J,kSrf,bi,bj) |
|
./pkg/seaice/seaice_calc_rhs.F |
95 |
SEAICE_CALC_RHS |
& - SIGN(SINWAT, _fCori(I,J,bi,bj))* 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_rhs.F |
96 |
SEAICE_CALC_RHS |
& ( DWATN(I ,J,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_rhs.F |
97 |
SEAICE_CALC_RHS |
& ( vVel(I ,J,kSrf,bi,bj)+vVel(I ,J+1,kSrf,bi,bj) ) |
|
./pkg/seaice/seaice_calc_rhs.F |
98 |
SEAICE_CALC_RHS |
& + DWATN(I-1,J,bi,bj) * 0.5 _d 0 * |
|
./pkg/seaice/seaice_calc_rhs.F |
99 |
SEAICE_CALC_RHS |
& ( vVel(I-1,J,kSrf,bi,bj)+vVel(I-1,J+1,kSrf,bi,bj) ) |
|
./pkg/seaice/seaice_calc_stressdiv.F |
82 |
SEAICE_CALC_STRESSDIV |
DO J=1,sNy |
|
./pkg/seaice/seaice_calc_stressdiv.F |
86 |
SEAICE_CALC_STRESSDIV |
& - sig11(i-1,J ) * _dyF(i-1,J,bi,bj) |
|
./pkg/seaice/seaice_calc_stressdiv.F |
93 |
SEAICE_CALC_STRESSDIV |
& + sig12(i+1,J ) * _dyU(i+1,J,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
106 |
SEAICE_CALC_VISCOSITIES |
& ( e12(I,J ,bi,bj) + e12(I+1,J ,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
107 |
SEAICE_CALC_VISCOSITIES |
& + e12(I,J+1,bi,bj) + e12(I+1,J+1,bi,bj) ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
116 |
SEAICE_CALC_VISCOSITIES |
e12Csq(i,j) = 0.25 _d 0 * recip_rA(I,J,bi,bj) * |
|
./pkg/seaice/seaice_calc_viscosities.F |
117 |
SEAICE_CALC_VISCOSITIES |
& ( rAz(I ,J ,bi,bj)*e12(I ,J ,bi,bj)**2 |
|
./pkg/seaice/seaice_calc_viscosities.F |
118 |
SEAICE_CALC_VISCOSITIES |
& + rAz(I+1,J ,bi,bj)*e12(I+1,J ,bi,bj)**2 |
|
./pkg/seaice/seaice_calc_viscosities.F |
119 |
SEAICE_CALC_VISCOSITIES |
& + rAz(I ,J+1,bi,bj)*e12(I ,J+1,bi,bj)**2 |
|
./pkg/seaice/seaice_calc_viscosities.F |
120 |
SEAICE_CALC_VISCOSITIES |
& + rAz(I+1,J+1,bi,bj)*e12(I+1,J+1,bi,bj)**2 ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
137 |
SEAICE_CALC_VISCOSITIES |
deltaC(I,J,bi,bj) = 0. _d 0 |
|
./pkg/seaice/seaice_calc_viscosities.F |
139 |
SEAICE_CALC_VISCOSITIES |
& deltaC(I,J,bi,bj) = SQRT(deltaCsq) |
|
./pkg/seaice/seaice_calc_viscosities.F |
141 |
SEAICE_CALC_VISCOSITIES |
deltaC(I,J,bi,bj) = SQRT(deltaCsq) |
|
./pkg/seaice/seaice_calc_viscosities.F |
149 |
SEAICE_CALC_VISCOSITIES |
deltaCreg = MAX(deltaC(I,J,bi,bj),SEAICE_deltaMin) |
|
./pkg/seaice/seaice_calc_viscosities.F |
156 |
SEAICE_CALC_VISCOSITIES |
zeta (I,J,bi,bj) = ZMAX(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
159 |
SEAICE_CALC_VISCOSITIES |
zeta (I,J,bi,bj) = HALF*( press0(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
160 |
SEAICE_CALC_VISCOSITIES |
& * ( 1. _d 0 + tnsFac(I,J,bi,bj) ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
163 |
SEAICE_CALC_VISCOSITIES |
zeta (I,J,bi,bj) = MIN(ZMAX(I,J,bi,bj),zeta(I,J,bi,bj)) |
|
./pkg/seaice/seaice_calc_viscosities.F |
165 |
SEAICE_CALC_VISCOSITIES |
zeta (I,J,bi,bj) = MAX(ZMIN(I,J,bi,bj),zeta(I,J,bi,bj)) |
|
./pkg/seaice/seaice_calc_viscosities.F |
167 |
SEAICE_CALC_VISCOSITIES |
zeta (I,J,bi,bj) = zeta(I,J,bi,bj)*HEFFM(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
168 |
SEAICE_CALC_VISCOSITIES |
eta (I,J,bi,bj) = recip_e2*zeta(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
170 |
SEAICE_CALC_VISCOSITIES |
press(I,J,bi,bj) = |
|
./pkg/seaice/seaice_calc_viscosities.F |
171 |
SEAICE_CALC_VISCOSITIES |
& ( press0(I,J,bi,bj)*( 1. _d 0 - SEAICEpressReplFac ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
172 |
SEAICE_CALC_VISCOSITIES |
& + TWO*zeta(I,J,bi,bj)*deltaC(I,J,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
173 |
SEAICE_CALC_VISCOSITIES |
& * SEAICEpressReplFac/( 1. _d 0 + tnsFac(I,J,bi,bj) ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
174 |
SEAICE_CALC_VISCOSITIES |
& ) * ( 1. _d 0 - tnsFac(I,J,bi,bj) ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
189 |
SEAICE_CALC_VISCOSITIES |
etaMax = ( 0.5 _d 0*press(I,J,bi,bj)-zeta(I,J,bi,bj)*ep |
|
./pkg/seaice/seaice_calc_viscosities.F |
191 |
SEAICE_CALC_VISCOSITIES |
eta(I,J,bi,bj) = MIN(eta(I,J,bi,bj),etaMax) |
|
./pkg/seaice/seaice_calc_viscosities.F |
199 |
SEAICE_CALC_VISCOSITIES |
sumNorm = maskC(I,J, k,bi,bj)+maskC(I-1,J, k,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
200 |
SEAICE_CALC_VISCOSITIES |
& + maskC(I,J-1,k,bi,bj)+maskC(I-1,J-1,k,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
202 |
SEAICE_CALC_VISCOSITIES |
etaZ(I,J,bi,bj) = sumNorm * |
|
./pkg/seaice/seaice_calc_viscosities.F |
203 |
SEAICE_CALC_VISCOSITIES |
& ( eta (I,J ,bi,bj) + eta (I-1,J ,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
204 |
SEAICE_CALC_VISCOSITIES |
& + eta (I,J-1,bi,bj) + eta (I-1,J-1,bi,bj) ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
205 |
SEAICE_CALC_VISCOSITIES |
zetaZ(I,J,bi,bj) = sumNorm * |
|
./pkg/seaice/seaice_calc_viscosities.F |
206 |
SEAICE_CALC_VISCOSITIES |
& ( zeta(I,J ,bi,bj) + zeta(I-1,J ,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
207 |
SEAICE_CALC_VISCOSITIES |
& + zeta(I,J-1,bi,bj) + zeta(I-1,J-1,bi,bj) ) |
|
./pkg/seaice/seaice_calc_viscosities.F |
214 |
SEAICE_CALC_VISCOSITIES |
DO J=1-OLy+1,sNy+OLy-1 |
|
./pkg/seaice/seaice_calc_viscosities.F |
216 |
SEAICE_CALC_VISCOSITIES |
maskZ = maskC(I,J, k,bi,bj)*maskC(I-1,J, k,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
217 |
SEAICE_CALC_VISCOSITIES |
& * maskC(I,J-1,k,bi,bj)*maskC(I-1,J-1,k,bi,bj) |
|
./pkg/seaice/seaice_calc_viscosities.F |
218 |
SEAICE_CALC_VISCOSITIES |
etaZ (I,J,bi,bj) = etaZ(I,J,bi,bj) * maskZ |
|
./pkg/seaice/seaice_calc_viscosities.F |
219 |
SEAICE_CALC_VISCOSITIES |
zetaZ(I,J,bi,bj) = zetaZ(I,J,bi,bj) * maskZ |
|
./pkg/seaice/seaice_diagnostics_state.F |
139 |
SEAICE_DIAGNOSTICS_STATE |
IF ( tensileStrFac(I,J,bi,bj) .EQ. 1. _d 0) THEN |
|
./pkg/seaice/seaice_diagnostics_state.F |
144 |
SEAICE_DIAGNOSTICS_STATE |
tIce(I,J,bi,bj) = press0(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
148 |
SEAICE_DIAGNOSTICS_STATE |
tIce(I,J,bi,bj) = tensileStrFac(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
149 |
SEAICE_DIAGNOSTICS_STATE |
& *press(I,J,bi,bj)/(1. _d 0 - tensileStrFac(I,J,bi,bj)) |
|
./pkg/seaice/seaice_diagnostics_state.F |
169 |
SEAICE_DIAGNOSTICS_STATE |
sigp = seaice_sigma1(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
170 |
SEAICE_DIAGNOSTICS_STATE |
sigm = seaice_sigma2(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
171 |
SEAICE_DIAGNOSTICS_STATE |
sig12(I,J) = 0.25 _d 0 * |
|
./pkg/seaice/seaice_diagnostics_state.F |
172 |
SEAICE_DIAGNOSTICS_STATE |
& ( seaice_sigma12(I, J, bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
173 |
SEAICE_DIAGNOSTICS_STATE |
& + seaice_sigma12(I+1,J, bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
174 |
SEAICE_DIAGNOSTICS_STATE |
& + seaice_sigma12(I+1,J+1,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
175 |
SEAICE_DIAGNOSTICS_STATE |
& + seaice_sigma12(I ,J+1,bi,bj) ) |
|
./pkg/seaice/seaice_diagnostics_state.F |
176 |
SEAICE_DIAGNOSTICS_STATE |
sigTmp = SQRT( sigm*sigm + 4. _d 0*sig12(I,J)*sig12(I,J) ) |
|
./pkg/seaice/seaice_diagnostics_state.F |
178 |
SEAICE_DIAGNOSTICS_STATE |
IF ( press0(I,J,bi,bj) .GT. 1. _d -13 ) |
|
./pkg/seaice/seaice_diagnostics_state.F |
179 |
SEAICE_DIAGNOSTICS_STATE |
& recip_prs = 1./press0(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
180 |
SEAICE_DIAGNOSTICS_STATE |
sig1(I,J) = 0.5*(sigp + sigTmp)*recip_prs |
|
./pkg/seaice/seaice_diagnostics_state.F |
181 |
SEAICE_DIAGNOSTICS_STATE |
sig2(I,J) = 0.5*(sigp - sigTmp)*recip_prs |
|
./pkg/seaice/seaice_diagnostics_state.F |
211 |
SEAICE_DIAGNOSTICS_STATE |
sigp = sig11(I,J) + sig22(I,J) |
|
./pkg/seaice/seaice_diagnostics_state.F |
212 |
SEAICE_DIAGNOSTICS_STATE |
sigm = sig11(I,J) - sig22(I,J) |
|
./pkg/seaice/seaice_diagnostics_state.F |
221 |
SEAICE_DIAGNOSTICS_STATE |
sigTmp = 2.*eta(I,J,bi,bj) * 0.25 _d 0 * |
|
./pkg/seaice/seaice_diagnostics_state.F |
222 |
SEAICE_DIAGNOSTICS_STATE |
& (e12(I,J,bi,bj) + e12(I+1,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
223 |
SEAICE_DIAGNOSTICS_STATE |
& +e12(I,J+1,bi,bj)+e12(I+1,J+1,bi,bj)) |
|
./pkg/seaice/seaice_diagnostics_state.F |
226 |
SEAICE_DIAGNOSTICS_STATE |
IF ( press0(I,J,bi,bj) .GT. 1. _d -13 ) |
|
./pkg/seaice/seaice_diagnostics_state.F |
227 |
SEAICE_DIAGNOSTICS_STATE |
& recip_prs = 1./press0(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
228 |
SEAICE_DIAGNOSTICS_STATE |
sig1(I,J) = 0.5*(sigp + sigTmp)*recip_prs |
|
./pkg/seaice/seaice_diagnostics_state.F |
229 |
SEAICE_DIAGNOSTICS_STATE |
sig2(I,J) = 0.5*(sigp - sigTmp)*recip_prs |
|
./pkg/seaice/seaice_diagnostics_state.F |
242 |
SEAICE_DIAGNOSTICS_STATE |
sigm = e11(I,J,bi,bj) - e22(I,J,bi,bj) |
|
./pkg/seaice/seaice_diagnostics_state.F |
244 |
SEAICE_DIAGNOSTICS_STATE |
& ( e12(I, J, bi,bj)**2 + e12(I+1,J, bi,bj)**2 |
|
./pkg/seaice/seaice_diagnostics_state.F |
245 |
SEAICE_DIAGNOSTICS_STATE |
& + e12(I+1,J+1,bi,bj)**2 + e12(I ,J+1,bi,bj)**2 ) |
|
./pkg/seaice/seaice_diagnostics_state.F |
248 |
SEAICE_DIAGNOSTICS_STATE |
sig1(I,J) = sqrt(sigm*sigm + sigTmp) |