Changeset 2728

Timestamp:

Dec 13, 2016, 10:03:48 AM (8 years ago)

Author:

fhourdin

Message:

Ajout de la routine tend_to_tke qui transforme des tendances
physiques de grande echelle en source d'energie cinetique
turbulent TKE

Externalisation de parametres d'albedo pour la neige pour
les calottes :
alb_vis_sno_lic, alb_nir_sno_lic = 0.77 par défaut,
valeur volontairement sousestimee au moment de CMIP3 pour
rechauffer les calottes.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• surf_landice_mod.F90 (modified) (4 diffs)
• tend_to_tke.F90 (added)

LMDZ5/trunk/libf/phylmd/surf_landice_mod.F90

@@ -30 +30 @@
     USE calcul_fluxs_mod
     USE phys_output_var_mod
+!FC
+    USE ioipsl_getin_p_mod, ONLY : getin_p
+
 #ifdef CPP_SISVAT
     USE surf_sisvat_mod,  ONLY : surf_sisvat
@@ -121 +124 @@
     CHARACTER (len = 80)                      :: abort_message
 
+
 !albedo SB >>>
     real,dimension(klon) :: alb1,alb2
@@ -127 +131 @@
 ! End definition
 !****************************************************************************************
+!FC 
+!FC
+   REAL,SAVE :: alb_vis_sno_lic
+  !$OMP THREADPRIVATE(alb_vis_sno_lic)
+   REAL,SAVE :: alb_nir_sno_lic
+  !$OMP THREADPRIVATE(alb_nir_sno_lic)
+  LOGICAL, SAVE :: firstcall = .TRUE.
+  !$OMP THREADPRIVATE(firstcall)
+!FC
+
+
+  IF (firstcall) THEN
+  alb_vis_sno_lic=0.77
+  CALL getin_p('alb_vis_sno_lic',alb_vis_sno_lic)
+           PRINT*, 'alb_vis_sno_lic',alb_vis_sno_lic
+  alb_nir_sno_lic=0.77
+  CALL getin_p('alb_nir_sno_lic',alb_nir_sno_lic)
+           PRINT*, 'alb_nir_sno_lic',alb_nir_sno_lic
+  firstcall=.false.
+  ENDIF
 !
 ! Initialize output variables
@@ -280 +304 @@
 
 ! Attantion: alb1 and alb2 are the same!
-    alb1(1:knon)  = 0.77
-    alb2(1:knon)  = alb1(1:knon)
+    alb1(1:knon)  = alb_vis_sno_lic
+    alb2(1:knon)  = alb_nir_sno_lic