Changeset 1424 for LMDZ4

Timestamp:

Sep 1, 2010, 11:20:28 AM (14 years ago)

Author:

jghattas

Message:

• Modified name of tracers in history files to their real name(from tracer.def) ; for exemple RN instead of trac01
• Corrected allocation of d_ts.

Location:

LMDZ4/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (1 diff)
• phys_output_mod.F90 (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -195 +195 @@
       allocate(d_t_lif(klon,klev),d_t_ec(klon,klev))
       allocate(d_u_lif(klon,klev),d_v_lif(klon,klev))
-      allocate(d_ts(klon,klev), d_tr(klon,klev,nbtr))
+      allocate(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
       allocate(topswad_aero(klon), solswad_aero(klon))
       allocate(topswai_aero(klon), solswai_aero(klon))

LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

@@ -1223 +1223 @@
      ENDIF
 
-!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
-    ENDDO
-!
-    if (nqtot>=3) THEN
+
+    IF (nqtot>=3) THEN
      DO iq=3,nqtot  
        iiq=niadv(iq)
-  CALL histdef3d (iff, o_trac(iq-2)%flag,o_trac(iq-2)%name,ttext(iiq), "-" )
+       o_trac(iq-2) = ctrl_out((/ 4, 5, 1, 1, 1 /),tname(iiq))
+       CALL histdef3d (iff, o_trac(iq-2)%flag,o_trac(iq-2)%name,'Tracer '//ttext(iiq), "-" )
      ENDDO
-    endif
+    ENDIF
 
         CALL histend(nid_files(iff))