Ignore:
Timestamp:
May 4, 2005, 5:11:29 PM (19 years ago)
Author:
Laurent Fairhead
Message:

Modifications faites à la physique pour la rendre parallele YM
Une branche de travail LMDZ4_par_0 a été créée provisoirement afin de tester
les modifs pleinement avant leurs inclusions dans le tronc principal
LF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4_par_0/libf/phylmd/ini_histins.h

    r524 r634  
    1010         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
    1111c
    12          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    13          DO i = 1, iim
    14             zx_lon(i,1) = rlon(i+1)
    15             zx_lon(i,jjmp1) = rlon(i+1)
    16          ENDDO
     12cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
     13cym         DO i = 1, iim
     14cym            zx_lon(i,1) = rlon(i+1)
     15cym            zx_lon(i,jjmp1) = rlon(i+1)
     16cym         ENDDO
    1717         DO ll=1,klev
    1818            znivsig(ll)=float(ll)
    1919         ENDDO
    20          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
    21          CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
    22      .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     20cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
     21cym         CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:),
     22cym     .                 1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     23cym     .                 nhori, nid_ins)
     24         CALL histbeg_phy("histins", itau_phy, zjulian, dtime,
    2325     .                 nhori, nid_ins)
    2426         write(*,*)'Inst ', itau_phy, zjulian
    … …  
    3032c
    3133         CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
    32      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     34     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    3335     .                "once", zsto,zout)
    3436c
    3537         CALL histdef(nid_ins, "aire", "Grid area", "-",
    36      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     38     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    3739     .                "once", zsto,zout)
    3840c
    … …  
    4042c
    4143        CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
    42      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     44     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4345     .                "inst(X)", zsto,zout)
    4446c
    4547         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K",
    46      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     48     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4749     .                "inst(X)", zsto,zout)
    4850c
    4951         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg",
    50      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     52     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    5153     .                "inst(X)", zsto,zout)
    5254c
    5355         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s",
    54      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     56     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    5557     .                "inst(X)", zsto,zout)
    5658c
    5759         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s",
    58      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     60     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    5961     .                "inst(X)", zsto,zout)
    6062c
    6163        CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
    62      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     64     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    6365     .                "inst(X)", zsto,zout)
    6466c
    6567         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day",
    66      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     68     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    6769     .                "inst(X)", zsto,zout)
    6870c
    6971         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day",
    70      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     72     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7173     .                "inst(X)", zsto,zout)
    7274c
    7375         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
    74      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     76     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7577     .                "inst(X)", zsto,zout)
    7678c
    7779         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
    78      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     80     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    7981     .                "inst(X)", zsto,zout)
    8082c
    8183         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",
    8284     .                "kg/(s*m2)",
    83      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     85     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    8486     .                "inst(X)", zsto,zout)
    8587c
    8688         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)",
    87      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     89     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    8890     .                "inst(X)", zsto,zout)
    8991c
    9092         CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
    91      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     93     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    9294     .                "inst(X)", zsto,zout)
    9395c
    9496         CALL histdef(nid_ins, "topl", "OLR", "W/m2",
    95      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     97     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    9698     .                "inst(X)", zsto,zout)
    9799c
    98100         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)",
    99      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     101     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    100102     .                "inst(X)", zsto,zout)
    101103c
    102104         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
    103      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     105     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    104106     .                "inst(X)", zsto,zout)
    105107c
    106108         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
    107      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     109     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    108110     .                "inst(X)", zsto,zout)
    109111c
    110112         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",
    111      .                "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32,
     113     .                "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    112114     .                "inst(X)", zsto,zout)
    113115c
    114116         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2",
    115      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     117     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    116118     .                "inst(X)", zsto,zout)
    117119c
    118120         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2",
    119      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     121     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    120122     .                "inst(X)", zsto,zout)
    121123c
    122124         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2",
    123      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     125     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    124126     .                "inst(X)", zsto,zout)
    125127c
    126128      CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s",
    127      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     129     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    128130     .                "inst(X)", zsto,zout)
    129131c
    130132      CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s",
    131      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     133     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    132134     .                "inst(X)", zsto,zout)
    133135c
    134136      CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s",
    135      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     137     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    136138     .                "inst(X)", zsto,zout)
    137139c
    138140      CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s",
    139      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     141     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    140142     .                "inst(X)", zsto,zout)
    141143
    … …  
    144146           call histdef(nid_ins, "pourc_"//clnsurf(nsrf),
    145147     $         "Fraction"//clnsurf(nsrf), "W/m2", 
    146      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     148     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    147149     $         "inst(X)", zsto,zout)
    148150
    149151           call histdef(nid_ins, "sens_"//clnsurf(nsrf),
    150152     $         "Sensible heat flux "//clnsurf(nsrf), "W/m2", 
    151      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     153     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    152154     $         "inst(X)", zsto,zout)
    153155c
    154156           call histdef(nid_ins, "tsol_"//clnsurf(nsrf),
    155157     $         "Surface Temperature"//clnsurf(nsrf), "W/m2", 
    156      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     158     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    157159     $         "inst(X)", zsto,zout)
    158160c
    159161           call histdef(nid_ins, "lat_"//clnsurf(nsrf),
    160162     $         "Latent heat flux "//clnsurf(nsrf), "W/m2", 
    161      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     163     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    162164     $         "inst(X)", zsto,zout)
    163165C
    164166           call histdef(nid_ins, "taux_"//clnsurf(nsrf),
    165167     $         "Zonal wind stress"//clnsurf(nsrf),"Pa",
    166      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     168     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    167169     $         "inst(X)", zsto,zout)
    168170
    169171           call histdef(nid_ins, "tauy_"//clnsurf(nsrf),
    170172     $         "Meridional xind stress "//clnsurf(nsrf), "Pa", 
    171      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     173     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    172174     $         "inst(X)", zsto,zout)
    173175c
    174176           call histdef(nid_ins, "albe_"//clnsurf(nsrf),
    175177     $         "Albedo "//clnsurf(nsrf), "-", 
    176      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     178     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    177179     $         "inst(X)", zsto,zout)
    178180c
    179181           call histdef(nid_ins, "rugs_"//clnsurf(nsrf),
    180182     $         "rugosite "//clnsurf(nsrf), "-", 
    181      $         iim,jjmp1,nhori, 1,1,1, -99, 32,
     183     $         iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    182184     $         "inst(X)", zsto,zout)
    183185CXXX
    184186         END DO
    185187         CALL histdef(nid_ins, "rugs", "rugosity", "-",
    186      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     188     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    187189     .                "inst(X)", zsto,zout)
    188190
    189191c
    190192         CALL histdef(nid_ins, "albs", "Surface albedo", "-",
    191      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     193     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    192194     .                "inst(X)", zsto,zout)
    193195         CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-",
    194      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     196     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    195197     .                "inst(X)", zsto,zout)
    196198c
    … …  
    199201c
    200202         CALL histdef(nid_ins, "temp", "Temperature", "K",
    201      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     203     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    202204     .                "inst(X)", zsto,zout)
    203205c
    204206         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
    205      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     207     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    206208     .                "inst(X)", zsto,zout)
    207209c
    208210         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s",
    209      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     211     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    210212     .                "inst(X)", zsto,zout)
    211213c
    212214         CALL histdef(nid_ins, "geop", "Geopotential height", "m",
    213      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     215     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    214216     .                "inst(X)", zsto,zout)
    215217c
    216218         CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
    217      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     219     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    218220     .                "inst(X)", zsto,zout)
    219221c
    220222         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
    221      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     223     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    222224     .                "inst(X)", zsto,zout)
    223225c
    224226         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s",
    225      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     227     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    226228     .                "inst(X)", zsto,zout)
    227229c
Note: See TracChangeset for help on using the changeset viewer.