Changeset 2237

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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• lmdz1d.F90 (modified) (1 diff)
• phyetat0.F90 (modified) (4 diffs)
• phyredem.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/lmdz1d.F90

@@ -668 +668 @@
         zpic = zpicinp
         ftsol=tsurf
+        nsw=6 ! on met le nb de bandes SW=6, pour initialiser
+              ! 6 albedo, mais on peut quand meme tourner avec
+              ! moins. Seules les 2 ou 4 premiers seront lus
         falb1 = albedo                            
         falb2 = albedo                            
+        falb_dir=albedo
+        falb_dif=albedo
         rugoro=rugos 
         t_ancien(1,:)=temp(:)

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:
 

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -51 +51 @@
   REAL tab_cntrl(length)
 
-  INTEGER isoil, nsrf
+  INTEGER isoil, nsrf,isw
   CHARACTER (len=7) :: str7
   CHARACTER (len=2) :: str2
@@ -142 +142 @@
   ENDDO
 
+! ================== Albedo =======================================
+  print*,'PHYREDEM NOUVEAU'
+  DO nsrf = 1, nbsrf
+     DO isw=1, nsw
+        IF (isw.LE.99 .AND. nsrf.LE.99) THEN
+           WRITE(str7, '(i2.2, "srf", i2.2)') isw, nsrf
+  print*,'PHYREDEM ',"A_dir_SW"//str7
+           CALL put_field("A_dir_SW"//str7, "Albedo direct du sol bande "//str7, &
+                falb_dir(:, isw, nsrf))
+           CALL put_field("A_dif_SW"//str7, "Albedo difus du sol bande "//str7, &
+                falb_dif(:, isw, nsrf))
+        ELSE
+           PRINT*, "Trop de couches"
+           call abort_gcm("phyredem", "", 1)
+        ENDIF
+     ENDDO
+  ENDDO
+
+! ================== Tsoil =======================================
   DO nsrf = 1, nbsrf
      DO isoil=1, nsoilmx