Ignore:
Timestamp:
Aug 17, 2006, 5:41:51 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Inclusion des modifs de Y. Meurdesoif pour la version V3
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/V3_test/libf/phylmd/ini_bilKP_ins.h

    r684 r704  
    1111         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
    1212c
    13          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    14          DO i = 1, iim
    15             zx_lon(i,1) = rlon(i+1)
    16             zx_lon(i,jjmp1) = rlon(i+1)
    17          ENDDO
     13cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
     14cym         DO i = 1, iim
     15cym            zx_lon(i,1) = rlon(i+1)
     16cym            zx_lon(i,jjmp1) = rlon(i+1)
     17cym         ENDDO
    1818         DO ll=1,klev
    1919            znivsig(ll)=float(ll)
    2020         ENDDO
    21          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
    22          write(*,*)'zx_lon = ',zx_lon(:,1)
    23          write(*,*)'zx_lat = ',zx_lat(1,:)
     21cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
     22cym         write(*,*)'zx_lon = ',zx_lon(:,1)
     23cym         write(*,*)'zx_lat = ',zx_lat(1,:)
    2424c
    2525cIM 280405 BEG
     
    3030          jmin_ins=1
    3131          jmax_ins=jjmp1
    32           do i=1,iim-1
    33              if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
    34              if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
    35           enddo
    36           do j=1,jjmp1
    37              if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
    38              if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
    39           enddo
     32cym          do i=1,iim-1
     33cym             if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i
     34cym             if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1
     35cym          enddo
     36cym          do j=1,jjmp1
     37cym             if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j
     38cym             if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j
     39cym          enddo
    4040c
    4141          print*,'On stoke le fichier bilKP instantanne sur ',
     
    4848cIM 280405 END
    4949c
    50          IF(1.EQ.0) THEN
    51          CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,
    52      .                zx_lat(1,:),
    53      .                1,iim,1,jjmp1, itau_phy, zjulian, dtime,
    54      .                nhori, nid_bilKPins)
     50cym         IF(1.EQ.0) THEN
     51cym         CALL histbeg("histbilKP_ins", iim,zx_lon(:,1), jjmp1,
     52cym     .                zx_lat(1,:),
     53cym     .                1,iim,1,jjmp1, itau_phy, zjulian, dtime,
     54cym     .                nhori, nid_bilKPins)
    5555         ENDIF
    5656c
     
    5858c
    5959cIM cf. AM 081204 BEG region
    60          CALL histbeg("histbilKP_ins", iim,zx_lon(:,1),
    61      .                 jjmp1,zx_lat(1,:),
    62      .                 imin_ins,imax_ins-imin_ins+1,
    63      .                 jmin_ins,jmax_ins-jmin_ins+1,
    64      .                 itau_phy, zjulian, dtime,
     60cym         CALL histbeg("histbilKP_ins", iim,zx_lon(:,1),
     61cym     .                 jjmp1,zx_lat(1,:),
     62cym     .                 imin_ins,imax_ins-imin_ins+1,
     63cym     .                 jmin_ins,jmax_ins-jmin_ins+1,
     64cym     .                 itau_phy, zjulian, dtime,
     65cym     .                 nhori, nid_bilKPins)
     66         CALL histbeg_phy("histbilKP_ins", itau_phy, zjulian, dtime,
    6567     .                 nhori, nid_bilKPins)
    6668cIM 081204 END
     
    7779         CALL histdef(nid_bilKPins,"ue",
    7880     .   "Zonal energy transport","-",
    79      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     81     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    8082     .                typeval, zsto,zout)
    8183c
    8284         CALL histdef(nid_bilKPins,"ve",
    8385     .   "Merid energy transport","-",
    84      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     86     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    8587     .                typeval, zsto,zout)
    8688c
    8789         CALL histdef(nid_bilKPins,"uq",
    8890     .   "Zonal humidity transport","-",
    89      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     91     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    9092     .                typeval, zsto,zout)
    9193c
    9294         CALL histdef(nid_bilKPins,"vq",
    9395     .   "Merid humidity transport","-",
    94      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     96     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    9597     .                typeval, zsto,zout)
    9698c
     
    99101         CALL histdef(nid_bilKPins, "temp",
    100102     .   "Air temperature", "K",
    101      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     103     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    102104     .                typeval, zsto,zout)
    103105c
    104106         CALL histdef(nid_bilKPins,"ovap",
    105107     .   "Specific humidity","Kg/Kg",
    106      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     108     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    107109     .                typeval, zsto,zout)
    108110c
    109111         CALL histdef(nid_bilKPins,"geop",
    110112     .   "Geopotential height", "m",
    111      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     113     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    112114     .                typeval, zsto,zout)
    113115c
    114116         CALL histdef(nid_bilKPins,"vitu",
    115117     .   "Zonal wind", "m/s",
    116      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     118     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    117119     .                typeval, zsto,zout)
    118120c
    119121         CALL histdef(nid_bilKPins,"vitv",
    120122     .   "Meridional wind", "m/s",
    121      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     123     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    122124     .                typeval, zsto,zout)
    123125c
    124126         CALL histdef(nid_bilKPins, "vitw",
    125127     .   "Vertical wind", "m/s",
    126      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     128     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    127129     .                typeval, zsto,zout)
    128130c
     
    130132     .   "Inter-Layer Air pressure",
    131133     .                "Pa",
    132      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     134     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    133135     .                typeval, zsto,zout)
    134136c
     
    136138     .   "Mean-Layer Air pressure",
    137139     .                "Pa",
    138      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     140     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    139141     .                typeval, zsto,zout)
    140142c
     
    142144     .   "Liquid water content",
    143145     .                "kg/kg",
    144      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     146     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    145147     .                typeval, zsto,zout)
    146148c
    147149         CALL histdef(nid_bilKPins, "dtdyn",
    148150     .   "Dynamics dT", "K/s",
    149      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     151     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    150152     .                typeval, zsto,zout)
    151153c
    152154         CALL histdef(nid_bilKPins, "dqdyn",
    153155     .   "Dynamics dQ", "Kg/Kg/s",
    154      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     156     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    155157     .                typeval, zsto,zout)
    156158c
    157159         CALL histdef(nid_bilKPins, "dtcon",
    158160     .   "Convection dT", "K/s",
    159      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     161     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    160162     .                typeval, zsto,zout)
    161163c
    162164         CALL histdef(nid_bilKPins, "ducon",
    163165     .   "Convection du", "m/s2",
    164      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     166     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    165167     .                typeval, zsto,zout)
    166168c
    167169         CALL histdef(nid_bilKPins, "dvcon",
    168170     .   "Convection dv", "m/s2",
    169      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     171     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    170172     .                typeval, zsto,zout)
    171173c
    172174         CALL histdef(nid_bilKPins,"dqcon",
    173175     .   "Convection dQ","Kg/Kg/s",
    174      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     176     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    175177     .                typeval, zsto,zout)
    176178c
    177179         CALL histdef(nid_bilKPins, "dtlsc",
    178180     .   "Condensation dT", "K/s",
    179      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     181     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    180182     .                typeval, zsto,zout)
    181183c
    182184         CALL histdef(nid_bilKPins,"dqlsc",
    183185     .   "Condensation dQ","Kg/Kg/s",
    184      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     186     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    185187     .                typeval, zsto,zout)
    186188c
    187189         CALL histdef(nid_bilKPins,"dtvdf",
    188190     .   "Boundary-layer dT","K/s",
    189      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     191     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    190192     .                typeval, zsto,zout)
    191193c
     
    193195     .   "Boundary-layer dQ",
    194196     .               "Kg/Kg/s",
    195      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     197     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    196198     .                typeval, zsto,zout)
    197199c
    198200         CALL histdef(nid_bilKPins,"dtajs",
    199201     .   "Ajustement sec dT","K/s",
    200      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     202     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    201203     .                typeval, zsto,zout)
    202204c
     
    204206     .   "Ajustement sec dQ",
    205207     .               "Kg/Kg/s",
    206      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     208     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    207209     .                typeval, zsto,zout)
    208210c
    209211         CALL histdef(nid_bilKPins,"dteva",
    210212     .   "Reevaporation dT","K/s",
    211      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     213     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    212214     .                typeval, zsto,zout)
    213215c
     
    215217     .   "Reevaporation dQ",
    216218     .                "Kg/Kg/s",
    217      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     219     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    218220     .                typeval, zsto,zout)
    219221
     
    221223         CALL histdef(nid_bilKPins, "dtswr",
    222224     .   "SW radiation dT", "K/s",
    223      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     225     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    224226     .                typeval, zsto,zout)
    225227c
    226228         CALL histdef(nid_bilKPins, "dtsw0",
    227229     .   "SW radiation dT", "K/s",
    228      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     230     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    229231     .                typeval, zsto,zout)
    230232c
    231233         CALL histdef(nid_bilKPins, "dtlwr",
    232234     .   "LW radiation dT", "K/s",
    233      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     235     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    234236     .                typeval, zsto,zout)
    235237c
    236238         CALL histdef(nid_bilKPins, "dtlw0",
    237239     .   "LW radiation dT", "K/s",
    238      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     240     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    239241     .                typeval, zsto,zout)
    240242c
    241243         CALL histdef(nid_bilKPins,"duvdf",
    242244     .   "Boundary-layer dU","m/s2",
    243      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     245     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    244246     .                typeval, zsto,zout)
    245247c
    246248         CALL histdef(nid_bilKPins,"dvvdf",
    247249     .   "Boundary-layer dV","m/s2",
    248      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     250     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    249251     .                typeval, zsto,zout)
    250252c
     
    253255         CALL histdef(nid_bilKPins, "duoli",
    254256     .   "Orography dU", "m/s2",
    255      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     257     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    256258     .                typeval, zsto,zout)
    257259c
    258260         CALL histdef(nid_bilKPins, "dvoli",
    259261     .   "Orography dV", "m/s2",
    260      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     262     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    261263     .                typeval, zsto,zout)
    262264c
     
    266268         CALL histdef(nid_bilKPins, "duphy",
    267269     .   "Physiq dU","-",
    268      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     270     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    269271     .                typeval, zsto,zout)
    270272C
    271273         CALL histdef(nid_bilKPins, "dvphy",
    272274     .   "Physiq dV","-",
    273      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     275     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    274276     .                typeval, zsto,zout)
    275277C
    276278         CALL histdef(nid_bilKPins, "dtphy",
    277279     .   "Physiq dT","-",
    278      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     280     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    279281     .                typeval, zsto,zout)
    280282C
    281283         CALL histdef(nid_bilKPins, "dqphy",
    282284     .   "Physiq dQ","-",
    283      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     285     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    284286     .                typeval, zsto,zout)
    285287C
    286288         CALL histdef(nid_bilKPins, "dqlphy",
    287289     .   "Physiq dQl","-",
    288      .                iim,jjmp1,nhori, klev,1,klev,nvert, 32,
     290     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    289291     .                typeval, zsto,zout)
    290292cIM 280405 BEG
     
    303305          CALL histdef(nid_bilKPins, "u"//bb2,
    304306     .                 "Zonal wind "//bb2//"mb","m/s",
    305      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     307     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    306308     .                typeval, zsto,zout)
    307309c
    308310          CALL histdef(nid_bilKPins, "v"//bb2,
    309311     .                 "Meridional wind "//bb2//"mb","m/s",
    310      .                iim,jjmp1,nhori, 1,1,1, -99, 32,
     312     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    311313     .                typeval, zsto,zout)
    312314c
Note: See TracChangeset for help on using the changeset viewer.