Changeset 887 for trunk/LMDZ.TITAN/libf

Timestamp:

Feb 18, 2013, 2:46:14 PM (12 years ago)

Author:

slebonnois

Message:

SL: add rcm1d tool in Venus and Titan physics to run the code in 1D, and associated small modifications

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• effg.F (modified) (1 diff)
• physiq.F (modified) (2 diffs)
• profile.F (added)
• rcm1d.F (added)

trunk/LMDZ.TITAN/libf/phytitan/effg.F

@@ -6 +6 @@
 !    (et dans physique, attention a clmain et autres...)
 
-!      EFFG = RG * (RA/(RA + Z ) )**2
+      EFFG = RG * (RA/(RA + Z ) )**2
 
 ! Pour l'instant:
-      EFFG = RG 
+!     EFFG = RG 
       RETURN
       END

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -780 +780 @@
       DO l=1,klev
          DO i=1,klon
-            zzlay(i,l)=zphi(i,l)/RG
+c           zzlay(i,l)=zphi(i,l)/RG
 c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE:
-c           zzlay(i,l)=RG*RA*RA/(RG*RA-zphi(i,l))-RA
+            zzlay(i,l)=RG*RA*RA/(RG*RA-zphi(i,l))-RA
          ENDDO
       ENDDO
@@ -1548 +1548 @@
 #endif
 
-c ---------- Sorties testphys1d -------------
-      
-      if (klon.eq.1) then
-        call writeg1d(klon,klev,t_seri,"Temp","Temperature")
-        call writeg1d(klon,1,ftsol,"tsurf","Surface Temp")
-      DO k = 1, klev
-      DO i = 1, klon
-c        tmpout(i,k) = real(heat(i,k))
-         tmpout(i,k) = heat(i,k)
-      ENDDO
-      ENDDO
-        call writeg1d(klon,klev,tmpout,
-     .                 "heat","Solar heating")
-      DO k = 1, klev
-      DO i = 1, klon
-c        tmpout(i,k) = real(dtrad(i,k))
-         tmpout(i,k) = dtrad(i,k)
-      ENDDO
-      ENDDO
-        call writeg1d(klon,klev,tmpout,
-     .                 "dtrad","IR cooling")
-        call writeg1d(klon,klev,lwnet,"lwnet","Net LW flux")
-        call writeg1d(klon,klev,swnet,"swnet","Net SW flux")
-        call writeg1d(klon,klev,fluxt,"flux_vdf","Turbulent flux")
-        call writeg1d(klon,klev,flux_ajs,"flux_ajs","Dry adjust. flux")
-        call writeg1d(klon,klev,flux_ec,"flux_ec","Ec flux")
-c       call writeg1d(klon,1,solsw,"surfsw","Net SW flux at surface")
-c       call writeg1d(klon,1,sollw,"surflw","Net LW flux at surface")
-        call writeg1d(klon,1,radsol,"surfnet","Net flux at surface")
-        call writeg1d(klon,klev,d_t_vdf,"dt_vdf","DT from clmain")
-        call writeg1d(klon,klev,d_t_ajs,"dt_ajs","DT from ajsec")
-        call writeg1d(klon,klev,d_t_ec,"dt_ec","DT from Ec")
-      endif
-      
 c====================================================================
 c Si c'est la fin, il faut conserver l'etat de redemarrage