Changeset 2069 for LMDZ5/trunk/libf

Timestamp:

Jun 20, 2014, 12:27:40 PM (11 years ago)

Author:

Laurent Fairhead

Message:

Modifications pour avoir 1+1=2 avec la nouvelle physique (testé avec NPv4.12)

---

Modifications needed to ensure 1+1=2 with the new physics (tested with NPv4.12)

Location:

LMDZ5/trunk/libf/phylmd

Files:

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

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -15 +15 @@
        solsw, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
        wake_deltat, wake_fip, wake_pe, wake_s, zgam, zmax0, zmea, zpic, zsig, &
-       zstd, zthe, zval
+       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
   USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
   USE infotrac, only: nbtr, type_trac, tname, niadv
@@ -997 +997 @@
   PRINT*, '(ecart-type) detr_therm:', xmin, xmax
 
+  CALL get_field("ALE_BL", ale_bl, found)
+  IF (.NOT. found) THEN
+     PRINT*, "phyetat0: Le champ <ale_bl> est absent"
+     PRINT*, "Depart legerement fausse. Mais je continue"
+     ale_bl=0.
+  ENDIF
+  xmin = 1.0E+20
+  xmax = -1.0E+20
+  xmin = MINval(ale_bl)
+  xmax = MAXval(ale_bl)
+  PRINT*, '(ecart-type) ale_bl:', xmin, xmax
+
+  CALL get_field("ALE_BL_TRIG", ale_bl_trig, found)
+  IF (.NOT. found) THEN
+     PRINT*, "phyetat0: Le champ <ale_bl_trig> est absent"
+     PRINT*, "Depart legerement fausse. Mais je continue"
+     ale_bl_trig=0.
+  ENDIF
+  xmin = 1.0E+20
+  xmax = -1.0E+20
+  xmin = MINval(ale_bl_trig)
+  xmax = MAXval(ale_bl_trig)
+  PRINT*, '(ecart-type) ale_bl_trig:', xmin, xmax
+
+  CALL get_field("ALP_BL", alp_bl, found)
+  IF (.NOT. found) THEN
+     PRINT*, "phyetat0: Le champ <alp_bl> est absent"
+     PRINT*, "Depart legerement fausse. Mais je continue"
+     alp_bl=0.
+  ENDIF
+  xmin = 1.0E+20
+  xmax = -1.0E+20
+  xmin = MINval(alp_bl)
+  xmax = MAXval(alp_bl)
+  PRINT*, '(ecart-type) alp_bl:', xmin, xmax
+
   ! Read and send field trs to traclmdz
 

LMDZ5/trunk/libf/phylmd/phyredem.F90

@@ -323 +323 @@
 
   CALL put_field("DETR_THERM", "DETR_THERM", detr_therm)
+
+  CALL put_field("ALE_BL", "ALE_BL", Ale_bl)
+
+  CALL put_field("ALE_BL_TRIG", "ALE_BL_TRIG", Ale_bl_trig)
+
+  CALL put_field("ALP_BL", "ALP_BL", Alp_bl)
 
   ! trs from traclmdz_mod

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1121 +1121 @@
         endif
 
-        do i = 1,klon
-           Ale_bl(i)=0.
-           Alp_bl(i)=0.
-        enddo
+!        do i = 1,klon
+!           Ale_bl(i)=0.
+!           Alp_bl(i)=0.
+!        enddo
 
         !================================================================================
@@ -1986 +1986 @@
         enddo
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
 
      endif