Changeset 10955 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfosm.F90
- Timestamp:
- 2019-05-10T11:15:50+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfosm.F90
r10954 r10955 941 941 ! Stable conditions 942 942 DO jk = 2, ibld(ji,jj) 943 zznd_d=gdepw _n(ji,jj,jk) / dstokes(ji,jj)943 zznd_d=gdepw(ji,jj,jk,Kmm) / dstokes(ji,jj) 944 944 ghamt(ji,jj,jk) = ghamt(ji,jj,jk) + 1.5 * EXP ( -0.9 * zznd_d ) & 945 945 & * ( 1.0 - EXP ( -4.0 * zznd_d ) ) * zsc_wth_1(ji,jj) … … 1066 1066 IF ( lconv(ji,jj) ) THEN 1067 1067 DO jk = 2, imld(ji,jj) 1068 zznd_ml=gdepw _n(ji,jj,jk) / zhml(ji,jj)1068 zznd_ml=gdepw(ji,jj,jk,Kmm) / zhml(ji,jj) 1069 1069 ghamt(ji,jj,jk) = ghamt(ji,jj,jk) + 0.3 * zsc_wth_1(ji,jj) & 1070 1070 & * ( -2.0 + 2.75 * ( ( 1.0 + 0.6 * zznd_ml**4 ) & … … 1188 1188 IF ( lconv(ji,jj) ) THEN 1189 1189 DO jk = 1, imld(ji,jj) - 1 1190 znd=gdepw _n(ji,jj,jk) / zhml(ji,jj)1190 znd=gdepw(ji,jj,jk,Kmm) / zhml(ji,jj) 1191 1191 ghamt(ji,jj,jk) = ghamt(ji,jj,jk) + zwth_ent(ji,jj) * znd 1192 1192 ghams(ji,jj,jk) = ghams(ji,jj,jk) + zws_ent(ji,jj) * znd
Note: See TracChangeset
for help on using the changeset viewer.