Changeset 14086 for NEMO/trunk/src/ICE

Timestamp:

2020-12-04T12:37:14+01:00 (3 years ago)

Author:

cetlod

Message:

Adding AGRIF branches into the trunk

File:

• NEMO/trunk/src/ICE/iceistate.F90 (modified) (19 diffs)

NEMO/trunk/src/ICE/iceistate.F90

@@ -18 +18 @@
    USE oce            ! dynamics and tracers variables
    USE dom_oce        ! ocean domain
-   USE sbc_oce , ONLY : sst_m, sss_m, ln_ice_embd
+   USE sbc_oce , ONLY : sst_m, sss_m, ln_ice_embd 
    USE sbc_ice , ONLY : tn_ice, snwice_mass, snwice_mass_b
    USE eosbn2         ! equation of state
@@ -39 +39 @@
 # if defined key_agrif
    USE agrif_oce
-   USE agrif_ice
-   USE agrif_ice_interp
-# endif
+   USE agrif_ice_interp 
+# endif   
 
    IMPLICIT NONE
@@ -91 +90 @@
       !!
       !! ** Method  :   This routine will put some ice where ocean
-      !!                is at the freezing point, then fill in ice
-      !!                state variables using prescribed initial
-      !!                values in the namelist
+      !!                is at the freezing point, then fill in ice 
+      !!                state variables using prescribed initial 
+      !!                values in the namelist            
       !!
       !! ** Steps   :   1) Set initial surface and basal temperatures
@@ -103 +102 @@
       !!              where there is no ice
       !!--------------------------------------------------------------------
-      INTEGER, INTENT(in) :: kt            ! time step
+      INTEGER, INTENT(in) :: kt            ! time step 
       INTEGER, INTENT(in) :: Kbb, Kmm, Kaa ! ocean time level indices
       !
@@ -129 +128 @@
       ! basal temperature (considered at freezing point)   [Kelvin]
       CALL eos_fzp( sss_m(:,:), t_bo(:,:) )
-      t_bo(:,:) = ( t_bo(:,:) + rt0 ) * tmask(:,:,1)
+      t_bo(:,:) = ( t_bo(:,:) + rt0 ) * tmask(:,:,1) 
       !
       ! surface temperature and conductivity
@@ -154 +153 @@
       e_i (:,:,:,:) = 0._wp
       e_s (:,:,:,:) = 0._wp
-
+      
       ! general fields
       a_i (:,:,:) = 0._wp
@@ -184 +183 @@
       IF( ln_iceini ) THEN
          !
-         IF( Agrif_Root() ) THEN
+#if defined key_agrif
+         IF ( ( Agrif_Root() ).OR.(.NOT.ln_init_chfrpar ) ) THEN
+#endif
             !                             !---------------!
             IF( nn_iceini_file == 1 )THEN ! Read a file   !
@@ -229 +230 @@
                IF( TRIM(si(jp_apd)%clrootname) == 'NOT USED' ) &
                   &     si(jp_apd)%fnow(:,:,1) = ( rn_apd_ini_n * zswitch + rn_apd_ini_s * (1._wp - zswitch) ) * tmask(:,:,1) & ! rn_apd = pond fraction => rn_apnd * a_i = pond conc.
-                  &                              * si(jp_ati)%fnow(:,:,1)
+                  &                              * si(jp_ati)%fnow(:,:,1) 
                !
                ! pond depth
@@ -248 +249 @@
                !
                ! change the switch for the following
-               WHERE( zat_i_ini(:,:) > 0._wp )   ;   zswitch(:,:) = tmask(:,:,1)
+               WHERE( zat_i_ini(:,:) > 0._wp )   ;   zswitch(:,:) = tmask(:,:,1) 
                ELSEWHERE                         ;   zswitch(:,:) = 0._wp
                END WHERE
@@ -256 +257 @@
                !                          !---------------!
                ! no ice if (sst - Tfreez) >= thresold
-               WHERE( ( sst_m(:,:) - (t_bo(:,:) - rt0) ) * tmask(:,:,1) >= rn_thres_sst )   ;   zswitch(:,:) = 0._wp
+               WHERE( ( sst_m(:,:) - (t_bo(:,:) - rt0) ) * tmask(:,:,1) >= rn_thres_sst )   ;   zswitch(:,:) = 0._wp 
                ELSEWHERE                                                                    ;   zswitch(:,:) = tmask(:,:,1)
                END WHERE
@@ -269 +270 @@
                   zt_su_ini(:,:) = rn_tsu_ini_n * zswitch(:,:)
                   ztm_s_ini(:,:) = rn_tms_ini_n * zswitch(:,:)
-                  zapnd_ini(:,:) = rn_apd_ini_n * zswitch(:,:) * zat_i_ini(:,:) ! rn_apd = pond fraction => rn_apd * a_i = pond conc.
+                  zapnd_ini(:,:) = rn_apd_ini_n * zswitch(:,:) * zat_i_ini(:,:) ! rn_apd = pond fraction => rn_apd * a_i = pond conc. 
                   zhpnd_ini(:,:) = rn_hpd_ini_n * zswitch(:,:)
                   zhlid_ini(:,:) = rn_hld_ini_n * zswitch(:,:)
@@ -295 +296 @@
                zhlid_ini(:,:) = 0._wp
             ENDIF
-
+            
             IF ( .NOT.ln_pnd_lids ) THEN
                zhlid_ini(:,:) = 0._wp
             ENDIF
-
+            
             !----------------!
             ! 3) fill fields !
@@ -323 +324 @@
             CALL tab_2d_1d( npti, nptidx(1:npti), h_ip_1d(1:npti)  , zhpnd_ini )
             CALL tab_2d_1d( npti, nptidx(1:npti), h_il_1d(1:npti)  , zhlid_ini )
-
+            
             ! allocate temporary arrays
             ALLOCATE( zhi_2d (npti,jpl), zhs_2d (npti,jpl), zai_2d (npti,jpl), &
@@ -377 +378 @@
             DO jl = 1, jpl
                DO_3D( 1, 1, 1, 1, 1, nlay_i )
-                  t_i (ji,jj,jk,jl) = zti_3d(ji,jj,jl)
+                  t_i (ji,jj,jk,jl) = zti_3d(ji,jj,jl) 
                   ztmelts          = - rTmlt * sz_i(ji,jj,jk,jl) + rt0 ! melting temperature in K
                   e_i(ji,jj,jk,jl) = zswitch(ji,jj) * v_i(ji,jj,jl) * r1_nlay_i * &
@@ -385 +386 @@
                END_3D
             END DO
-
-#if  defined key_agrif
+            !
+#if defined key_agrif
          ELSE
-
-            Agrif_SpecialValue    = -9999.
-            Agrif_UseSpecialValue = .TRUE.
-            CALL Agrif_init_variable(tra_iceini_id,procname=interp_tra_ice)
-            use_sign_north = .TRUE.
-            sign_north = -1.
-            CALL Agrif_init_variable(u_iceini_id  ,procname=interp_u_ice)
-            CALL Agrif_init_variable(v_iceini_id  ,procname=interp_v_ice)
-            Agrif_SpecialValue    = 0._wp
-            use_sign_north = .FALSE.
-            Agrif_UseSpecialValue = .FALSE.
-        ! lbc ????
-   ! Here we know : a_i, v_i, v_s, sv_i, oa_i, a_ip, v_ip, v_il, t_su, e_s, e_i
-            CALL ice_var_glo2eqv
-            CALL ice_var_zapsmall
-            CALL ice_var_agg(2)
+            CALL  agrif_istate_ice
+         ENDIF
 #endif
-         ENDIF ! Agrif_Root
-         !
          ! Melt ponds
          WHERE( a_i > epsi10 )   ;   a_ip_eff(:,:,:) = a_ip(:,:,:) / a_i(:,:,:)
@@ -413 +398 @@
          v_ip(:,:,:) = h_ip(:,:,:) * a_ip(:,:,:)
          v_il(:,:,:) = h_il(:,:,:) * a_ip(:,:,:)
-
+         
          ! specific temperatures for coupled runs
          tn_ice(:,:,:) = t_su(:,:,:)
@@ -423 +408 @@
             WHERE( at_i(:,:) > rn_amax_2d(:,:) )   a_i(:,:,jl) = a_i(:,:,jl) * rn_amax_2d(:,:) / at_i(:,:)
          END DO
-         at_i(:,:) = SUM( a_i, dim=3 )
+        at_i(:,:) = SUM( a_i, dim=3 )
          !
       ENDIF ! ln_iceini
@@ -456 +441 @@
       !!-------------------------------------------------------------------
       !!                   ***  ROUTINE ice_istate_init  ***
-      !!
-      !! ** Purpose :   Definition of initial state of the ice
-      !!
-      !! ** Method  :   Read the namini namelist and check the parameter
+      !!        
+      !! ** Purpose :   Definition of initial state of the ice 
+      !!
+      !! ** Method  :   Read the namini namelist and check the parameter 
       !!              values called at the first timestep (nit000)
       !!
@@ -500 +485 @@
          WRITE(numout,*) '      max ocean temp. above Tfreeze with initial ice   rn_thres_sst   = ', rn_thres_sst
          IF( ln_iceini .AND. nn_iceini_file == 0 ) THEN
-            WRITE(numout,*) '      initial snw thickness in the north-south         rn_hts_ini     = ', rn_hts_ini_n,rn_hts_ini_s
+            WRITE(numout,*) '      initial snw thickness in the north-south         rn_hts_ini     = ', rn_hts_ini_n,rn_hts_ini_s 
             WRITE(numout,*) '      initial ice thickness in the north-south         rn_hti_ini     = ', rn_hti_ini_n,rn_hti_ini_s
             WRITE(numout,*) '      initial ice concentr  in the north-south         rn_ati_ini     = ', rn_ati_ini_n,rn_ati_ini_s