Ignore:
Timestamp:
Apr 6, 2016, 8:51:34 AM (9 years ago)
Author:
emillour
Message:

Venus and Titan GCMs:
Updates in the physics to keep up with updates in LMDZ5 (up to
LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:

  • Adapted makelmdz and makelmdz_fcm script to stop if trying to compile 1d model or newstart or start2archive in parallel.
  • got rid of references to "dimensions.h" in physics. Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead. Only partially done for Titan, because of many hard-coded commons; a necessary first step will be to clean these up (using modules).

EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.TITAN/libf/phytitan/ini_histins.h

    r1056 r1530  
    2121
    2222         CALL histdef(nid_ins, "phis", "Surface geop. height", "-",
    23      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     23     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2424     .                "once",  zsto,zout)
    2525
    2626         CALL histdef(nid_ins, "aire", "Grid area", "-",
    27      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     27     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    2828     .                "once",  zsto,zout)
    2929
    3030ccccccc axe Ls
    3131         CALL histdef(nid_ins, "ls", "Solar longitude", "degrees",
    32      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     32     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3333     .                "inst(X)", zsto,zout)
    3434
    … …  
    3636
    3737         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K",
    38      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     38     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    3939     .                "inst(X)", zsto,zout)
    4040
    4141         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa",
    42      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     42     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4343     .                "inst(X)", zsto,zout)
    4444
    4545c        CALL histdef(nid_ins, "ue", "Zonal energy transport", "-",
    46 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     46c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    4747c    .                "inst(X)", zsto,zout)
    4848
    4949c        CALL histdef(nid_ins, "ve", "Merid energy transport", "-",
    50 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     50c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    5151c    .                "inst(X)", zsto,zout)
    5252
    … …  
    5959
    6060         CALL histdef(nid_ins, "temp", "Air temperature", "K",
    61      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     61     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6262     .                "inst(X)", zsto,zout)
    6363
    6464         CALL histdef(nid_ins, "pres", "Air pressure", "Pa",
    65      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     65     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    6666     .                "inst(X)", zsto,zout)
    6767
    6868         CALL histdef(nid_ins, "geop", "Geopotential height", "m",
    69      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     69     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7070     .                "inst(X)", zsto,zout)
    7171
    7272         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s",
    73      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     73     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7474     .                "inst(X)", zsto,zout)
    7575
    7676         CALL histdef(nid_ins, "vitv", "Meridional wind", "m/s",
    77      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     77     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    7878     .                "inst(X)", zsto,zout)
    7979
    8080         CALL histdef(nid_ins, "vitw", "Vertical wind", "Pa/s",
    81      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     81     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    8282     .                "inst(X)", zsto,zout)
    8383
    8484         CALL histdef(nid_ins, "tops", "Solar rad. at TOA", "W/m2",
    85      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     85     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    8686     .                "inst(X)", zsto,zout)
    8787
    8888c        CALL histdef(nid_ins, "duvdf", "Boundary-layer dU", "m/s2",
    89 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     89c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9090c    .                "inst(X)", zsto,zout)
    9191
    9292c        CALL histdef(nid_ins, "dudyn", "Dynamics dU", "m/s2",
    93 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     93c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    9494c    .                "inst(X)", zsto,zout)
    9595
    … …  
    105105           DO iq=1,nmicro
    106106         CALL histdef(nid_ins, tname(iq), ttext(iq), "n/m2",
    107      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     107     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    108108     .                "inst(X)", zsto,zout)
    109109           ENDDO
    … …  
    112112           DO iq=nmicro+1,nqmax
    113113         CALL histdef(nid_ins, tname(iq), ttext(iq), "ppm",
    114      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     114     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    115115     .                "inst(X)", zsto,zout)
    116116           ENDDO
    … …  
    123123
    124124         CALL histdef(nid_ins, "topl", "IR rad. at TOA", "W/m2",
    125      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     125     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    126126     .                "inst(X)", zsto,zout)
    127127
    128128         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2",
    129      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     129     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    130130     .                "inst(X)", zsto,zout)
    131131
    132132         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2",
    133      .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     133     .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    134134     .                "inst(X)", zsto,zout)
    135135
    … …  
    137137
    138138         CALL histdef(nid_ins, "SWnet", "Net SW flux","W/m2",
    139      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     139     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    140140     .                32, "inst(X)", zsto,zout)
    141141
    142142         CALL histdef(nid_ins, "LWnet", "Net LW flux","W/m2",
    143      .                iim,jj_nb,nhori, klev,1,klev,nvert,
     143     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert,
    144144     .                32, "inst(X)", zsto,zout)
    145145
    … …  
    149149           write(str2,'(i2.2)') k
    150150         CALL histdef(nid_ins,"thv"//str2,"Haze Opa Vis",
    151      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     151     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    152152     .                "ins(X)",zsto,zout)
    153153         ENDDO
    … …  
    156156           write(str2,'(i2.2)') k
    157157         CALL histdef(nid_ins,"thi"//str2,"Haze Opa IR",
    158      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     158     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    159159     .                "ins(X)",zsto,zout)
    160160         ENDDO
    … …  
    165165           write(str2,'(i2.2)') k
    166166         CALL histdef(nid_ins,"khv"//str2,"Haze ext Vis ",
    167      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     167     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    168168     .                "ins(X)",zsto,zout)
    169169         ENDDO
    … …  
    172172           write(str2,'(i2.2)') k
    173173         CALL histdef(nid_ins,"khi"//str2,"Haze ext IR ",
    174      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     174     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    175175     .                "ins(X)",zsto,zout)
    176176         ENDDO
    … …  
    181181           write(str2,'(i2.2)') k
    182182         CALL histdef(nid_ins,"tgv"//str2,"Haze Opa Vis",
    183      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     183     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    184184     .                "ins(X)",zsto,zout)
    185185         ENDDO
    … …  
    188188           write(str2,'(i2.2)') k
    189189         CALL histdef(nid_ins,"tgi"//str2,"Haze Opa IR",
    190      .                "--",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     190     .                "--",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    191191     .                "ins(X)",zsto,zout)
    192192         ENDDO
    … …  
    197197           write(str2,'(i2.2)') k
    198198         CALL histdef(nid_ins,"kgv"//str2,"Haze ext Vis ",
    199      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     199     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    200200     .                "ins(X)",zsto,zout)
    201201         ENDDO
    … …  
    204204           write(str2,'(i2.2)') k
    205205         CALL histdef(nid_ins,"kgi"//str2,"Haze ext IR ",
    206      .                "m-1",iim,jj_nb,nhori,klev,1,klev,nvert,32,
     206     .                "m-1",nbp_lon,jj_nb,nhori,klev,1,klev,nvert,32,
    207207     .                "ins(X)",zsto,zout)
    208208         ENDDO
    … …  
    214214
    215215         CALL histdef(nid_ins, "dtdyn", "Dynamics dT", "K/s",
    216      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     216     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    217217     .                "inst(X)", zsto,zout)
    218218
    219219         CALL histdef(nid_ins, "dtphy", "Physics dT", "K/s",
    220      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     220     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    221221     .                "inst(X)", zsto,zout)
    222222
    223223         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s",
    224      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     224     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    225225     .                "inst(X)", zsto,zout)
    226226
    227227         CALL histdef(nid_ins, "dtajs", "Dry adjust. dT", "K/s",
    228      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     228     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    229229     .                "inst(X)", zsto,zout)
    230230
    231231         CALL histdef(nid_ins, "dtswr", "SW radiation dT", "K/s",
    232      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     232     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    233233     .                "inst(X)", zsto,zout)
    234234
    235235         CALL histdef(nid_ins, "dtlwr", "LW radiation dT", "K/s",
    236      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     236     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    237237     .                "inst(X)", zsto,zout)
    238238
    239239c        CALL histdef(nid_ins, "dtec", "Cinetic dissip dT", "K/s",
    240 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     240c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    241241c    .                "inst(X)", zsto,zout)
    242242
    243243c        CALL histdef(nid_ins, "dvvdf", "Boundary-layer dV", "m/s2",
    244 c    .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
     244c    .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
    245245c    .                "inst(X)", zsto,zout)
    246246
    … …  
    253253c        call histdef(nid_ins, "taux",
    254254c    $         "Zonal wind stress", "Pa", 
    255 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     255c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    256256c    $         "inst(X)", zsto,zout)
    257257
    258258c        call histdef(nid_ins, "tauy",
    259259c    $         "Meridional xind stress", "Pa", 
    260 c    $         iim,jj_nb,nhori, 1,1,1, nvert, 32,
     260c    $         nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    261261c    $         "inst(X)", zsto,zout)
    262262
    263263c        CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-",
    264 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     264c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    265265c    .                "inst(X)", zsto,zout)
    266266
    267267c        CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-",
    268 c    .                iim,jj_nb,nhori, 1,1,1, nvert, 32,
     268c    .                nbp_lon,jj_nb,nhori, 1,1,1, nvert, 32,
    269269c    .                "inst(X)", zsto,zout)
    270270
Note: See TracChangeset for help on using the changeset viewer.