Changeset 306

Timestamp:

Sep 24, 2011, 1:23:28 AM (13 years ago)

Author:

slebonnois

Message:

SLebonnois: correction de bugs dans la physique Titan:

• effg.F : Z doit etre en km, donc conversion
• optc*_1pt_2.F : On utilise cfffv11 et plus optfrac Du coup, les fichiers input testag* ne sont plus necessaires.
• phytrac.F : passage de la tendance aerosols en intensif dans tous les cas

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• cooling.F (modified) (6 diffs)
• effg.F (modified) (1 diff)
• ini_histday.h (modified) (2 diffs)
• ini_histmth.h (modified) (2 diffs)
• optci.F (modified) (3 diffs)
• optci_1pt_2.F (modified) (1 diff)
• optcv_1pt_2.F (modified) (1 diff)
• optfrac.old (moved) (moved from trunk/LMDZ.TITAN/libf/phytitan/optfrac.F)
• physiq.F (modified) (2 diffs)
• phytrac.F (modified) (2 diffs)
• tgmdat.F (modified) (1 diff)
• write_histday.h (modified) (2 diffs)
• write_histmth.h (modified) (2 diffs)

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

@@ -46 +46 @@
 c     spectral ranges. This common is initialized by radtitan.
 c   
+c     COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI
+c     This common is initialized by tgmdat.
+c
 c=======================================================================
 c-----------------------------------------------------------------------
@@ -81 +84 @@
       COMMON /IRTAUS/ dtaui,dtauip
 
+      COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI
+      REAL CSUBP,RSFI,RSFV,F0PI
+
 c   Local:
 c   ------
@@ -101 +107 @@
 c   -----
 
-      REAL RHOP,CSUBP,UBARI,RSFI
+      REAL RHOP,UBARI
       DATA RHOP/1.E4/      ! CONVERSION FROM PRESSURE TO MASS
       DATA UBARI/0.5/      ! MEAN COSINE FOR 2-STREAM
-      DATA RSFI/0.0/       ! SURFACE ALBEDO
       DATA WNOI/
      &    11.500,    20.000,    31.250,    50.000,    75.000,
@@ -129 +134 @@
 
 
-      save RHOP,UBARI,RSFI,WNOI,DWNI
+      save RHOP,UBARI,WNOI,DWNI
 
 c-----------------------------------------------------------------------
@@ -135 +140 @@
 c   Initialisations:
 c   ----------------
-
-      CSUBP = RCPD*1e4   ! HEAT CAPACITY OF TITANS ATMOSPHERE in CGS
 
       UBARI2=1./1.66
@@ -234 +237 @@
           ENDDO
 c          BSURF=PLNCK( WAVEN, TEMP(ig,NL), DW)
-           BSURF=.2*PLNCK
+           BSURF=.2*PLNCK*(1-RSFI)
         FUPI(ig,NL)=BSURF * 2.*RPI*UBARI + RSFI*FDI(ig,NL)
         FUPIS(ig,NL,K)=BSURF*2.*RPI*UBARI+RSFI*FDIS(ig,NL,K)

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

@@ -1 +1 @@
       FUNCTION EFFG(Z)
 #include "YOMCST.h"
-      EFFG = RG * (RA/(RA + Z ) )**2
+! RA en m, Z en km...
+      EFFG = RG * (RA/(RA + Z*1000. ) )**2
       RETURN
       END

trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h

@@ -113 +113 @@
      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
      .                    "ave(X)", zsto,zout)
+
+             if (clouds.eq.1) then
+
              CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
@@ -127 +130 @@
 c--------------
 c ----- SATURATION ESP NUAGES
-             if (clouds.eq.1) then
                CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,

trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h

@@ -117 +117 @@
      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
      .                    "ave(X)", zsto,zout)
+
+             if (clouds.eq.1) then
+
              CALL histdef(nid_mth, "qnoy","nb tot noy" , "n/m2",
      .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
@@ -131 +134 @@
 c--------------
 c ----- SATURATION ESP NUAGES
-             if (clouds.eq.1) then
                CALL histdef(nid_mth,"ch4sat", "saturation CH4", "--",
      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,

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

@@ -11 +11 @@
       REAL    ykim(klon,klev,nqtot)
       real    qaer(klon,klev,nqtot)
-      integer nmicro
+      integer nmicro,IPRINT
 c   ---------
 
@@ -93 +93 @@
 
       integer iopti,iwarning     ! iopti: premier appel, une seule boucle sur les l.d'o.
-      integer ig,seulmtunpt
+      integer ig,seulmtunpt,iout
       save iopti,iwarning,seulmtunpt
       data iopti,iwarning,seulmtunpt/0,0,0/
@@ -265 +265 @@
 c        endif
         
+        iout=0
+c       if ((microfi.eq.0).or.(ig.eq.(klon/2+16))) iout=1
         if (seulmtunpt.eq.0) then
           call optci_1pt2(zqaer_1pt,rmcbar(ig,:),xfbar(ig,:,:),
-     &                   iopti,IPRINT)
+     &                   iopti,iout)
            iopti = 1
         endif

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

@@ -192 +192 @@
              if (iopti.eq.0) then
 
-c       CALL OPTFRAC(XMONO,10000./WNOI(K)
-c     &                         ,QEXT,QSCT,QABS,QBAR)
-
-
        CALL CFFFV11(1.e-2/WNOI(K),REALI(K),XIMGI(K),RF(inq),2.
      &                ,XMONO,QSCT,QEXT,QABS,QBAR)
-
 
               QF1(inq,K)=QEXT*XRULE

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

@@ -148 +148 @@
 
             if(ioptv.eq.0.and.J.eq.1) then
-
-c        CALL OPTFRAC(XMONO,10000./WNOV(K)
-c     &                        ,QEXT,QSCT,QABS,QBAR)
 
         CALL CFFFV11(1.e-2/WNOV(K),REALV(K),XIMGV(K),RF(inq),2.

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

@@ -413 +413 @@
 
       COMMON /TAUD/   TAUHID,TAUCID,TAUGID,TAUHVD,TAUCVD,TAUGVD
+      COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI
+      REAL CSUBP,RSFI,RSFV,F0PI
+
 * common relatifs au nuages
       real rmcbar(ngrid,NLAYER),xfbar(ngrid,NLAYER,4) 
@@ -992 +995 @@
 c
       DO i = 1, klon
-            dlw(i) = - 4.0*RSIGMA*ftsol(i)**3 
+            dlw(i) = - 4.0*(1-RSFI)*RSIGMA*ftsol(i)**3 
       ENDDO
 

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

@@ -728 +728 @@
                d_tr_mph(i,l,iq) = (qaer(i,l,iq)-qaer0(i,l,iq))/
      &                            ptimestep
+c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
+               d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)
              ENDDO
            ENDDO
@@ -737 +739 @@
              DO i = 1, klon
                d_tr_mph(i,l,iq)=(qaer(i,l,iq)-qaer0(i,l,iq))/ptimestep
+c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
+               d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)
              ENDDO
            ENDDO
          ENDDO
-
-         do iq=1,nmicro
-          DO l=1,llm
-           DO i = 1, klon
-c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
-             d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)
-           ENDDO
-          ENDDO
-         enddo
 
        ENDIF   ! microfi

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

@@ -30 +30 @@
       DATA RHOP/1.E4/
 C RSF IS THE SURFACE REFLECTANCE FOR VIS AND IR
-      DATA RSFV,RSFI/0.0,0.0/
+      DATA RSFV,RSFI/0.15,0.05/
 C FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2
       DATA F0PI/1.5E4/

trunk/LMDZ.TITAN/libf/phytitan/write_histday.h

@@ -112 +112 @@
        CALL histwrite(nid_day,"qaer",itau_w,zx_tmp_3d,
      .                                   iim*jjmp1*klev,ndex3d)
+
+             if (clouds.eq.1) then
+
 c    -------   NB NOY TOT
                do i=1,klon
@@ -150 +153 @@
 c --------------
 c ----- SATURATION ESP NUAGES
-             if (clouds.eq.1) then
-
        CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satch4,zx_tmp_3d)
        CALL histwrite(nid_day,"ch4sat", itau_w, zx_tmp_3d,

trunk/LMDZ.TITAN/libf/phytitan/write_histmth.h

@@ -110 +110 @@
        CALL histwrite(nid_mth,"qaer",itau_w,zx_tmp_3d,
      .                                   iim*jjmp1*klev,ndex3d)
+
+             if (clouds.eq.1) then
+
 c    -------   NB NOY TOT
                do i=1,klon
@@ -148 +151 @@
 c --------------
 c ----- SATURATION ESP NUAGES
-             if (clouds.eq.1) then
-
        CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satch4,zx_tmp_3d)
        CALL histwrite(nid_mth,"ch4sat", itau_w, zx_tmp_3d,