Changeset 13885
- Timestamp:
- 2020-11-26T15:01:11+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13787_OSMOSIS_IMMERSE/src/OCE/ZDF/zdfosm.F90
r13867 r13885 639 639 CALL zdf_osm_calculate_dhdt( zdhdt, zddhdt ) 640 640 DO_2D( 0, 0, 0, 0 ) 641 zhbl_t(ji,jj) = hbl(ji,jj) + (zdhdt(ji,jj) - ww(ji,jj,ibld(ji,jj)))* rn_Dt ! certainly need w nhere, so subtract it641 zhbl_t(ji,jj) = hbl(ji,jj) + (zdhdt(ji,jj) - ww(ji,jj,ibld(ji,jj)))* rn_Dt ! certainly need ww here, so subtract it 642 642 ! adjustment to represent limiting by ocean bottom 643 643 IF ( zhbl_t(ji,jj) >= gdepw(ji, jj, mbkt(ji,jj) + 1, Kmm ) ) THEN … … 2713 2713 IF( id1 > 0 ) THEN ! 'wn' exists; read 2714 2714 CALL iom_get( numror, jpdom_auto, 'wn', ww, ldxios = lrxios ) 2715 WRITE(numout,*) ' ===>>>> : w wread from restart file'2715 WRITE(numout,*) ' ===>>>> : wn read from restart file' 2716 2716 ELSE 2717 2717 ww(:,:,:) = 0._wp 2718 WRITE(numout,*) ' ===>>>> : w wnot in restart file, set to zero initially'2718 WRITE(numout,*) ' ===>>>> : wn not in restart file, set to zero initially' 2719 2719 END IF 2720 2720 … … 2746 2746 IF( TRIM(cdrw) == 'WRITE') THEN !* Write hbl into the restart file, then return 2747 2747 IF(lwp) WRITE(numout,*) '---- osm-rst ----' 2748 CALL iom_rstput( kt, nitrst, numrow, 'w w' , ww, ldxios = lwxios )2748 CALL iom_rstput( kt, nitrst, numrow, 'wn' , ww, ldxios = lwxios ) 2749 2749 CALL iom_rstput( kt, nitrst, numrow, 'hbl' , hbl, ldxios = lwxios ) 2750 2750 CALL iom_rstput( kt, nitrst, numrow, 'dh' , dh, ldxios = lwxios )
Note: See TracChangeset
for help on using the changeset viewer.