Changeset 8505 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90

Timestamp:

2017-09-06T19:02:43+02:00 (7 years ago)

Author:

clem

Message:

changes in style - part5 - start changing init routines

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90 (modified) (7 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90

@@ -13 +13 @@
    !!   alb_init      : initialisation of albedo computation
    !!----------------------------------------------------------------------
-   USE ice     , ONLY : jpl   ! number of ice category
+   USE ice, ONLY: jpl ! sea-ice: number of categories
    USE phycst         ! physical constants
    !
@@ -23 +23 @@
    PRIVATE
 
-   PUBLIC   ice_alb   ! routine called in iceforcing.F90 and iceupdate.F90
+   PUBLIC   ice_alb_init   ! called in icestp
+   PUBLIC   ice_alb        ! called in iceforcing.F90 and iceupdate.F90
 
    REAL(wp), PUBLIC, PARAMETER ::   rn_alb_oce = 0.066   !: ocean or lead albedo (Pegau and Paulson, Ann. Glac. 2001)
 
-   INTEGER  ::   albd_init = 0       ! control flag for initialization  
    REAL(wp) , PARAMETER ::   rc1    = 0.05    ! snow thickness (only for nn_ice_alb=0)
    REAL(wp) , PARAMETER ::   rc2    = 0.10    !  "        "
@@ -33 +33 @@
    REAL(wp) , PARAMETER ::   r1_c1 = 1. / rc1
    REAL(wp) , PARAMETER ::   r1_c2 = 1. / rc2
- 
-   !                             !!* namelist namsbc_alb *
-   INTEGER  ::   nn_ice_alb
-   REAL(wp) ::   rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt, rn_alb_dpnd
+   !
+   ! ** albedo namelist (namicealb)
+   INTEGER  ::   nn_ice_alb       ! type of albedo scheme: 0: Shine & Henderson-Sellers (JGR 1985)
+   !                                      !                         1: "home made" based on Brandt et al. (JClim 2005)
+   !                                      !                             and Grenfell & Perovich (JGR 2004)
+   !                                      !                         2: Same as 1 but with melt ponds
+   REAL(wp) ::   rn_alb_sdry      ! dry snow albedo
+   REAL(wp) ::   rn_alb_smlt      ! melting snow albedo
+   REAL(wp) ::   rn_alb_idry      ! dry ice albedo
+   REAL(wp) ::   rn_alb_imlt      ! bare puddled ice albedo
+   REAL(wp) ::   rn_alb_dpnd      ! ponded ice albedo
 
    !!----------------------------------------------------------------------
@@ -102 +109 @@
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zafrac_pnd      ! relative ice fraction (effective)
       !!---------------------------------------------------------------------
-
-      IF( albd_init == 0 )   CALL alb_init      ! initialization 
 
       !-----------------------------------------------------
@@ -299 +304 @@
 
 
-   SUBROUTINE alb_init
+   SUBROUTINE ice_alb_init
       !!----------------------------------------------------------------------
       !!                 ***  ROUTINE alb_init  ***
@@ -311 +316 @@
       NAMELIST/namicealb/ nn_ice_alb, rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt, rn_alb_dpnd
       !!----------------------------------------------------------------------
-      !
-      albd_init = 1                     ! indicate that the initialization has been done
       !
       REWIND( numnam_ice_ref )              ! Namelist namicealb in reference namelist : Albedo parameters
@@ -336 +339 @@
       ENDIF
       !
-   END SUBROUTINE alb_init
+   END SUBROUTINE ice_alb_init
 
 #else