Ignore:
Timestamp:
Mar 6, 2002, 3:58:31 PM (23 years ago)
Author:
lmdz
Message:

Inclusion des modifs de D. Hauglustaine pour la version 1 de INCA
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/phylmd/physiq.F

    r303 r344  
    11c $Header$
    22c
    3       SUBROUTINE physiq (nlon,nlev,nqmax  ,
    4      .            debut,lafin,rjourvrai,rjour_ecri,gmtime,pdtphys,
    5      .            paprs,pplay,pphi,pphis,paire,presnivs,clesphy0,
    6      .            u,v,t,qx,
    7      .            omega,
    8      .            d_u, d_v, d_t, d_qx, d_ps)
     3      SUBROUTINE physiq (nlon,
     4     .                   nlev,
     5     .                   nqmax,
     6     .                   debut,
     7     .                   lafin,
     8     .                   rjourvrai,
     9     .                   rjour_ecri,
     10     .                   gmtime,
     11     .                   pdtphys,
     12     .                   paprs,
     13     .                   pplay,
     14     .                   pphi,
     15     .                   pphis,
     16     .                   paire,
     17     .                   presnivs,
     18     .                   clesphy0,
     19     .                   u,
     20     .                   v,
     21     .                   t,
     22     .                   qx,
     23     .                   omega,
     24#ifdef INCA_CH4
     25     .                   flxmass_w,
     26#endif
     27     .                   d_u,
     28     .                   d_v,
     29     .                   d_t,
     30     .                   d_qx,
     31     .                   d_ps)
    932      USE ioipsl
     33#ifdef INCA
     34      USE chemshut
     35      USE obs_pos
     36#endif
    1037      IMPLICIT none
    1138c======================================================================
    … …  
    7198      INTEGER npas, nexca, itimestep
    7299      logical rnpb
     100#ifdef INCA
     101      parameter(rnpb=.false.)
     102#else
    73103      parameter(rnpb=.true.)
     104#endif
    74105      PARAMETER (npas=1440)
    75106      PARAMETER (nexca=48)
    … …  
    174205      REAL omega(klon,klev)
    175206
     207      REAL flxmass_w(klon,klev)
     208
    176209      REAL d_u(klon,klev)
    177210      REAL d_v(klon,klev)
    … …  
    326359      REAL frac_impa(klon,klev) ! fractions d'aerosols lessivees (impaction)
    327360      REAL frac_nucl(klon,klev) ! idem (nucleation)
     361
     362#ifdef INCA
     363      REAL          :: calday
     364#endif
     365
    328366cAA
    329367      REAL rain_fall(klon) ! pluie
    … …  
    403441c Variables locales
    404442c
     443      REAL rhcl(klon,klev)    ! humiditi relative ciel clair
    405444      REAL dialiq(klon,klev)  ! eau liquide nuageuse
    406445      REAL diafra(klon,klev)  ! fraction nuageuse
    … …  
    14011440         ENDIF
    14021441c
     1442
     1443#ifdef INCA
     1444           calday = FLOAT(MOD(NINT(xjour),360)) + gmtime
     1445           print *, 'initial time ', xjour, calday
     1446#ifdef INCAINFO
     1447           PRINT *, 'Appel CHEMINI ...'
     1448#endif
     1449           CALL chemini( rpi,
     1450     $                   rg,
     1451     $                   ra,
     1452     $                   paire,
     1453     $                   rlat,
     1454     $                   rlon,
     1455     $                   calday,
     1456     $                   tracnam,
     1457     $                   natsnam,
     1458     $                   mxoutflds,
     1459     $                   outinst,
     1460     $                   outtimav,
     1461     $                   klon,
     1462     $                   nqmax)
     1463#ifdef INCAINFO
     1464           PRINT *, 'OK.'
     1465#endif
     1466#endif
     1467
    14031468      ENDIF
    14041469c
    … …  
    18481913     .           pfrac_impa, pfrac_nucl, pfrac_1nucl,
    18491914     .           frac_impa, frac_nucl,
    1850      .           prfl, psfl)
     1915     .           prfl, psfl, rhcl)
    18511916      DO k = 1, klev
    18521917      DO i = 1, klon
    … …  
    19291994      ENDDO
    19301995      ENDDO
     1996
     1997#ifdef INCA
     1998           calday = FLOAT(julien) + gmtime
     1999
     2000#ifdef INCA_AER
     2001      call AEROSOL_METEO_CALC(calday,pdtphys,pplay,paprs,t,pmflxr,pmflxs,
     2002     &   prfl,psfl,pctsrf,paire)
     2003#endif
     2004
     2005#ifdef INCAINFO
     2006           PRINT *, 'Appel CHEMHOOK_BEGIN ...'
     2007#endif
     2008           CALL chemhook_begin (calday,
     2009     $                          pctsrf(1,3),
     2010     $                          rlat,
     2011     $                          rlon,
     2012     $                          paire,
     2013     $                          paprs,
     2014     $                          pplay,
     2015     $                          ycoefh,
     2016     $                          pphi,
     2017     $                          t_seri,
     2018     $                          u,
     2019     $                          v,
     2020     $                          wo,
     2021     $                          q_seri,
     2022     $                          zxtsol,
     2023     $                          zxsnow,
     2024     $                          solsw,
     2025     $                          albsol,
     2026     $                          rain_fall,
     2027     $                          snow_fall,
     2028     $                          itop_con,
     2029     $                          ibas_con,
     2030     $                          cldfra,
     2031     $                          iim,
     2032     $                          jjm,
     2033     $                          tr_seri)
     2034#ifdef INCAINFO
     2035           PRINT *, 'OK.'
     2036#endif
     2037#endif
     2038
    19312039c
    19322040c Calculer les parametres optiques des nuages et quelques
    … …  
    21242232C
    21252233      call phytrac (rnpb,
     2234     I                   itap, julien, gmtime,
    21262235     I                   debut,lafin,
    21272236     I                   nqmax-2,
    21282237     I                   nlon,nlev,dtime,
    2129      I                   t,paprs,pplay,
     2238     I                   u,v,t,paprs,pplay,
    21302239     I                   pmfu, pmfd, pen_u, pde_u, pen_d, pde_d,
    21312240     I                   ycoefh,yu1,yv1,ftsol,pctsrf,rlat,
    21322241     I                   frac_impa, frac_nucl,
    2133      I                   rlon,presnivs,paire,pphis,
     2242     I                   rlon,paire,pphis,pphi,
     2243     I                   albsol,
     2244     I                   qx(1,1,1), rhcl,
     2245     I                   cldfra, rneb, diafra, cldliq, itop_con,
     2246     I                   ibas_con,
     2247     I                   pmflxr,pmflxs,prfl,psfl,flxmass_w,
    21342248     O                   tr_seri)
    21352249      ENDIF
    … …  
    28222936cc         CALL ecrirega(84,d_q_lsc)
    28232937      ENDIF
     2938
     2939#ifdef INCA
     2940#ifdef INCAINFO
     2941           PRINT *, 'Appel CHEMHOOK_END ...'
     2942#endif
     2943           CALL chemhook_end (calday,
     2944     $                        dtime,
     2945     $                        pplay,
     2946     $                        t_seri,
     2947     $                        tr_seri,
     2948     $                        nbtr,
     2949     $                        paprs,
     2950     $                        anne_ini,
     2951     $                        day_ini,
     2952     $                        xjour)
     2953#ifdef INCAINFO
     2954           PRINT *, 'OK.'
     2955#endif
     2956#endif
     2957
    28242958c
    28252959c Convertir les incrementations en tendances
Note: See TracChangeset for help on using the changeset viewer.