- Timestamp:
- 2020-01-14T13:40:47+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/ZDF/zdfosm.F90
r12312 r12317 2210 2210 CHARACTER(len=*), INTENT(in) :: cdrw ! "READ"/"WRITE" flag 2211 2211 2212 INTEGER :: id1, id2 ! iom enquiry index2212 INTEGER :: id1, id2, id3 ! iom enquiry index 2213 2213 INTEGER :: ji, jj, jk ! dummy loop indices 2214 2214 INTEGER :: iiki, ikt ! local integer … … 2238 2238 CALL iom_get( numror, jpdom_autoglo, 'dh', dh, ldxios = lrxios ) 2239 2239 WRITE(numout,*) ' ===>>>> : hbl & dh read from restart file' 2240 IF( ln_osm_mle ) THEN 2241 id3 = iom_varid( numror, 'hmle' , ldstop = .FALSE. ) 2242 IF( id3 > 0) THEN 2243 CALL iom_get( numror, jpdom_autoglo, 'hmle' , hmle , ldxios = lrxios ) 2244 WRITE(numout,*) ' ===>>>> : hmle read from restart file' 2245 ELSE 2246 WRITE(numout,*) ' ===>>>> : hmle not found, set to hbl' 2247 hmle(:,:) = hbl(:,:) ! Initialise MLE depth. 2248 END IF 2249 END IF 2240 2250 RETURN 2241 2251 ELSE ! 'hbl' & 'dh' not in restart file, recalculate … … 2249 2259 IF( TRIM(cdrw) == 'WRITE') THEN !* Write hbli into the restart file, then return 2250 2260 IF(lwp) WRITE(numout,*) '---- osm-rst ----' 2251 CALL iom_rstput( kt, nitrst, numrow, 'wn' , wn , ldxios = lwxios ) 2252 CALL iom_rstput( kt, nitrst, numrow, 'hbl' , hbl , ldxios = lwxios ) 2253 CALL iom_rstput( kt, nitrst, numrow, 'dh' , dh, ldxios = lwxios ) 2261 CALL iom_rstput( kt, nitrst, numrow, 'wn' , wn, ldxios = lwxios ) 2262 CALL iom_rstput( kt, nitrst, numrow, 'hbl' , hbl, ldxios = lwxios ) 2263 CALL iom_rstput( kt, nitrst, numrow, 'dh' , dh, ldxios = lwxios ) 2264 IF( ln_osm_mle ) THEN 2265 CALL iom_rstput( kt, nitrst, numrow, 'hmle', hmle, ldxios = lwxios ) 2266 END IF 2254 2267 RETURN 2255 2268 END IF … … 2285 2298 END DO 2286 2299 2287 IF( ln_osm_mle ) hmle(:,:) = hbl(:,:) ! Initialise MLE depth.2288 2289 2300 WRITE(numout,*) ' ===>>>> : hbl computed from stratification' 2301 2302 IF( ln_osm_mle ) THEN 2303 hmle(:,:) = hbl(:,:) ! Initialise MLE depth. 2304 WRITE(numout,*) ' ===>>>> : hmle set = to hbl' 2305 END IF 2306 2290 2307 wn(:,:,:) = 0._wp 2291 2308 WRITE(numout,*) ' ===>>>> : wn not in restart file, set to zero initially'
Note: See TracChangeset
for help on using the changeset viewer.