Changeset 12317


Ignore:
Timestamp:
2020-01-14T13:40:47+01:00 (2 weeks ago)
Author:
agn
Message:

hmle is now restartable

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/ZDF/zdfosm.F90

    r12312 r12317  
    22102210     CHARACTER(len=*), INTENT(in) ::   cdrw   ! "READ"/"WRITE" flag 
    22112211 
    2212      INTEGER ::   id1, id2   ! iom enquiry index 
     2212     INTEGER ::   id1, id2, id3   ! iom enquiry index 
    22132213     INTEGER  ::   ji, jj, jk     ! dummy loop indices 
    22142214     INTEGER  ::   iiki, ikt ! local integer 
     
    22382238           CALL iom_get( numror, jpdom_autoglo, 'dh', dh, ldxios = lrxios  ) 
    22392239           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 
    22402250           RETURN 
    22412251        ELSE                      ! 'hbl' & 'dh' not in restart file, recalculate 
     
    22492259     IF( TRIM(cdrw) == 'WRITE') THEN     !* Write hbli into the restart file, then return 
    22502260        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 
    22542267        RETURN 
    22552268     END IF 
     
    22852298     END DO 
    22862299 
    2287      IF( ln_osm_mle ) hmle(:,:) = hbl(:,:)            ! Initialise MLE depth. 
    2288  
    22892300     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 
    22902307     wn(:,:,:) = 0._wp 
    22912308     WRITE(numout,*) ' ===>>>> :  wn not in restart file, set to zero initially' 
Note: See TracChangeset for help on using the changeset viewer.