- Timestamp:
- 2018-07-29T11:23:51+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/DIA/dia25h.F90
r9939 r10009 94 94 ! ------------------------- ! 95 95 cnt_25h = 1 ! sets the first value of sum at timestep 1 (note - should strictly be at timestep zero so before values used where possible) 96 tn_25h (:,:,:) = tsb 97 sn_25h (:,:,:) = tsb 98 sshn_25h(:,:) = ssh b(:,:)99 un_25h (:,:,:) = ub 100 vn_25h (:,:,:) = vb 101 wn_25h (:,:,:) = wn 102 avt_25h (:,:,:) = avt 103 avm_25h (:,:,:) = avm 96 tn_25h (:,:,:) = tsb(:,:,:,jp_tem) 97 sn_25h (:,:,:) = tsb(:,:,:,jp_sal) 98 sshn_25h(:,:) = ssh(:,:,Nbb) 99 un_25h (:,:,:) = ub (:,:,:) 100 vn_25h (:,:,:) = vb (:,:,:) 101 wn_25h (:,:,:) = wn (:,:,:) 102 avt_25h (:,:,:) = avt(:,:,:) 103 avm_25h (:,:,:) = avm(:,:,:) 104 104 IF( ln_zdftke ) THEN 105 105 en_25h(:,:,:) = en(:,:,:) … … 156 156 ENDIF 157 157 158 tn_25h (:,:,:) = tn_25h (:,:,:) + tsn 159 sn_25h (:,:,:) = sn_25h (:,:,:) + tsn 160 sshn_25h(:,:) = sshn_25h(:,:) + ssh n(:,:)161 un_25h (:,:,:) = un_25h (:,:,:) + un 162 vn_25h (:,:,:) = vn_25h (:,:,:) + vn 163 wn_25h (:,:,:) = wn_25h (:,:,:) + wn 164 avt_25h (:,:,:) = avt_25h (:,:,:) + avt 165 avm_25h (:,:,:) = avm_25h (:,:,:) + avm 158 tn_25h (:,:,:) = tn_25h (:,:,:) + tsn(:,:,:,jp_tem) 159 sn_25h (:,:,:) = sn_25h (:,:,:) + tsn(:,:,:,jp_sal) 160 sshn_25h(:,:) = sshn_25h(:,:) + ssh(:,:,Nnn) 161 un_25h (:,:,:) = un_25h (:,:,:) + un (:,:,:) 162 vn_25h (:,:,:) = vn_25h (:,:,:) + vn (:,:,:) 163 wn_25h (:,:,:) = wn_25h (:,:,:) + wn (:,:,:) 164 avt_25h (:,:,:) = avt_25h (:,:,:) + avt(:,:,:) 165 avm_25h (:,:,:) = avm_25h (:,:,:) + avm(:,:,:) 166 166 IF( ln_zdftke ) THEN 167 167 en_25h(:,:,:) = en_25h (:,:,:) + en(:,:,:) … … 206 206 zmdi=1.e+20 !missing data indicator for masking 207 207 ! write tracers (instantaneous) 208 zw3d(:,:,:) = tn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 209 CALL iom_put("temper25h", zw3d) ! potential temperature 210 zw3d(:,:,:) = sn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 211 CALL iom_put( "salin25h", zw3d ) ! salinity 212 zw2d(:,:) = sshn_25h(:,:)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1)) 213 CALL iom_put( "ssh25h", zw2d ) ! sea surface 208 zw3d(:,:,:) = tn_25h (:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put("temper25h" , zw3d ) ! potential temperature 209 zw3d(:,:,:) = sn_25h (:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put( "salin25h" , zw3d ) ! salinity 210 zw2d(:,:) = sshn_25h(:,:) *ssmask(:,:) + zmdi*(1.0-tmask(:,:,1)) ; CALL iom_put( "ssh25h" , zw2d ) ! sea surface 214 211 ! Write velocities (instantaneous) 215 zw3d(:,:,:) = un_25h(:,:,:)*umask(:,:,:) + zmdi*(1.0-umask(:,:,:)) 216 CALL iom_put("vozocrtx25h", zw3d) ! i-current 217 zw3d(:,:,:) = vn_25h(:,:,:)*vmask(:,:,:) + zmdi*(1.0-vmask(:,:,:)) 218 CALL iom_put("vomecrty25h", zw3d ) ! j-current 219 zw3d(:,:,:) = wn_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 220 CALL iom_put("vomecrtz25h", zw3d ) ! k-current 212 zw3d(:,:,:) = un_25h (:,:,:)* umask(:,:,:) + zmdi*(1.0-umask(:,:,:)) ; CALL iom_put("vozocrtx25h", zw3d ) ! i-current 213 zw3d(:,:,:) = vn_25h (:,:,:)* vmask(:,:,:) + zmdi*(1.0-vmask(:,:,:)) ; CALL iom_put("vomecrty25h", zw3d ) ! j-current 214 zw3d(:,:,:) = wn_25h (:,:,:)* wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put("vomecrtz25h", zw3d ) ! k-current 221 215 ! Write vertical physics 222 zw3d(:,:,:) = avt_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 223 CALL iom_put("avt25h", zw3d ) ! diffusivity 224 zw3d(:,:,:) = avm_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 225 CALL iom_put("avm25h", zw3d) ! viscosity 226 IF( ln_zdftke ) THEN 227 zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 228 CALL iom_put("tke25h", zw3d) ! tke 229 ENDIF 230 IF( ln_zdfgls ) THEN 231 zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 232 CALL iom_put("tke25h", zw3d) ! tke 233 zw3d(:,:,:) = rmxln_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 234 CALL iom_put( "mxln25h",zw3d) 216 zw3d(:,:,:) = avt_25h (:,:,:)* wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put("avt25h" , zw3d ) ! diffusivity 217 zw3d(:,:,:) = avm_25h (:,:,:)* wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put("avm25h" , zw3d ) ! viscosity 218 IF( ln_zdftke ) THEN 219 zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put("tke25h" , zw3d ) ! tke 220 ENDIF 221 IF( ln_zdfgls ) THEN 222 zw3d(:,:,:) = en_25h (:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put("tke25h" , zw3d ) ! tke 223 zw3d(:,:,:) = rmxln_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) ; CALL iom_put( "mxln25h", zw3d ) 235 224 ENDIF 236 225 ! 237 226 ! After the write reset the values to cnt=1 and sum values equal current value 238 tn_25h (:,:,:) = tsn 239 sn_25h (:,:,:) = tsn 240 sshn_25h(:,:) = ssh n(:,:)241 un_25h (:,:,:) = un 242 vn_25h (:,:,:) = vn 243 wn_25h (:,:,:) = wn 244 avt_25h (:,:,:) = avt 245 avm_25h (:,:,:) = avm 227 tn_25h (:,:,:) = tsn(:,:,:,jp_tem) 228 sn_25h (:,:,:) = tsn(:,:,:,jp_sal) 229 sshn_25h(:,:) = ssh(:,:,Nnn) 230 un_25h (:,:,:) = un (:,:,:) 231 vn_25h (:,:,:) = vn (:,:,:) 232 wn_25h (:,:,:) = wn (:,:,:) 233 avt_25h (:,:,:) = avt(:,:,:) 234 avm_25h (:,:,:) = avm(:,:,:) 246 235 IF( ln_zdftke ) THEN 247 236 en_25h(:,:,:) = en(:,:,:)
Note: See TracChangeset
for help on using the changeset viewer.