Changeset 12377 for NEMO/trunk/src/OCE/DIA/dia25h.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/OCE/DIA/dia25h.F90
r11536 r12377 39 39 CONTAINS 40 40 41 SUBROUTINE dia_25h_init 41 SUBROUTINE dia_25h_init( Kbb ) 42 42 !!--------------------------------------------------------------------------- 43 43 !! *** ROUTINE dia_25h_init *** … … 47 47 !! ** Method : Read namelist 48 48 !!--------------------------------------------------------------------------- 49 INTEGER, INTENT(in) :: Kbb ! Time level index 50 ! 49 51 INTEGER :: ios ! Local integer output status for namelist read 50 52 INTEGER :: ierror ! Local integer for memory allocation … … 53 55 !!---------------------------------------------------------------------- 54 56 ! 55 REWIND ( numnam_ref ) ! Read Namelist nam_dia25h in reference namelist : 25hour mean diagnostics56 57 READ ( numnam_ref, nam_dia25h, IOSTAT=ios, ERR= 901 ) 57 58 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_dia25h in reference namelist' ) 58 REWIND( numnam_cfg ) ! Namelist nam_dia25h in configuration namelist 25hour diagnostics59 59 READ ( numnam_cfg, nam_dia25h, IOSTAT = ios, ERR = 902 ) 60 60 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nam_dia25h in configuration namelist' ) … … 95 95 ! ------------------------- ! 96 96 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) 97 tn_25h (:,:,:) = ts b (:,:,:,jp_tem)98 sn_25h (:,:,:) = ts b (:,:,:,jp_sal)99 sshn_25h(:,:) = ssh b(:,:)100 un_25h (:,:,:) = u b (:,:,:)101 vn_25h (:,:,:) = v b (:,:,:)97 tn_25h (:,:,:) = ts (:,:,:,jp_tem,Kbb) 98 sn_25h (:,:,:) = ts (:,:,:,jp_sal,Kbb) 99 sshn_25h(:,:) = ssh(:,:,Kbb) 100 un_25h (:,:,:) = uu (:,:,:,Kbb) 101 vn_25h (:,:,:) = vv (:,:,:,Kbb) 102 102 avt_25h (:,:,:) = avt (:,:,:) 103 103 avm_25h (:,:,:) = avm (:,:,:) … … 116 116 117 117 118 SUBROUTINE dia_25h( kt )118 SUBROUTINE dia_25h( kt, Kmm ) 119 119 !!---------------------------------------------------------------------- 120 120 !! *** ROUTINE dia_25h *** … … 125 125 !!---------------------------------------------------------------------- 126 126 INTEGER, INTENT(in) :: kt ! ocean time-step index 127 INTEGER, INTENT(in) :: Kmm ! ocean time level index 127 128 !! 128 129 INTEGER :: ji, jj, jk … … 150 151 ! wn_25h could not be initialised in dia_25h_init, so we do it here instead 151 152 IF( kt == nn_it000 ) THEN 152 wn_25h(:,:,:) = w n(:,:,:)153 wn_25h(:,:,:) = ww(:,:,:) 153 154 ENDIF 154 155 … … 161 162 ENDIF 162 163 163 tn_25h (:,:,:) = tn_25h (:,:,:) + ts n (:,:,:,jp_tem)164 sn_25h (:,:,:) = sn_25h (:,:,:) + ts n (:,:,:,jp_sal)165 sshn_25h(:,:) = sshn_25h(:,:) + ssh n(:,:)166 un_25h (:,:,:) = un_25h (:,:,:) + u n (:,:,:)167 vn_25h (:,:,:) = vn_25h (:,:,:) + v n (:,:,:)168 wn_25h (:,:,:) = wn_25h (:,:,:) + w n(:,:,:)164 tn_25h (:,:,:) = tn_25h (:,:,:) + ts (:,:,:,jp_tem,Kmm) 165 sn_25h (:,:,:) = sn_25h (:,:,:) + ts (:,:,:,jp_sal,Kmm) 166 sshn_25h(:,:) = sshn_25h(:,:) + ssh(:,:,Kmm) 167 un_25h (:,:,:) = un_25h (:,:,:) + uu (:,:,:,Kmm) 168 vn_25h (:,:,:) = vn_25h (:,:,:) + vv (:,:,:,Kmm) 169 wn_25h (:,:,:) = wn_25h (:,:,:) + ww (:,:,:) 169 170 avt_25h (:,:,:) = avt_25h (:,:,:) + avt (:,:,:) 170 171 avm_25h (:,:,:) = avm_25h (:,:,:) + avm (:,:,:) … … 245 246 ! 246 247 ! After the write reset the values to cnt=1 and sum values equal current value 247 tn_25h (:,:,:) = ts n (:,:,:,jp_tem)248 sn_25h (:,:,:) = ts n (:,:,:,jp_sal)249 sshn_25h(:,:) = ssh n(:,:)250 un_25h (:,:,:) = u n (:,:,:)251 vn_25h (:,:,:) = v n (:,:,:)252 wn_25h (:,:,:) = w n(:,:,:)248 tn_25h (:,:,:) = ts (:,:,:,jp_tem,Kmm) 249 sn_25h (:,:,:) = ts (:,:,:,jp_sal,Kmm) 250 sshn_25h(:,:) = ssh(:,:,Kmm) 251 un_25h (:,:,:) = uu (:,:,:,Kmm) 252 vn_25h (:,:,:) = vv (:,:,:,Kmm) 253 wn_25h (:,:,:) = ww (:,:,:) 253 254 avt_25h (:,:,:) = avt (:,:,:) 254 255 avm_25h (:,:,:) = avm (:,:,:)
Note: See TracChangeset
for help on using the changeset viewer.