Changeset 2237 for LMDZ5/trunk/libf/phylmd/phyetat0.F90

Timestamp:

Mar 20, 2015, 7:53:17 AM (9 years ago)

Author:

fhourdin

Message:

Prise en compte des nouveaux alebedo dans les fichiers de redémarrage.
Retour à 1+1=2

Taking into account new albedos in restart files. 1+1=2

File:

• LMDZ5/trunk/libf/phylmd/phyetat0.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -11 +11 @@
   USE phys_state_var_mod, ONLY : ancien_ok, clwcon, detr_therm, dtime, &
        du_gwd_rando, dv_gwd_rando, entr_therm, f0, falb1, falb2, fm_therm, &
+       falb_dir, falb_dif, &
        ftsol, pbl_tke, pctsrf, q_ancien, radpas, radsol, rain_fall, ratqs, &
        rlat, rlon, rnebcon, rugoro, sig1, snow_fall, solaire_etat0, sollw, sollwdown, &
@@ -69 +70 @@
   INTEGER length
   PARAMETER (length=100)
-  INTEGER it, iiq
+  INTEGER it, iiq, isw
   REAL tab_cntrl(length), tabcntr0(length)
   CHARACTER*7 str7
@@ -243 +244 @@
   ENDIF
 
+!===================================================================
+  ! Lecture des albedo difus et direct
+
+  DO nsrf = 1, nbsrf
+     DO isw=1, nsw
+        IF (isw.GT.99 .AND. nsrf.GT.99) THEN
+           PRINT*, "Trop de bandes SW ou sous-mailles"
+           call abort_gcm("phyetat0", "", 1)
+        ENDIF
+        WRITE(str7, '(i2.2, "srf", i2.2)') isw, nsrf
+
+        CALL get_field('A_dir_SW'//str7, falb_dir(:, isw, nsrf), found)
+        IF (.NOT. found) THEN
+           PRINT*, "phyetat0: Le champ <A_dir_SW"//str7//"> est absent"
+           PRINT*, "          Il prend donc la valeur de surface"
+           DO i=1, klon
+              falb_dir(i, isw, nsrf)=0.2
+           ENDDO
+        ENDIF
+        CALL get_field('A_dif_SW'//str7, falb_dif(:, isw, nsrf), found)
+        IF (.NOT. found) THEN
+           PRINT*, "phyetat0: Le champ <A_dif_SW"//str7//"> est absent"
+           PRINT*, "          Il prend donc la valeur de surface"
+           DO i=1, klon
+              falb_dif(i, isw, nsrf)=0.2
+           ENDDO
+        ENDIF
+     ENDDO
+  ENDDO
+
+!===================================================================
   ! Lecture des temperatures du sol profond:
 
@@ -264 +296 @@
   ENDDO
 
+!===================================================================
   ! Lecture de l'humidite de l'air juste au dessus du sol: