Changeset 2580


Ignore:
Timestamp:
Jun 29, 2016, 2:23:19 PM (8 years ago)
Author:
idelkadi
Message:

Implementation du simulateur AIRS:
Le but du simulateur est de permettre la comparaison de proprietes macro-et microphysiques des nuages hauts de LMDZ avec celles restituees par les observations du satellite AIRS (Atmospheric IR Sounder). La methode est decrite dans Hendricks et al. Meteorol. Z., 2010

Location:
LMDZ5/trunk
Files:
1 added
8 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/DefLists/field_def_lmdz.xml

    r2563 r2580  
    837837    </field_group>
    838838
     839    <field_group id="fields_AIRS" domain_ref="dom_glo" operation="instant" freq_op="3h">
     840      <field id="alt_tropo" long_name="Tropopause pressure" unit="hPa" />
     841      <field id="map_prop_hc" long_name="Proportion of high clouds" unit=" "/>
     842      <field id="map_prop_hist" long_name="Proportion of high ice semi-transp clouds" unit=" "/>
     843      <field id="map_emis_hc" long_name="Emissivity of high clouds" unit=" "/>
     844      <field id="map_iwp_hc" long_name="Ice water path of high clouds" unit="g/m2"/>
     845      <field id="map_deltaz_hc" long_name="geom thickness of high clouds" unit="m"/>
     846      <field id="map_pcld_hc" long_name="cloud pressure of high clouds" unit="hPa"/>
     847      <field id="map_tcld_hc" long_name="cloud temperature of high clouds" unit="K"/>
     848      <field id="map_emis_hist" long_name="Emissivity of high ice st clouds" unit=" "/>
     849      <field id="map_iwp_hist" long_name="Ice water path of high ice st clouds" unit="g/m2"/>
     850      <field id="map_deltaz_hist" long_name="geom thickness of high ice st clouds" unit="m"/>
     851      <field id="map_rad_hist" long_name="ice crystals radius in high ice st clouds" unit="m"/>
     852      <field id="map_emis_Cb" long_name="Emissivity of high Cb clouds" unit=" "/>
     853      <field id="map_pcld_Cb" long_name="cloud pressure of high Cb clouds" unit="hPa"/>
     854      <field id="map_tcld_Cb" long_name="cloud temperature of high Cb clouds" unit="K"/>
     855      <field id="map_emis_Anv" long_name="Emissivity of high Anv clouds" unit=" "/>
     856      <field id="map_pcld_Anv" long_name="cloud pressure of high Anv clouds" unit="hPa"/>
     857      <field id="map_tcld_Anv" long_name="cloud temperature of high Anv clouds" unit="K"/>
     858      <field id="map_emis_ThCi" long_name="Emissivity of high ThCi clouds" unit=" "/>
     859      <field id="map_pcld_ThCi" long_name="cloud pressure of high ThCi clouds" unit="hPa"/>
     860      <field id="map_tcld_ThCi" long_name="cloud temperature of high ThCi clouds" unit="K"/>
     861      <field id="map_ntot" long_name="total AIRS cloud fraction" unit=" "/>
     862      <field id="map_hc" long_name="high clouds AIRS cloud fraction" unit=" "/>
     863      <field id="map_hist" long_name="high clouds ice st AIRS cloud fraction" unit=" "/>
     864      <field id="map_Cb" long_name="high clouds Cb AIRS cloud fraction" unit=" "/>
     865      <field id="map_ThCi" long_name="high clouds ThCi AIRS cloud fraction" unit=" "/>
     866      <field id="map_Anv" long_name="high clouds Anv AIRS cloud fraction" unit=" "/>
     867    </field_group>
     868
    839869</field_definition>
    840870

  • LMDZ5/trunk/DefLists/file_def_histhf_lmdz.xml

    r2563 r2580  
    411411            </field_group>
    412412
     413            <field_group operation="instant" freq_op="1ts" >
     414                 <field field_ref="alt_tropo" level="10" />
     415                 <field field_ref="map_prop_hc" level="10" />
     416                 <field field_ref="map_prop_hist" level="10" />
     417                 <field field_ref="map_emis_hc" level="10" />
     418                 <field field_ref="map_iwp_hc" level="10" />
     419                 <field field_ref="map_deltaz_hc" level="10" />
     420                 <field field_ref="map_pcld_hc" level="10" />
     421                 <field field_ref="map_tcld_hchc" level="10" />
     422                 <field field_ref="map_emis_hist" level="10" />
     423                 <field field_ref="map_iwp_hist" level="10" />
     424                 <field field_ref="map_deltaz_hist" level="10" />
     425                 <field field_ref="map_rad_hist" level="10" />
     426                 <field field_ref="map_emis_Cb" level="10" />
     427                 <field field_ref="map_pcld_Cb" level="10" />
     428                 <field field_ref="map_tcld_Cb" level="10" />
     429                 <field field_ref="map_emis_Anv" level="10" />
     430                 <field field_ref="map_pcld_Anv" level="10" />
     431                 <field field_ref="map_tcld_Anv" level="10" />
     432                 <field field_ref="map_emis_ThCi" level="10" />
     433                 <field field_ref="map_pcld_ThCi" level="10" />
     434                 <field field_ref="map_tcld_ThCi" level="10" />
     435                 <field field_ref="map_ntot" level="10" />
     436                 <field field_ref="map_hc" level="10" />
     437                 <field field_ref="map_hist" level="10" />
     438                 <field field_ref="map_Cb" level="10" />
     439                 <field field_ref="map_ThCi" level="10" />
     440                 <field field_ref="map_Anv" level="10" />
     441            </field_group>
     442
    413443            <!-- VARS 3D -->
    414444            <field_group operation="average" freq_op="1ts" axis_ref="presnivs">

  • LMDZ5/trunk/libf/phylmd/clesphys.h

    r2561 r2580  
    7878       REAL ecrit_LES
    7979       REAL freq_ISCCP, ecrit_ISCCP
    80        REAL freq_COSP
     80       REAL freq_COSP, freq_AIRS
    8181       LOGICAL :: ok_cosp,ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP
     82       LOGICAL :: ok_airs
    8283       INTEGER :: ip_ebil_phy, iflag_rrtm, iflag_ice_thermo, NSW, iflag_albedo
    8384       LOGICAL :: ok_chlorophyll
    … …  
    102103     &     , pasphys            , freq_outNMC, freq_calNMC              &
    103104     &     , lonmin_ins, lonmax_ins, latmin_ins, latmax_ins             &
    104      &     , freq_ISCCP, ecrit_ISCCP, freq_COSP                         &
     105     &     , freq_ISCCP, ecrit_ISCCP, freq_COSP, freq_AIRS              &
    105106     &     , cvl_corr                                                   &
    106107     &     , qsol0,albsno0,evap0                                        &
    … …  
    121122     &     , lev_histins, lev_histLES, lev_histdayNMC, levout_histNMC   &
    122123     &     , ok_histNMC                                                 &
    123      &     , type_run, ok_regdyn, ok_cosp                               &
     124     &     , type_run, ok_regdyn, ok_cosp, ok_airs                      &
    124125     &     , ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP                     &
    125126     &     , ip_ebil_phy                                                &

  • LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

    r2561 r2580  
    9898    REAL,SAVE           :: bl95_b0_omp, bl95_b1_omp
    9999    REAL,SAVE           :: freq_ISCCP_omp, ecrit_ISCCP_omp
    100     REAL,SAVE           :: freq_COSP_omp
     100    REAL,SAVE           :: freq_COSP_omp, freq_AIRS_omp
    101101    real,SAVE           :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
    102102    real,SAVE           :: tau_cld_cv_omp, coefw_cld_cv_omp
    … …  
    190190    REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3)
    191191    CHARACTER*4, SAVE :: type_run_omp
    192     LOGICAL,SAVE :: ok_cosp_omp
     192    LOGICAL,SAVE :: ok_cosp_omp, ok_airs_omp
    193193    LOGICAL,SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
    194194    REAL,SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
    … …  
    462462    call getin('freq_COSP', freq_COSP_omp)
    463463
     464    !Config Key  = freq_AIRS
     465    !Config Desc = Frequence d'appel du simulateur AIRS en secondes;
     466    !              par defaut 10800, i.e. 3 heures
     467    !Config Def  = 10800.
     468    !Config Help = Used in ini_histdayAIRS.h
     469    !
     470    freq_AIRS_omp = 10800.
     471    call getin('freq_AIRS', freq_AIRS_omp)
     472
    464473    !
    465474    !Config Key  = ip_ebil_phy
    … …  
    16091618    ok_cosp_omp = .false.
    16101619    call getin('ok_cosp',ok_cosp_omp)
     1620
     1621    !
     1622    !Config Key  = ok_airs
     1623    !Config Desc =
     1624    !Config Def  = .false.
     1625    !Config Help =
     1626    !
     1627    ok_airs_omp = .false.
     1628    call getin('ok_airs',ok_airs_omp)
    16111629
    16121630    !
    … …  
    21112129    ecrit_ISCCP = ecrit_ISCCP_omp
    21122130    freq_COSP = freq_COSP_omp
     2131    freq_AIRS = freq_AIRS_omp
    21132132    ok_ade = ok_ade_omp
    21142133    ok_aie = ok_aie_omp
    … …  
    21612180    type_run = type_run_omp
    21622181    ok_cosp = ok_cosp_omp
     2182    ok_airs = ok_airs_omp
     2183
    21632184    ok_mensuelCOSP = ok_mensuelCOSP_omp
    21642185    ok_journeCOSP = ok_journeCOSP_omp
    … …  
    23082329    write(lunout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
    23092330    write(lunout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
     2331    write(lunout,*)' Frequence appel simulateur AIRS, freq_AIRS =', freq_AIRS
    23102332    write(lunout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
    23112333    write(lunout,*)' Excentricite = ',R_ecc
    … …  
    24232445    write(lunout,*)' type_run = ',type_run
    24242446    write(lunout,*)' ok_cosp = ',ok_cosp
     2447    write(lunout,*)' ok_airs = ',ok_airs
     2448
    24252449    write(lunout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP
    24262450    write(lunout,*)' ok_journeCOSP = ',ok_journeCOSP

  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r2562 r2580  
    4141
    4242!!! 2D
     43
     44! Marine
     45
     46  TYPE(ctrl_out), SAVE :: o_alt_tropo = ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     47  'alt_tropo','Tropopause pressure','hPa',&
     48   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     49    "inst(X)", "inst(X)","inst(X)" /))
     50
     51  TYPE(ctrl_out), SAVE :: o_map_prop_hc = ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     52  'map_prop_hc','Proportion of high clouds',' ',&
     53   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     54    "inst(X)", "inst(X)","inst(X)" /))
     55
     56  TYPE(ctrl_out), SAVE :: o_map_prop_hist = &
     57  ctrl_out((/1,1,1,1,1,1,10,10,10/),&
     58  'map_prop_hist','Proportion of high ice semi-transp clouds',' ',&
     59   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     60    "inst(X)", "inst(X)","inst(X)" /))
     61
     62  TYPE(ctrl_out), SAVE :: o_map_emis_hc = &
     63  ctrl_out((/1,1,1,1,1,1,10,10,10/),&
     64  'map_emis_hc','Emissivity of high clouds',' ',&
     65   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     66    "inst(X)", "inst(X)","inst(X)" /))
     67
     68  TYPE(ctrl_out), SAVE :: o_map_iwp_hc = &
     69  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     70  'map_iwp_hc','Ice water path of high clouds','g/m2',&
     71   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     72    "inst(X)", "inst(X)","inst(X)" /))
     73
     74  TYPE(ctrl_out), SAVE :: o_map_deltaz_hc = &
     75  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     76  'map_deltaz_hc','geom thickness of high clouds','m',&
     77   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     78    "inst(X)", "inst(X)","inst(X)" /))
     79
     80  TYPE(ctrl_out), SAVE :: o_map_pcld_hc = &
     81  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     82  'map_pcld_hc','cloud pressure of high clouds','hPa',&
     83   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     84    "inst(X)", "inst(X)","inst(X)" /))
     85
     86   TYPE(ctrl_out), SAVE :: o_map_tcld_hc = &
     87  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     88  'map_tcld_hc','cloud temperature of high clouds','K',&
     89   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     90    "inst(X)", "inst(X)","inst(X)" /))
     91
     92
     93  TYPE(ctrl_out), SAVE :: o_map_emis_hist = &
     94  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     95  'map_emis_hist','Emissivity of high ice st clouds',' ',&
     96   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     97    "inst(X)", "inst(X)","inst(X)" /))
     98
     99  TYPE(ctrl_out), SAVE :: o_map_iwp_hist = &
     100  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     101  'map_iwp_hist','Ice water path of high ice st clouds','g/m2',&
     102   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     103    "inst(X)", "inst(X)","inst(X)" /))
     104
     105  TYPE(ctrl_out), SAVE :: o_map_deltaz_hist = &
     106  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     107  'map_deltaz_hist','geom thickness of high ice st clouds','m',&
     108   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     109    "inst(X)", "inst(X)","inst(X)" /))
     110
     111  TYPE(ctrl_out), SAVE :: o_map_rad_hist = &
     112  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     113  'map_rad_hist','ice crystals radius in high ice st clouds','µm',&
     114   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     115    "inst(X)", "inst(X)","inst(X)" /))
     116
     117
     118 TYPE(ctrl_out), SAVE :: o_map_emis_Cb = &
     119  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     120  'map_emis_Cb','Emissivity of high Cb clouds',' ',&
     121   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     122    "inst(X)", "inst(X)","inst(X)" /))
     123
     124 TYPE(ctrl_out), SAVE :: o_map_pcld_Cb = &
     125  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     126  'map_pcld_Cb','cloud pressure of high Cb clouds','hPa',&
     127   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     128    "inst(X)", "inst(X)","inst(X)" /))
     129
     130 TYPE(ctrl_out), SAVE :: o_map_tcld_Cb = &
     131  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     132  'map_tcld_Cb','cloud temperature of high Cb clouds','K',&
     133   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     134    "inst(X)", "inst(X)","inst(X)" /))
     135
     136
     137 TYPE(ctrl_out), SAVE :: o_map_emis_Anv = &
     138  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     139  'map_emis_Anv','Emissivity of high Anv clouds',' ',&
     140   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     141    "inst(X)", "inst(X)","inst(X)" /))
     142
     143 TYPE(ctrl_out), SAVE :: o_map_pcld_Anv = &
     144  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     145  'map_pcld_Anv','cloud pressure of high Anv clouds','hPa',&
     146   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     147    "inst(X)", "inst(X)","inst(X)" /))
     148
     149  TYPE(ctrl_out), SAVE :: o_map_tcld_Anv = &
     150  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     151  'map_tcld_Anv','cloud temperature of high Anv clouds','K',&
     152   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     153    "inst(X)", "inst(X)","inst(X)" /))
     154
     155  TYPE(ctrl_out), SAVE :: o_map_emis_ThCi = &
     156  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     157  'map_emis_ThCi','Emissivity of high ThCi clouds',' ',&
     158   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     159    "inst(X)", "inst(X)","inst(X)" /))
     160
     161  TYPE(ctrl_out), SAVE :: o_map_pcld_ThCi = &
     162  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     163  'map_pcld_ThCi','cloud pressure of high ThCi clouds','hPa',&
     164   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     165    "inst(X)", "inst(X)","inst(X)" /))
     166
     167  TYPE(ctrl_out), SAVE :: o_map_tcld_ThCi = &
     168  ctrl_out((/10,10,1,10,10,10,10,10,10/),&
     169  'map_tcld_ThCi','cloud temperature of high ThCi clouds','K',&
     170   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     171    "inst(X)", "inst(X)","inst(X)" /))
     172
     173   TYPE(ctrl_out), SAVE :: o_map_ntot = &
     174  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     175  'map_ntot','total AIRS cloud fraction',' ',&
     176   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     177    "inst(X)", "inst(X)","inst(X)" /))
     178
     179  TYPE(ctrl_out), SAVE :: o_map_hc = &
     180  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     181  'map_hc','high clouds AIRS cloud fraction',' ',&
     182   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     183    "inst(X)", "inst(X)","inst(X)" /))
     184
     185  TYPE(ctrl_out), SAVE :: o_map_hist = &
     186  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     187  'map_hist','high clouds ice st AIRS cloud fraction',' ',&
     188   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     189    "inst(X)", "inst(X)","inst(X)" /))
     190
     191  TYPE(ctrl_out), SAVE :: o_map_Cb = &
     192  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     193  'map_Cb','high clouds Cb AIRS cloud fraction',' ',&
     194   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     195    "inst(X)", "inst(X)","inst(X)" /))
     196
     197 TYPE(ctrl_out), SAVE :: o_map_ThCi = &
     198  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     199  'map_ThCi','high clouds ThCi AIRS cloud fraction',' ',&
     200   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     201    "inst(X)", "inst(X)","inst(X)" /))
     202
     203 TYPE(ctrl_out), SAVE :: o_map_Anv = &
     204  ctrl_out((/1,1,1,1,1,10,10,10,10/),&
     205  'map_Anv','high clouds Anv AIRS cloud fraction',' ',&
     206   (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)",&
     207    "inst(X)", "inst(X)","inst(X)" /))
     208
     209
     210! Fin Marine
     211
    43212  TYPE(ctrl_out), SAVE :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10, 11, 11, 11 /), &
    44213    'flat', 'Latent heat flux', 'W/m2', (/ ('', i=1, 9) /))

  • LMDZ5/trunk/libf/phylmd/phys_output_var_mod.F90

    r2538 r2580  
    2424  REAL, ALLOCATABLE :: bils_latent(:) ! bilan de chaleur au sol
    2525  !$OMP THREADPRIVATE(bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
     26
     27! Marine
     28! Variables de sortie du simulateur AIRS
     29
     30  REAL, ALLOCATABLE :: map_prop_hc(:),map_prop_hist(:)
     31  REAL, ALLOCATABLE :: alt_tropo(:)
     32  REAL, ALLOCATABLE :: &
     33     & map_emis_hc(:),map_iwp_hc(:),map_deltaz_hc(:),&
     34     & map_pcld_hc(:),map_tcld_hc(:)
     35  REAL, ALLOCATABLE :: map_emis_hist(:),map_iwp_hist(:),map_deltaz_hist(:),&
     36     & map_rad_hist(:)
     37  REAL, ALLOCATABLE :: map_ntot(:),map_hc(:),map_hist(:)
     38  REAL, ALLOCATABLE :: map_Cb(:),map_ThCi(:),map_Anv(:)
     39  REAL, ALLOCATABLE :: map_emis_Cb(:),map_pcld_Cb(:),map_tcld_Cb(:)
     40  REAL, ALLOCATABLE :: map_emis_ThCi(:),map_pcld_ThCi(:),map_tcld_ThCi(:)
     41  REAL, ALLOCATABLE :: map_emis_Anv(:),map_pcld_Anv(:),map_tcld_Anv(:)
    2642
    2743
    … …  
    103119    allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
    104120
     121! Marine
     122! Variables de sortie simulateur AIRS
     123
     124      allocate (map_prop_hc(klon),map_prop_hist(klon))
     125      allocate (alt_tropo(klon))
     126      allocate &
     127     & (map_emis_hc(klon),map_iwp_hc(klon),map_deltaz_hc(klon))
     128      allocate (map_pcld_hc(klon),map_tcld_hc(klon))
     129      allocate (map_emis_hist(klon),map_iwp_hist(klon),map_deltaz_hist(klon))
     130      allocate (map_rad_hist(klon))
     131      allocate (map_ntot(klon),map_hc(klon),map_hist(klon))
     132      allocate (map_Cb(klon),map_ThCi(klon),map_Anv(klon))
     133      allocate (map_emis_Cb(klon),map_pcld_Cb(klon),map_tcld_Cb(klon))
     134      allocate (map_emis_ThCi(klon),map_pcld_ThCi(klon),map_tcld_ThCi(klon))
     135      allocate (map_emis_Anv(klon),map_pcld_Anv(klon),map_tcld_Anv(klon))
     136
    105137    IF (ok_hines) allocate(zustr_gwd_hines(klon), zvstr_gwd_hines(klon))
    106138    IF (.not.ok_hines.and.ok_gwd_rando) &
    … …  
    118150    deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
    119151
     152! Marine
     153! Variables de sortie simulateur AIRS
     154
     155      deallocate (map_prop_hc,map_prop_hist)
     156      deallocate (alt_tropo)
     157      deallocate &
     158     & (map_emis_hc,map_iwp_hc,map_deltaz_hc)
     159      deallocate (map_pcld_hc,map_tcld_hc)
     160      deallocate (map_emis_hist,map_iwp_hist,map_deltaz_hist)
     161      deallocate (map_rad_hist)
     162      deallocate (map_ntot,map_hc,map_hist)
     163      deallocate (map_Cb,map_ThCi,map_Anv)
     164      deallocate (map_emis_Cb,map_pcld_Cb,map_tcld_Cb)
     165      deallocate (map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi)
     166      deallocate (map_emis_Anv,map_pcld_Anv,map_tcld_Anv)
     167
    120168  END SUBROUTINE phys_output_var_end
    121169

  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2561 r2580  
    168168         o_sens_prec_sol_oce, o_sens_prec_sol_sic, &
    169169         o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
    170          o_lat_prec_sol_oce, o_lat_prec_sol_sic
     170         o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
     171! Marine
     172         o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
     173         o_map_deltaz_hc, o_map_pcld_hc, o_map_tcld_hc, &
     174         o_map_emis_hist, o_map_iwp_hist, o_map_deltaz_hist, &
     175         o_map_rad_hist, &
     176         o_map_emis_Cb, o_map_pcld_Cb, o_map_tcld_Cb, &
     177         o_map_emis_ThCi, o_map_pcld_ThCi, o_map_tcld_ThCi, &
     178         o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &
     179         o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &
     180         o_alt_tropo
     181
    171182
    172183    USE phys_state_var_mod, only: pctsrf, paire_ter, rain_fall, snow_fall, &
    … …  
    264275         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
    265276         zustr_gwd_front, zvstr_gwd_front,     &
    266          sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
     277         sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
     278! Marine
     279         map_prop_hc, map_prop_hist, &
     280         map_emis_hc,map_iwp_hc,map_deltaz_hc,&
     281         map_pcld_hc,map_tcld_hc,&
     282         map_emis_hist,map_iwp_hist,map_deltaz_hist,&
     283         map_rad_hist,&
     284         map_ntot,map_hc,map_hist,&
     285         map_Cb,map_ThCi,map_Anv,&
     286         map_emis_Cb,map_pcld_Cb,map_tcld_Cb,&
     287         map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&
     288         map_emis_Anv,map_pcld_Anv,map_tcld_Anv, &
     289         alt_tropo
     290
    267291 
    268292
    … …  
    380404       CALL histwrite_phy(o_aireTER, paire_ter)
    381405!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     406! Simulateur AIRS
     407     IF (ok_airs) then
     408      CALL histwrite_phy(o_alt_tropo,alt_tropo)
     409 
     410      CALL histwrite_phy(o_map_prop_hc,map_prop_hc)
     411      CALL histwrite_phy(o_map_prop_hist,map_prop_hist)
     412
     413      CALL histwrite_phy(o_map_emis_hc,map_emis_hc)
     414      CALL histwrite_phy(o_map_iwp_hc,map_iwp_hc)
     415      CALL histwrite_phy(o_map_deltaz_hc,map_deltaz_hc)
     416      CALL histwrite_phy(o_map_pcld_hc,map_pcld_hc)
     417      CALL histwrite_phy(o_map_tcld_hc,map_tcld_hc)
     418
     419      CALL histwrite_phy(o_map_emis_hist,map_emis_hist)
     420      CALL histwrite_phy(o_map_iwp_hist,map_iwp_hist)
     421      CALL histwrite_phy(o_map_deltaz_hist,map_deltaz_hist)
     422
     423      CALL histwrite_phy(o_map_ntot,map_ntot)
     424      CALL histwrite_phy(o_map_hc,map_hc)
     425      CALL histwrite_phy(o_map_hist,map_hist)
     426
     427      CALL histwrite_phy(o_map_Cb,map_Cb)
     428      CALL histwrite_phy(o_map_ThCi,map_ThCi)
     429      CALL histwrite_phy(o_map_Anv,map_Anv)
     430
     431      CALL histwrite_phy(o_map_emis_Cb,map_emis_Cb)
     432      CALL histwrite_phy(o_map_pcld_Cb,map_pcld_Cb)
     433      CALL histwrite_phy(o_map_tcld_Cb,map_tcld_Cb)
     434
     435      CALL histwrite_phy(o_map_emis_ThCi,map_emis_ThCi)
     436      CALL histwrite_phy(o_map_pcld_ThCi,map_pcld_ThCi)
     437      CALL histwrite_phy(o_map_tcld_ThCi,map_tcld_ThCi)
     438
     439      CALL histwrite_phy(o_map_emis_Anv,map_emis_Anv)
     440      CALL histwrite_phy(o_map_pcld_Anv,map_pcld_Anv)
     441      CALL histwrite_phy(o_map_tcld_Anv,map_tcld_Anv)
     442     ENDIF
     443
    382444       CALL histwrite_phy(o_flat, zxfluxlat)
    383445       CALL histwrite_phy(o_ptstar, ptstar)

  • LMDZ5/trunk/libf/phylmd/physiq_mod.F90

    r2569 r2580  
    39483948#endif
    39493949    ENDIF  !ok_cosp
     3950
     3951
     3952! Marine
     3953
     3954  IF (ok_airs) then
     3955
     3956  IF (itap.eq.1.or.MOD(itap,NINT(freq_airs/dtime)).EQ.0) THEN
     3957  write(*,*) 'je vais appeler simu_airs, ok_airs, freq_airs=', &
     3958     & ok_airs, freq_airs
     3959  call simu_airs(itap,rneb, t_seri, cldemi, fiwc, ref_ice, pphi, pplay, paprs,&
     3960     & map_prop_hc,map_prop_hist,&
     3961     & map_emis_hc,map_iwp_hc,map_deltaz_hc,map_pcld_hc,map_tcld_hc,&
     3962     & map_emis_Cb,map_pcld_Cb,map_tcld_Cb,&
     3963     & map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&
     3964     & map_emis_Anv,map_pcld_Anv,map_tcld_Anv,&
     3965     & map_emis_hist,map_iwp_hist,map_deltaz_hist,map_rad_hist,&
     3966     & map_ntot,map_hc,map_hist,&
     3967     & map_Cb,map_ThCi,map_Anv,&
     3968     & alt_tropo )
     3969  ENDIF
     3970
     3971  ENDIF  ! ok_airs
     3972
     3973
    39503974    ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    39513975    !AA
Note: See TracChangeset for help on using the changeset viewer.