Changeset 1565 for LMDZ5/trunk/libf/phylmd/radlwsw.F90

Timestamp:

Aug 31, 2011, 2:53:29 PM (13 years ago)

Author:

jghattas

Message:

Added interface with chemestry model REPROBUS :

• Compile LMDZ together with Reprobus code (dependecies in both directions) and cpp key REPROBUS :

./makelmdz_fcm -ext_src my_path_to_reprobus -cpp REPROBUS ...

• For running, add type_trac=repr in run.def.

/Marion Marchand, JG

File:

• LMDZ5/trunk/libf/phylmd/radlwsw.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/radlwsw.F90

@@ -30 +30 @@
 
   USE DIMPHY
-  use assert_m, only: assert
+  USE assert_m, ONLY : assert
+  USE infotrac, ONLY : type_trac
+#ifdef REPROBUS
+  USE CHEM_REP, ONLY : solaireTIME, ok_SUNTIME, ndimozon
+#endif
 
   !======================================================================
@@ -229 +233 @@
   !
   PSCT = solaire/zdist/zdist
+
+  IF (type_trac == 'repr') THEN
+#ifdef REPROBUS
+     if(ok_SUNTIME) PSCT = solaireTIME/zdist/zdist
+     print*,'Constante solaire: ',PSCT*zdist*zdist
+#endif
+  END IF
+
   DO j = 1, nb_gr
     iof = kdlon*(j-1)
@@ -281 +293 @@
       ENDDO
     ENDDO
+
+    IF (type_trac == 'repr') THEN
+#ifdef REPROBUS
+       ndimozon = size(wo, 3)
+       CALL RAD_INTERACTIF(POZON,iof)
+#endif
+    END IF
+
     !
     DO k = 1, kflev+1