Changeset 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_do.F90

Timestamp:

2017-12-26T17:32:56+01:00 (7 years ago)

Author:

gm

Message:

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_do.F90 (modified) (10 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_do.F90

@@ -39 +39 @@
    PUBLIC   ice_thd_do_init   ! called by ice_stp
 
-   ! ** namelist (namthd_do) **
-   REAL(wp) ::   rn_hinew         ! thickness for new ice formation (m)
-   LOGICAL  ::   ln_frazil        ! use of frazil ice collection as function of wind (T) or not (F)
-   REAL(wp) ::   rn_maxfraz       ! maximum portion of frazil ice collecting at the ice bottom
-   REAL(wp) ::   rn_vfraz         ! threshold drift speed for collection of bottom frazil ice
-   REAL(wp) ::   rn_Cfraz         ! squeezing coefficient for collection of bottom frazil ice
+   !                          !!** namelist (namthd_do) **
+   REAL(wp) ::   rn_hinew      ! thickness for new ice formation (m)
+   LOGICAL  ::   ln_frazil     ! use of frazil ice collection as function of wind (T) or not (F)
+   REAL(wp) ::   rn_maxfraz    ! maximum portion of frazil ice collecting at the ice bottom
+   REAL(wp) ::   rn_vfraz      ! threshold drift speed for collection of bottom frazil ice
+   REAL(wp) ::   rn_Cfraz      ! squeezing coefficient for collection of bottom frazil ice
 
    !!----------------------------------------------------------------------
@@ -78 +78 @@
       !!               update h_s_1d, h_i_1d      
       !!------------------------------------------------------------------------
-      INTEGER  ::   ji,jj,jk,jl      ! dummy loop indices
+      INTEGER  ::   ji, jj, jk, jl   ! dummy loop indices
       INTEGER  ::   iter     !   -       -
       REAL(wp) ::   ztmelts, zfrazb, zweight, zde                          ! local scalars
       REAL(wp) ::   zgamafr, zvfrx, zvgx, ztaux, ztwogp, zf                     !   -      -
       REAL(wp) ::   ztenagm, zvfry, zvgy, ztauy, zvrel2, zfp, zsqcd , zhicrit   !   -      -
-
+      !
       REAL(wp) ::   zQm          ! enthalpy exchanged with the ocean (J/m2, >0 towards ocean)
       REAL(wp) ::   zEi          ! sea ice specific enthalpy (J/kg)
       REAL(wp) ::   zEw          ! seawater specific enthalpy (J/kg)
       REAL(wp) ::   zfmdt        ! mass flux x time step (kg/m2, >0 towards ocean)
-     
+      !
       REAL(wp) ::   zv_newfra
-  
+      !
       INTEGER , DIMENSION(jpij) ::   jcat        ! indexes of categories where new ice grows
       REAL(wp), DIMENSION(jpij) ::   zswinew     ! switch for new ice or not
-
+      !
       REAL(wp), DIMENSION(jpij) ::   zv_newice   ! volume of accreted ice
       REAL(wp), DIMENSION(jpij) ::   za_newice   ! fractional area of accreted ice
@@ -104 +104 @@
       REAL(wp), DIMENSION(jpij) ::   zv_frazb    ! accretion of frazil ice at the ice bottom
       REAL(wp), DIMENSION(jpij) ::   zvrel_1d    ! relative ice / frazil velocity (1D vector)
-
+      !
       REAL(wp), DIMENSION(jpij,jpl) ::   zv_b      ! old volume of ice in category jl
       REAL(wp), DIMENSION(jpij,jpl) ::   za_b      ! old area of ice in category jl
-
+      !
       REAL(wp), DIMENSION(jpij,nlay_i,jpl) ::   ze_i_2d !: 1-D version of e_i
-
+      !
       REAL(wp), DIMENSION(jpi,jpj) ::   zvrel     ! relative ice / frazil velocity
-
+      !
       REAL(wp) :: zcai = 1.4e-3_wp                     ! ice-air drag (clem: should be dependent on coupling/forcing used)
       !!-----------------------------------------------------------------------!
 
-      IF( ln_icediachk )   CALL ice_cons_hsm(0, 'icethd_do', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft)
+      IF( ln_icediachk )   CALL ice_cons_hsm( 0, 'icethd_do', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft )
 
       CALL ice_var_agg(1)
@@ -141 +141 @@
 
       IF( ln_frazil ) THEN
-
+         !
          !--------------------
          ! Physical constants
          !--------------------
          ht_i_new(:,:) = 0._wp
-
+         !
          zhicrit = 0.04 ! frazil ice thickness
          ztwogp  = 2. * rau0 / ( grav * 0.3 * ( rau0 - rhoic ) ) ! reduced grav
          zsqcd   = 1.0 / SQRT( 1.3 * zcai ) ! 1/SQRT(airdensity*drag)
          zgamafr = 0.03
-
+         !
          DO jj = 2, jpjm1
             DO ji = 2, jpim1
@@ -204 +204 @@
                      iter = iter + 1
                   END DO
-
+                  !
                ENDIF ! end of selection of pixels where ice forms
-
+               !
             END DO 
          END DO 
@@ -222 +222 @@
       !-------------------------------------
       ! This occurs if open water energy budget is negative (cooling) and there is no landfast ice
-      npti = 0 ; nptidx(:) = 0
+      npti = 0   ;   nptidx(:) = 0
       DO jj = 1, jpj
          DO ji = 1, jpi
@@ -431 +431 @@
                h_i_old (ji,nlay_i+1) = zv_newfra
                eh_i_old(ji,nlay_i+1) = ze_newice(ji) * zv_newfra
-            ENDDO
+            END DO
             ! --- Ice enthalpy remapping --- !
             CALL ice_thd_ent( ze_i_2d(1:npti,:,jl) ) 
-         ENDDO
+         END DO
 
          !-----------------
@@ -474 +474 @@
       !
    END SUBROUTINE ice_thd_do
+
 
    SUBROUTINE ice_thd_do_init
@@ -487 +488 @@
       !! ** input   :   Namelist namthd_do
       !!-------------------------------------------------------------------
-      INTEGER  ::   ios   ! Local integer output status for namelist read
+      INTEGER  ::   ios   ! Local integer 
       !!
       NAMELIST/namthd_do/ rn_hinew, ln_frazil, rn_maxfraz, rn_vfraz, rn_Cfraz
@@ -494 +495 @@
       REWIND( numnam_ice_ref )              ! Namelist namthd_do in reference namelist : Ice thermodynamics
       READ  ( numnam_ice_ref, namthd_do, IOSTAT = ios, ERR = 901)
-901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namthd_do in reference namelist', lwp )
-
+901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namthd_do in reference namelist', lwp )
       REWIND( numnam_ice_cfg )              ! Namelist namthd_do in configuration namelist : Ice thermodynamics
       READ  ( numnam_ice_cfg, namthd_do, IOSTAT = ios, ERR = 902 )
-902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namthd_do in configuration namelist', lwp )
-      IF(lwm) WRITE ( numoni, namthd_do )
-      !
+902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namthd_do in configuration namelist', lwp )
+      IF(lwm) WRITE( numoni, namthd_do )
       !
       IF(lwp) THEN                          ! control print
+         WRITE(numout,*)
          WRITE(numout,*) 'ice_thd_do_init: Ice growth in open water'
          WRITE(numout,*) '~~~~~~~~~~~~~~~'