Changeset 11858 for NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/ABL/sbcabl.F90

Timestamp:

2019-11-05T15:32:44+01:00 (5 years ago)

Author:

gsamson

Message:

dev_r11265_ABL: add write/read restart options for abl (#2131)

File:

• NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/ABL/sbcabl.F90 (modified) (4 diffs)

NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/src/ABL/sbcabl.F90

@@ -16 +16 @@
    USE par_abl        ! abl parameters
    USE ablmod
+   USE ablrst
 
    USE phycst         ! physical constants
@@ -67 +68 @@
       REAL(wp)           ::   zcff,zcff1
       LOGICAL            ::   lluldl
-      NAMELIST/namsbc_abl/ cn_dir , cn_dom, ln_hpgls_frc, ln_geos_winds,          &
-         &                 nn_dyn_restore,                                        &
+      NAMELIST/namsbc_abl/ cn_dir, cn_dom, cn_ablrst_in, cn_ablrst_out,           &
+         &                 cn_ablrst_indir, cn_ablrst_outdir,                     &
+         &                 ln_hpgls_frc, ln_geos_winds, nn_dyn_restore,           &
          &                 rn_ldyn_min , rn_ldyn_max, rn_ltra_min, rn_ltra_max,   &
          &                 nn_amxl, rn_cm, rn_ct, rn_ce, rn_ceps, rn_Rod, rn_Ric, &
@@ -267 +269 @@
       CALL sbc_blk_init
 
-      ! initialize ABL from data or restart
-      !!GS  disabled for now
-      !IF( ln_rstart ) THEN
-      !  CALL ctl_stop( 'STOP', 'sbc_abl_init: restart mode not supported yet' )
-      !ELSE
-
-      CALL fld_read( nit000, nn_fsbc, sf ) ! input fields provided at the first time-step
-
       ! Initialize the time index for now time (nt_n) and after time (nt_a)
       nt_n = 1 + MOD( nit000  , 2)
       nt_a = 1 + MOD( nit000+1, 2)
 
-       u_abl(:,:,:,nt_n      ) = sf(jp_wndi)%fnow(:,:,:)
-       v_abl(:,:,:,nt_n      ) = sf(jp_wndj)%fnow(:,:,:)
-      tq_abl(:,:,:,nt_n,jp_ta) = sf(jp_tair)%fnow(:,:,:)
-      tq_abl(:,:,:,nt_n,jp_qa) = sf(jp_humi)%fnow(:,:,:)
-
-      tke_abl(:,:,:,nt_n     ) = tke_min
-      avm_abl(:,:,:          ) = avm_bak
-      avt_abl(:,:,:          ) = avt_bak
-      mxl_abl(:,:,:          ) = mxl_min 
-      pblh   (:,:            ) = ghw_abl( 3 )  !<-- assume that the pbl contains 3 grid points 
-      u_abl  (:,:,:,nt_a     ) = 0._wp
-      v_abl  (:,:,:,nt_a     ) = 0._wp
-      tq_abl (:,:,:,nt_a,:   ) = 0._wp
-      tke_abl(:,:,:,nt_a     ) = 0._wp
-      !ENDIF
-      !!GS restart case not supported
+      ! initialize ABL from data or restart
+      IF( ln_rstart ) THEN
+         CALL abl_rst_read
+      ELSE
+         CALL fld_read( nit000, nn_fsbc, sf ) ! input fields provided at the first time-step
+
+          u_abl(:,:,:,nt_n      ) = sf(jp_wndi)%fnow(:,:,:)
+          v_abl(:,:,:,nt_n      ) = sf(jp_wndj)%fnow(:,:,:)
+         tq_abl(:,:,:,nt_n,jp_ta) = sf(jp_tair)%fnow(:,:,:)
+         tq_abl(:,:,:,nt_n,jp_qa) = sf(jp_humi)%fnow(:,:,:)
+   
+         tke_abl(:,:,:,nt_n     ) = tke_min
+         avm_abl(:,:,:          ) = avm_bak
+         avt_abl(:,:,:          ) = avt_bak
+         mxl_abl(:,:,:          ) = mxl_min 
+         pblh   (:,:            ) = ghw_abl( 3 )  !<-- assume that the pbl contains 3 grid points 
+         u_abl  (:,:,:,nt_a     ) = 0._wp
+         v_abl  (:,:,:,nt_a     ) = 0._wp
+         tq_abl (:,:,:,nt_a,:   ) = 0._wp
+         tke_abl(:,:,:,nt_a     ) = 0._wp
+      ENDIF
      
    END SUBROUTINE sbc_abl_init
 
-  
 
    SUBROUTINE sbc_abl( kt )
@@ -379 +377 @@
          &            sst_m, zsen, zevp                                )
 
+      CALL abl_rst_opn( kt )                       ! Open abl restart file (if necessary) 
+      IF( lrst_abl ) CALL abl_rst_write( kt )      ! -- abl restart file 
+
 #if defined key_si3
       ! Avoid a USE abl in icesbc module