Changeset 6802
- Timestamp:
- 2016-07-13T13:04:12+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_25hr_mean_assim_bkg/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90
r6775 r6802 176 176 CALL iom_open( c_asmbkg, inum, ldwrt = .TRUE., kiolib = jprstlib) 177 177 ! 178 IF( nitbkg_r == nit000 - 1 ) THEN ! Treat special case when nitbkg = 0179 zdate = REAL( ndastp )180 #if defined key_zdftke181 ! lk_zdftke=T : Read turbulent kinetic energy ( en )182 IF(lwp) WRITE(numout,*) ' Reading TKE (en) from restart...'183 CALL tke_rst( nit000, 'READ' ) ! lk_zdftke=T : Read turbulent kinetic energy ( en )184 185 #endif186 ELSE187 zdate = REAL( ndastp )188 ENDIF189 178 ! 190 179 ! Write the information 191 CALL iom_rstput( kt, nitbkg_r, inum, 'rdastp' , zdate )192 193 180 IF ( ln_avgbkg ) THEN 194 CALL iom_rstput( kt, nitbkg_r, inum, 'un' , un_tavg ) 195 CALL iom_rstput( kt, nitbkg_r, inum, 'vn' , vn_tavg ) 196 CALL iom_rstput( kt, nitbkg_r, inum, 'tn' , tn_tavg ) 197 CALL iom_rstput( kt, nitbkg_r, inum, 'sn' , sn_tavg ) 198 CALL iom_rstput( kt, nitbkg_r, inum, 'sshn' , sshn_tavg) 199 #if defined key_zdftke 200 CALL iom_rstput( kt, nitbkg_r, inum, 'en' , en_tavg ) 201 #endif 202 CALL iom_rstput( kt, nitbkg_r, inum, 'gcx' , gcx_tavg) 203 CALL iom_rstput( kt, nitbkg_r, inum, 'avt' , avt_tavg) 181 IF( nitavgbkg_r == nit000 - 1 ) THEN ! Treat special case when nitavgbkg = 0 182 zdate = REAL( ndastp ) 183 #if defined key_zdftke 184 ! lk_zdftke=T : Read turbulent kinetic energy ( en ) 185 IF(lwp) WRITE(numout,*) ' Reading TKE (en) from restart...' 186 CALL tke_rst( nit000, 'READ' ) ! lk_zdftke=T : Read turbulent kinetic energy ( en ) 187 188 #endif 189 ELSE 190 zdate = REAL( ndastp ) 191 ENDIF 192 CALL iom_rstput( kt, nitavgbkg_r, inum, 'rdastp' , zdate ) 193 CALL iom_rstput( kt, nitavgbkg_r, inum, 'un' , un_tavg ) 194 CALL iom_rstput( kt, nitavgbkg_r, inum, 'vn' , vn_tavg ) 195 CALL iom_rstput( kt, nitavgbkg_r, inum, 'tn' , tn_tavg ) 196 CALL iom_rstput( kt, nitavgbkg_r, inum, 'sn' , sn_tavg ) 197 CALL iom_rstput( kt, nitavgbkg_r, inum, 'sshn' , sshn_tavg) 198 #if defined key_zdftke 199 CALL iom_rstput( kt, nitavgbkg_r, inum, 'en' , en_tavg ) 200 #endif 201 CALL iom_rstput( kt, nitavgbkg_r, inum, 'gcx' , gcx_tavg) 202 CALL iom_rstput( kt, nitavgbkg_r, inum, 'avt' , avt_tavg) 204 203 ! 205 204 ELSE 205 IF( nitbkg_r == nit000 - 1 ) THEN ! Treat special case when nitbkg = 0 206 zdate = REAL( ndastp ) 207 #if defined key_zdftke 208 ! lk_zdftke=T : Read turbulent kinetic energy ( en ) 209 IF(lwp) WRITE(numout,*) ' Reading TKE (en) from restart...' 210 CALL tke_rst( nit000, 'READ' ) ! lk_zdftke=T : Read turbulent kinetic energy ( en ) 211 212 #endif 213 ELSE 214 zdate = REAL( ndastp ) 215 ENDIF 216 CALL iom_rstput( kt, nitbkg_r, inum, 'rdastp' , zdate ) 206 217 CALL iom_rstput( kt, nitbkg_r, inum, 'un' , un ) 207 218 CALL iom_rstput( kt, nitbkg_r, inum, 'vn' , vn )
Note: See TracChangeset
for help on using the changeset viewer.