Ignore:
Timestamp:
Dec 19, 2002, 5:46:39 PM (22 years ago)
Author:
lmdzadmin
Message:

Convergence avec la version de Ionela dec 2002

YOMCST.? : suppression RI0 (IM)
albedo.F : facteur 1.2 sur le nouveau calcul (IM)
clesphys.h : rajout de différentes ctes (concentration des gaz) (IM)
clmain.F : separation des flux LW, SW (JLD)

remplace qsurf par yqsol (IM)

conf_phys.F90 : rajout de différentes ctes (gaz + orbite) (IM)
convect3.F : DPINV+SIGD*0.5*(EVAP(1)+EVAP(2)) (SBL)
cv3_routines.F:
cvparam3.h : compatibilite avec conema3 TEMPORAIRE (FH)
phyetat0.F : lecture de co2_ppm et solaire pour tests de coherence
phyredem.F : co2_ppm et solaire passé en common
physiq.F : separation flux LW, SW

rajout diagnostiques (slp, w500)
suppression iflag_con = 4
clwcon0=qcondc (FH)
position dU "ENDIF ! ok_cvl"

radlwsw.F : passage des concentrations gaz dans un common (IM)

PEMIS(i) = 1.0 (JLD pour cohérence ORCHIDEE)

stdlevvar.F90 :
suphec.F : suppression init. des ctes orbitales (IM)

nouvelles E/S (ini_hist..., write_hist...)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h

    r418 r433  
    5353     .                "ave(X)", zsto,zout)
    5454c
    55          CALL histdef(nid_day, "tter", "Surface Temperature", "K",
    56      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    57      .                "ave(X)", zsto,zout)
    58 c
    59          CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
    60      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    61      .                "ave(X)", zsto,zout)
    62 c
    63          CALL histdef(nid_day, "toce", "Surface Temperature", "K",
    64      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    65      .                "ave(X)", zsto,zout)
    66 c
    67          CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
    68      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    69      .                "ave(X)", zsto,zout)
    70 c
    71 cccIM
    7255c
    7356         CALL histdef(nid_day, "t2m", "Temperature 2m", "K",
    … …  
    7558     .                "ave(X)", zsto,zout)
    7659c
    77          CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
    78      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    79      .                "ave(X)", zsto,zout)
    80 c
    81          CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
    82      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    83      .                "ave(X)", zsto,zout)
    84 c
    85          CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
    86      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    87      .                "ave(X)", zsto,zout)
    88 c
    89          CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
    90      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    91      .                "ave(X)", zsto,zout)
    92 c
     60         CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",
     61     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     62     .                "ave(X)", zsto,zout)
     63c
     64         CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
     65     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     66     .                "ave(X)", zsto,zout)
     67c
     68         CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
     69     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     70     .                "ave(X)", zsto,zout)
     71c
     72         CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
     73     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     74     .                "ave(X)", zsto,zout)
     75c
     76         CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
     77     .                , "kg/s",
     78     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     79     .                "ave(X)", zsto,zout)
     80c
     81         CALL histdef(nid_day, "snow", "Snow fall", "kg/s",
     82     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     83     .                "ave(X)", zsto,zout)
     84c
     85         CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
     86     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     87     .                "ave(X)", zsto,zout)
     88c
     89         CALL histdef(nid_day, "evap", "Evaporation", "kg/s",
     90     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     91     .                "ave(X)", zsto,zout)
     92c
     93         CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
     94     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     95     .                "ave(X)", zsto,zout)
     96c
     97         CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
     98     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     99     .                "ave(X)", zsto,zout)
     100c
     101         CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
     102     .                "W/m2",
     103     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     104     .                "ave(X)", zsto,zout)
     105c
     106         CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
     107     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     108     .                "ave(X)", zsto,zout)
     109c
     110         CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
     111     .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     112     .                "ave(X)", zsto,zout)
     113c
     114         CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
     115     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     116     .                "ave(X)", zsto,zout)
     117c
     118         CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
     119     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     120     .                "ave(X)", zsto,zout)
     121c
     122         CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",
     123     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     124     .                "ave(X)", zsto,zout)
     125c
     126c        CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",
     127c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     128c    .                "ave(X)", zsto,zout)
     129c
     130c        CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",
     131c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     132c    .                "ave(X)", zsto,zout)
     133c
     134c        CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
     135c    .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     136c    .                "ave(X)", zsto,zout)
     137c
     138         CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
     139     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     140     .                "ave(X)", zsto,zout)
     141c
     142         CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
     143     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     144     .                "ave(X)", zsto,zout)
     145c
     146         CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
     147     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     148     .                "ave(X)", zsto,zout)
     149c
     150         CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
     151     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     152     .                "ave(X)", zsto,zout)
     153c
     154         CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",
     155     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     156     .                "ave(X)", zsto,zout)
     157c
     158         CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2",
     159     .                iim,jjmp1,nhori, 1,1,1,-99,
     160     .                32, "ave(X)", zsto,zout)
     161c
     162         CALL histdef(nid_day, "SWupSFC", "SWup at surface","W/m2",
     163     .                iim,jjmp1,nhori, 1,1,1,-99,
     164     .                32, "ave(X)", zsto,zout)
     165c
     166         CALL histdef(nid_day, "SWdnTOA", "SWdn at TOA","W/m2",
     167     .                iim,jjmp1,nhori, 1,1,1,-99,
     168     .                32, "ave(X)", zsto,zout)
     169c
     170         CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2",
     171     .                iim,jjmp1,nhori, 1,1,1,-99,
     172     .                32, "ave(X)", zsto,zout)
     173c
     174c  Champs dynamiques sur niveaux de pression
     175
     176         CALL histdef(nid_day, "u850", "Zonal wind 850mb", "m/s",
     177     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     178     .                "ave(X)", zsto,zout)
     179
     180         CALL histdef(nid_day, "v850", "Meridional wind 850mb", "m/s",
     181     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     182     .                "ave(X)", zsto,zout)
     183c
     184         CALL histdef(nid_day, "u500", "Zonal wind 500mb", "m/s",
     185     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     186     .                "ave(X)", zsto,zout)
     187
     188         CALL histdef(nid_day, "v500", "Meridional wind 500mb", "m/s",
     189     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     190     .                "ave(X)", zsto,zout)
     191
     192         CALL histdef(nid_day, "u200", "Zonal wind 200mb", "m/s",
     193     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     194     .                "ave(X)", zsto,zout)
     195
     196         CALL histdef(nid_day, "v200", "Meridional wind 200mb", "m/s",
     197     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     198     .                "ave(X)", zsto,zout)
     199
     200         CALL histdef(nid_day,"phi500", "Geopotentiel à 500mb", "m2/s2",
     201     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     202     .                "ave(X)", zsto,zout)
     203
     204         CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa",
     205     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     206     .                "ave(X)", zsto,zout)
     207
     208
     209c=================================================================
     210c   CI DESSOUS, SORTIES A UTILISER EN DEBUG
     211c=================================================================
     212      if ( lev_histday.gt.1 ) then
     213
     214c Champs retires momentannéement en attendant un hypothetique
     215c debugage
     216
    93217         CALL histdef(nid_day, "t2m_min", "Temp. 2m min.",
    94218     .                "K",
    … …  
    101225     .                t2maxcels, zsto,zout)
    102226c
     227c=================================================================
     228c=================================================================
     229c=================================================================
     230c   INITIALISATION DES CHAMPS SUR LES SOUS SURFACES
     231c=================================================================
     232
     233         CALL histdef(nid_day, "tter", "Surface Temperature", "K",
     234     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     235     .                "ave(X)", zsto,zout)
     236c
     237         CALL histdef(nid_day, "tlic", "Surface Temperature", "K",
     238     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     239     .                "ave(X)", zsto,zout)
     240c
     241         CALL histdef(nid_day, "toce", "Surface Temperature", "K",
     242     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     243     .                "ave(X)", zsto,zout)
     244c
     245         CALL histdef(nid_day, "tsic", "Surface Temperature", "K",
     246     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     247     .                "ave(X)", zsto,zout)
     248c
     249cccIM
     250c
     251         CALL histdef(nid_day, "t2mter", "Temp.terre 2m", "K",
     252     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     253     .                "ave(X)", zsto,zout)
     254c
     255         CALL histdef(nid_day, "t2mlic", "Temp.lic 2m", "K",
     256     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     257     .                "ave(X)", zsto,zout)
     258c
     259         CALL histdef(nid_day, "t2moce", "Temp.oce 2m", "K",
     260     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     261     .                "ave(X)", zsto,zout)
     262c
     263         CALL histdef(nid_day, "t2msic", "Temp.sic 2m", "K",
     264     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     265     .                "ave(X)", zsto,zout)
     266c
    103267         CALL histdef(nid_day, "t2mter_min", "Temp.terre 2m min.",
    104268     .                "K",
    … …  
    111275     .                t2maxcels, zsto,zout)
    112276c
    113          CALL histdef(nid_day, "q2m", "Specific humidity", "Kg/Kg",
    114      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    115      .                "ave(X)", zsto,zout)
    116 c
    117          CALL histdef(nid_day, "u10m", "Vent zonal 10m", "m/s",
    118      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    119      .                "ave(X)", zsto,zout)
    120 c
    121277         CALL histdef(nid_day, "u10mter", "Vent zonal ter 10m", "m/s",
    122278     .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    … …  
    135291     .                "ave(X)", zsto,zout)
    136292c
    137          CALL histdef(nid_day, "v10m", "Vent meridien 10m", "m/s",
    138      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    139      .                "ave(X)", zsto,zout)
    140 c
    141293         CALL histdef(nid_day, "v10mter", "Vent meridien ter 10m",
    142294     .                "m/s", iim,jjmp1,nhori, 1,1,1, -99, 32,
    … …  
    155307     .                "ave(X)", zsto,zout)
    156308c
    157          CALL histdef(nid_day, "psol", "Surface Pressure", "Pa",
    158      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    159      .                "ave(X)", zsto,zout)
    160 c
    161          CALL histdef(nid_day, "precip","Precipitation Totale liq+sol"
    162      .                , "kg/s",
    163      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    164      .                "ave(X)", zsto,zout)
    165 c
    166          CALL histdef(nid_day, "snow", "Snow fall", "kg/s",
    167      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    168      .                "ave(X)", zsto,zout)
    169 c
    170          CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2",
    171      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    172      .                "ave(X)", zsto,zout)
    173 c
    174          CALL histdef(nid_day, "evap", "Evaporation", "kg/s",
    175      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    176      .                "ave(X)", zsto,zout)
    177 c
    178          CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",
    179      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    180      .                "ave(X)", zsto,zout)
    181 c
    182          CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2",
    183      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    184      .                "ave(X)", zsto,zout)
    185 c
    186          CALL histdef(nid_day, "sols", "Net Solar rad. at surf.",
    187      .                "W/m2",
    188      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    189      .                "ave(X)", zsto,zout)
    190 c
    191          CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2",
    192      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    193      .                "ave(X)", zsto,zout)
    194 c
    195          CALL histdef(nid_day, "solldown", "Down. IR rad. at surface",
    196      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
    197      .                "ave(X)", zsto,zout)
    198 c
    199          CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2",
    200      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    201      .                "ave(X)", zsto,zout)
    202 c
    203          CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2",
    204      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    205      .                "ave(X)", zsto,zout)
    206 c
    207          CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2",
    208      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    209      .                "ave(X)", zsto,zout)
    210 c
    211          CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa",
    212      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    213      .                "ave(X)", zsto,zout)
    214 c
    215          CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa",
    216      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    217      .                "ave(X)", zsto,zout)
    218309c
    219310CXXX PB flux pour chaque sous surface
    … …  
    264355         END DO
    265356           
    266          CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-",
    267      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    268      .                "ave(X)", zsto,zout)
    269 c
    270          CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-",
    271      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    272      .                "ave(X)", zsto,zout)
    273 c
    274          CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-",
    275      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    276      .                "ave(X)", zsto,zout)
    277 c
    278          CALL histdef(nid_day, "cldh", "High-level cloudiness", "-",
    279      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    280      .                "ave(X)", zsto,zout)
    281 c
    282          CALL histdef(nid_day, "cldt", "Total cloudiness", "-",
    283      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    284      .                "ave(X)", zsto,zout)
    285 c
    286          CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-",
    287      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
    288      .                "ave(X)", zsto,zout)
    289 c
     357c=================================================================
     358c   FIN DES CHAMPS SUR LES SOUS SURFACES
     359c=================================================================
     360c=================================================================
     361c   FIN INITIALISATION DES CHAMPS 3D
     362c=================================================================
     363
    290364c Champs 3D:
    291365c
    … …  
    337411         CALL histend(nid_day)
    338412c
     413c=================================================================
     414c   FIN INITIALISATION DES CHAMPS 3D
     415c=================================================================
     416
     417
     418c=================================================================
     419c=================================================================
     420      endif !  lev_histday.gt.1
     421c=================================================================
     422
     423         CALL histend(nid_day)
     424c
    339425         ndex2d = 0
    340426         ndex3d = 0
    341427c
     428c=================================================================
    342429      ENDIF ! fin de test sur ok_journe
Note: See TracChangeset for help on using the changeset viewer.