Ignore:
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:
12 edited
1 moved

Legend:

Unmodified
Added
Removed

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

    r175 r306  
    4646c     spectral ranges. This common is initialized by radtitan.
    4747c   
     48c     COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI
     49c     This common is initialized by tgmdat.
     50c
    4851c=======================================================================
    4952c-----------------------------------------------------------------------
    … …  
    8184      COMMON /IRTAUS/ dtaui,dtauip
    8285
     86      COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI
     87      REAL CSUBP,RSFI,RSFV,F0PI
     88
    8389c   Local:
    8490c   ------
    … …  
    101107c   -----
    102108
    103       REAL RHOP,CSUBP,UBARI,RSFI
     109      REAL RHOP,UBARI
    104110      DATA RHOP/1.E4/      ! CONVERSION FROM PRESSURE TO MASS
    105111      DATA UBARI/0.5/      ! MEAN COSINE FOR 2-STREAM
    106       DATA RSFI/0.0/       ! SURFACE ALBEDO
    107112      DATA WNOI/
    108113     &    11.500,    20.000,    31.250,    50.000,    75.000,
    … …  
    129134
    130135
    131       save RHOP,UBARI,RSFI,WNOI,DWNI
     136      save RHOP,UBARI,WNOI,DWNI
    132137
    133138c-----------------------------------------------------------------------
    … …  
    135140c   Initialisations:
    136141c   ----------------
    137 
    138       CSUBP = RCPD*1e4   ! HEAT CAPACITY OF TITANS ATMOSPHERE in CGS
    139142
    140143      UBARI2=1./1.66
    … …  
    234237          ENDDO
    235238c          BSURF=PLNCK( WAVEN, TEMP(ig,NL), DW)
    236            BSURF=.2*PLNCK
     239           BSURF=.2*PLNCK*(1-RSFI)
    237240        FUPI(ig,NL)=BSURF * 2.*RPI*UBARI + RSFI*FDI(ig,NL)
    238241        FUPIS(ig,NL,K)=BSURF*2.*RPI*UBARI+RSFI*FDIS(ig,NL,K)

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

    r175 r306  
    11      FUNCTION EFFG(Z)
    22#include "YOMCST.h"
    3       EFFG = RG * (RA/(RA + Z ) )**2
     3! RA en m, Z en km...
     4      EFFG = RG * (RA/(RA + Z*1000. ) )**2
    45      RETURN
    56      END

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

    r175 r306  
    113113     .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    114114     .                    "ave(X)", zsto,zout)
     115
     116             if (clouds.eq.1) then
     117
    115118             CALL histdef(nid_day, "qnoy","nb tot noy" , "n/m2",
    116119     .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    … …  
    127130c--------------
    128131c ----- SATURATION ESP NUAGES
    129              if (clouds.eq.1) then
    130132               CALL histdef(nid_day,"ch4sat", "saturation CH4", "--",
    131133     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,

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

    r175 r306  
    117117     .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    118118     .                    "ave(X)", zsto,zout)
     119
     120             if (clouds.eq.1) then
     121
    119122             CALL histdef(nid_mth, "qnoy","nb tot noy" , "n/m2",
    120123     .                    iim,jjmp1,nhori, klev,1,klev,nvert, 32,
    … …  
    131134c--------------
    132135c ----- SATURATION ESP NUAGES
    133              if (clouds.eq.1) then
    134136               CALL histdef(nid_mth,"ch4sat", "saturation CH4", "--",
    135137     .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,

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

    r175 r306  
    1111      REAL    ykim(klon,klev,nqtot)
    1212      real    qaer(klon,klev,nqtot)
    13       integer nmicro
     13      integer nmicro,IPRINT
    1414c   ---------
    1515
    … …  
    9393
    9494      integer iopti,iwarning     ! iopti: premier appel, une seule boucle sur les l.d'o.
    95       integer ig,seulmtunpt
     95      integer ig,seulmtunpt,iout
    9696      save iopti,iwarning,seulmtunpt
    9797      data iopti,iwarning,seulmtunpt/0,0,0/
    … …  
    265265c        endif
    266266       
     267        iout=0
     268c       if ((microfi.eq.0).or.(ig.eq.(klon/2+16))) iout=1
    267269        if (seulmtunpt.eq.0) then
    268270          call optci_1pt2(zqaer_1pt,rmcbar(ig,:),xfbar(ig,:,:),
    269      &                   iopti,IPRINT)
     271     &                   iopti,iout)
    270272           iopti = 1
    271273        endif

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

    r175 r306  
    192192             if (iopti.eq.0) then
    193193
    194 c       CALL OPTFRAC(XMONO,10000./WNOI(K)
    195 c     &                         ,QEXT,QSCT,QABS,QBAR)
    196 
    197 
    198194       CALL CFFFV11(1.e-2/WNOI(K),REALI(K),XIMGI(K),RF(inq),2.
    199195     &                ,XMONO,QSCT,QEXT,QABS,QBAR)
    200 
    201196
    202197              QF1(inq,K)=QEXT*XRULE

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

    r175 r306  
    148148
    149149            if(ioptv.eq.0.and.J.eq.1) then
    150 
    151 c        CALL OPTFRAC(XMONO,10000./WNOV(K)
    152 c     &                        ,QEXT,QSCT,QABS,QBAR)
    153150
    154151        CALL CFFFV11(1.e-2/WNOV(K),REALV(K),XIMGV(K),RF(inq),2.

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

    r175 r306  
    413413
    414414      COMMON /TAUD/   TAUHID,TAUCID,TAUGID,TAUHVD,TAUCVD,TAUGVD
     415      COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI
     416      REAL CSUBP,RSFI,RSFV,F0PI
     417
    415418* common relatifs au nuages
    416419      real rmcbar(ngrid,NLAYER),xfbar(ngrid,NLAYER,4)
    … …  
    992995c
    993996      DO i = 1, klon
    994             dlw(i) = - 4.0*RSIGMA*ftsol(i)**3
     997            dlw(i) = - 4.0*(1-RSFI)*RSIGMA*ftsol(i)**3
    995998      ENDDO
    996999

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

    r176 r306  
    728728               d_tr_mph(i,l,iq) = (qaer(i,l,iq)-qaer0(i,l,iq))/
    729729     &                            ptimestep
     730c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
     731               d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)
    730732             ENDDO
    731733           ENDDO
    … …  
    737739             DO i = 1, klon
    738740               d_tr_mph(i,l,iq)=(qaer(i,l,iq)-qaer0(i,l,iq))/ptimestep
     741c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
     742               d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)
    739743             ENDDO
    740744           ENDDO
    741745         ENDDO
    742 
    743          do iq=1,nmicro
    744           DO l=1,llm
    745            DO i = 1, klon
    746 c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
    747              d_tr_mph(i,l,iq) = d_tr_mph(i,l,iq)*RG/delp(i,l)
    748            ENDDO
    749           ENDDO
    750          enddo
    751746
    752747       ENDIF   ! microfi

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

    r3 r306  
    3030      DATA RHOP/1.E4/
    3131C RSF IS THE SURFACE REFLECTANCE FOR VIS AND IR
    32       DATA RSFV,RSFI/0.0,0.0/
     32      DATA RSFV,RSFI/0.15,0.05/
    3333C FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2
    3434      DATA F0PI/1.5E4/

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

    r175 r306  
    112112       CALL histwrite(nid_day,"qaer",itau_w,zx_tmp_3d,
    113113     .                                   iim*jjmp1*klev,ndex3d)
     114
     115             if (clouds.eq.1) then
     116
    114117c    -------   NB NOY TOT
    115118               do i=1,klon
    … …  
    150153c --------------
    151154c ----- SATURATION ESP NUAGES
    152              if (clouds.eq.1) then
    153 
    154155       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satch4,zx_tmp_3d)
    155156       CALL histwrite(nid_day,"ch4sat", itau_w, zx_tmp_3d,

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

    r175 r306  
    110110       CALL histwrite(nid_mth,"qaer",itau_w,zx_tmp_3d,
    111111     .                                   iim*jjmp1*klev,ndex3d)
     112
     113             if (clouds.eq.1) then
     114
    112115c    -------   NB NOY TOT
    113116               do i=1,klon
    … …  
    148151c --------------
    149152c ----- SATURATION ESP NUAGES
    150              if (clouds.eq.1) then
    151 
    152153       CALL gr_fi_ecrit(klev,klon,iim,jjmp1, satch4,zx_tmp_3d)
    153154       CALL histwrite(nid_mth,"ch4sat", itau_w, zx_tmp_3d,
Note: See TracChangeset for help on using the changeset viewer.