Changeset 1376 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_output_mod.F90

Timestamp:

May 5, 2010, 3:23:18 PM (14 years ago)

Author:

musat

Message:

Output all tracers defined in .def in hist files with dynamic
declaration of LMDZ atmospheric tracers' output levels
Add 6 pseudo-water tracers with and without transport by boundary layer
IM

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_output_mod.F90 (modified) (5 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_output_mod.F90

@@ -11 +11 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-MODULE phys_output_mod
+MODULE phys_output_mod 
 
   IMPLICIT NONE
@@ -525 +525 @@
   type(ctrl_out),save :: o_dvlif        = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvlif')
 
-! Attention a refaire correctement
-  type(ctrl_out),save,dimension(2) :: o_trac         = (/ ctrl_out((/ 4, 10, 10, 10, 10 /),'trac01'), &
-                                                     ctrl_out((/ 4, 10, 10, 10, 10 /),'trac02') /)
     CONTAINS
 
@@ -539 +536 @@
        type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
        ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
-       new_aod, aerosol_couple)   
+       new_aod, aerosol_couple, o_trac)   
 
 
@@ -602 +599 @@
   real, dimension(nfiles), save     :: phys_out_latmin        = (/ -90., -90., -90., -90., -90. /)
   real, dimension(nfiles), save     :: phys_out_latmax        = (/ 90., 90., 90., 90., 90. /)
-  
-  
-
-! 
+
+!IM definition dynamique flag o_trac pour sortie traceurs
+  INTEGER :: nq
+  type(ctrl_out) :: o_trac(nqtot)
+  CHARACTER(len=8) :: solsym(nqtot)
+
    print*,'Debut phys_output_mod.F90'
 ! Initialisations (Valeurs par defaut
@@ -1220 +1219 @@
      ENDIF
 
-      if (nqtot>=3) THEN
-!Attention    DO iq=3,nqtot
-    DO iq=3,4  
+!IM traceurs dynamiques
+    DO nq=1,nqtot
+      IF(nq.LT.10) THEN
+       WRITE(solsym(nq),'(i1)') nq
+       o_trac(nq)    =  ctrl_out((/ 4, 5, 1, 1, 1 /),'trac0'//TRIM(solsym(nq)))
+      ELSE
+       WRITE(solsym(nq),'(i2)') nq
+       o_trac(nq)    =  ctrl_out((/ 4, 5, 1, 1, 1 /),'trac'//TRIM(solsym(nq)))
+      ENDIF
+     WRITE(*,*) 'nq, o_trac(nq)=',nq, o_trac(nq)
+    ENDDO
+!
+    if (nqtot>=3) THEN
+     DO iq=3,nqtot  
        iiq=niadv(iq)
-! CALL histdef3d (iff, o_trac%flag,'o_'//tnom(iq)%name,ttext(iiq), "-" )
   CALL histdef3d (iff, o_trac(iq-2)%flag,o_trac(iq-2)%name,ttext(iiq), "-" )
-    ENDDO
-      endif
+     ENDDO
+    endif
 
         CALL histend(nid_files(iff))