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

Timestamp:

Jun 20, 2014, 12:27:40 PM (10 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)

File:

• LMDZ5/trunk/libf/phylmd/phyetat0.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