Changeset 2075

Timestamp:

Jul 3, 2014, 2:27:38 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Paramétrisation de la fréquence de couplage océan/atmosphère

---

Parametrising the ocean/atmosphere coupling timestep

Location:

LMDZ5/trunk/libf/phylmd

Files:

• conf_phys_m.F90 (modified) (3 diffs)
• cpl_mod.F90 (modified) (1 diff)
• surface_data.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -99 +99 @@
     real,SAVE           :: ratqshaut_omp
     real,SAVE           :: tau_ratqs_omp
+    REAL, SAVE          :: t_coupl_omp
     integer,SAVE        :: iflag_radia_omp
     integer,SAVE        :: iflag_rrtm_omp
@@ -228 +229 @@
        WRITE(lunout,*)'You have to update your parameter file physiq.def to succed running'
        CALL abort_gcm('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1)
+    END IF
+
+    !Config Key  = t_coupl
+    !Config Desc = Pas de temps du couplage atm/oce en sec.
+    !Config Def  = 86400
+    !Config Help = This is only for testing purpose
+    !
+    t_coupl_omp = 86400.
+    call getin('t_coupl', t_coupl_omp)
+    IF (t_coupl_omp == 0) THEN
+       WRITE(lunout,*)'ERROR! Timestep of coupling between atmosphere and ocean'
+       WRITE(lunout,*)'cannot be zero.'
+       CALL abort_gcm('conf_phys','t_coupl = 0.',1)
     END IF
 
@@ -1823 +1837 @@
     type_ocean = type_ocean_omp
     version_ocean = version_ocean_omp
+    t_coupl = t_coupl_omp
 
     ok_veget=.true.

LMDZ5/trunk/libf/phylmd/cpl_mod.F90

@@ -132 +132 @@
      
     npas = itaufin/ iphysiq
-    nexca = 86400 / dtime
+!    nexca = 86400 / dtime
+    nexca = t_coupl / dtime
     WRITE(lunout,*)' ##### Ocean couple #####'
     WRITE(lunout,*)' Valeurs des pas de temps'

LMDZ5/trunk/libf/phylmd/surface_data.F90

@@ -26 +26 @@
   !$OMP THREADPRIVATE(version_ocean)
 
+  ! Pas de temps couplage atm/oce (en secondes)
+  REAL, SAVE             :: t_coupl
+  !$OMP THREADPRIVATE(t_coupl)
+
 END MODULE surface_data