Ignore:
Timestamp:
Jun 25, 2014, 1:19:59 PM (11 years ago)
Author:
emillour
Message:

Common dynamics:
Some updates to keep up with LMDZ5 Earth model evolution (up to LMDZ5 rev 1955).
Main change is the introduction of a "dyn3d_common" directory
to store files common to dyn3d and dyn3dpar.
See file "DOC/chantiers/commit_importants.log" for detailed list
of changes. These changes do not change results on test cases.
EM

Location:
trunk/LMDZ.COMMON/libf/dyn3d_common
Files:
1 added
129 moved

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d_common/defrun.F

    r1299 r1300  
    370370        WRITE(tapeout,*)' La valeur de dzoomx passee par run.def est dif
    371371     *ferente de celle lue sur le fichier  start '
    372         CALL ABORT
     372        CALL ABORT_gcm("defrun", "", 1)
    373373       ENDIF
    374374
    … …  
    376376        WRITE(tapeout,*)' La valeur de dzoomy passee par run.def est dif
    377377     *ferente de celle lue sur le fichier  start '
    378         CALL ABORT
     378        CALL ABORT_gcm("defrun", "", 1)
    379379       ENDIF
    380380
    … …  
    382382        WRITE(6,*)' La valeur de taux passee par run.def est differente
    383383     *  de celle lue sur le fichier  start '
    384         CALL ABORT
     384        CALL ABORT_gcm("defrun", "", 1)
    385385       ENDIF
    386386
    … …  
    388388        WRITE(6,*)' La valeur de tauy passee par run.def est differente
    389389     *  de celle lue sur le fichier  start '
    390         CALL ABORT
     390        CALL ABORT_gcm("defrun", "", 1)
    391391       ENDIF
    392392

  • trunk/LMDZ.COMMON/libf/dyn3d_common/dynetat0.F

    r1299 r1300  
    8282        write(lunout,*)'dynetat0: Pb d''ouverture du fichier start.nc'
    8383        write(lunout,*)trim(nf90_strerror(ierr))
    84         CALL ABORT
     84        CALL ABORT_gcm("dynetat0", "", 1)
    8585      ENDIF
    8686
    … …  
    9090         write(lunout,*)"dynetat0: Le champ <controle> est absent"
    9191         write(lunout,*)trim(nf90_strerror(ierr))
    92          CALL abort
     92         CALL ABORT_gcm("dynetat0", "", 1)
    9393      ENDIF
    9494      ierr = nf90_get_var(nid, nvarid, tab_cntrl)
    … …  
    9696         write(lunout,*)"dynetat0: Lecture echoue pour <controle>"
    9797         write(lunout,*)trim(nf90_strerror(ierr))
    98          CALL abort
     98         CALL ABORT_gcm("dynetat0", "", 1)
    9999      ENDIF
    100100
    … …  
    190190         write(lunout,*)"dynetat0: Le champ <rlonu> est absent"
    191191         write(lunout,*)trim(nf90_strerror(ierr))
    192          CALL abort
     192         CALL ABORT_gcm("dynetat0", "", 1)
    193193      ENDIF
    194194      ierr = nf90_get_var(nid, nvarid, rlonu)
    … …  
    196196         write(lunout,*)"dynetat0: Lecture echouee pour <rlonu>"
    197197         write(lunout,*)trim(nf90_strerror(ierr))
    198          CALL abort
     198         CALL ABORT_gcm("dynetat0", "", 1)
    199199      ENDIF
    200200
    … …  
    203203         write(lunout,*)"dynetat0: Le champ <rlatu> est absent"
    204204         write(lunout,*)trim(nf90_strerror(ierr))
    205          CALL abort
     205         CALL ABORT_gcm("dynetat0", "", 1)
    206206      ENDIF
    207207      ierr = nf90_get_var(nid, nvarid, rlatu)
    … …  
    209209         write(lunout,*)"dynetat0: Lecture echouee pour <rlatu>"
    210210         write(lunout,*)trim(nf90_strerror(ierr))
    211          CALL abort
     211         CALL ABORT_gcm("dynetat0", "", 1)
    212212      ENDIF
    213213
    … …  
    216216         write(lunout,*)"dynetat0: Le champ <rlonv> est absent"
    217217         write(lunout,*)trim(nf90_strerror(ierr))
    218          CALL abort
     218         CALL ABORT_gcm("dynetat0", "", 1)
    219219      ENDIF
    220220      ierr = nf90_get_var(nid, nvarid, rlonv)
    … …  
    222222         write(lunout,*)"dynetat0: Lecture echouee pour <rlonv>"
    223223         write(lunout,*)trim(nf90_strerror(ierr))
    224          CALL abort
     224         CALL ABORT_gcm("dynetat0", "", 1)
    225225      ENDIF
    226226
    … …  
    229229         write(lunout,*)"dynetat0: Le champ <rlatv> est absent"
    230230         write(lunout,*)trim(nf90_strerror(ierr))
    231          CALL abort
     231         CALL ABORT_gcm("dynetat0", "", 1)
    232232      ENDIF
    233233      ierr = nf90_get_var(nid, nvarid, rlatv)
    … …  
    235235         write(lunout,*)"dynetat0: Lecture echouee pour rlatv"
    236236         write(lunout,*)trim(nf90_strerror(ierr))
    237          CALL abort
     237         CALL ABORT_gcm("dynetat0", "", 1)
    238238      ENDIF
    239239
    … …  
    242242         write(lunout,*)"dynetat0: Le champ <cu> est absent"
    243243         write(lunout,*)trim(nf90_strerror(ierr))
    244          CALL abort
     244         CALL ABORT_gcm("dynetat0", "", 1)
    245245      ENDIF
    246246      ierr = nf90_get_var(nid, nvarid, cu)
    … …  
    248248         write(lunout,*)"dynetat0: Lecture echouee pour <cu>"
    249249         write(lunout,*)trim(nf90_strerror(ierr))
    250          CALL abort
     250         CALL ABORT_gcm("dynetat0", "", 1)
    251251      ENDIF
    252252
    … …  
    255255         write(lunout,*)"dynetat0: Le champ <cv> est absent"
    256256         write(lunout,*)trim(nf90_strerror(ierr))
    257          CALL abort
     257         CALL ABORT_gcm("dynetat0", "", 1)
    258258      ENDIF
    259259      ierr = nf90_get_var(nid, nvarid, cv)
    … …  
    261261         write(lunout,*)"dynetat0: Lecture echouee pour <cv>"
    262262         write(lunout,*)trim(nf90_strerror(ierr))
    263          CALL abort
     263         CALL ABORT_gcm("dynetat0", "", 1)
    264264      ENDIF
    265265
    … …  
    268268         write(lunout,*)"dynetat0: Le champ <aire> est absent"
    269269         write(lunout,*)trim(nf90_strerror(ierr))
    270          CALL abort
     270         CALL ABORT_gcm("dynetat0", "", 1)
    271271      ENDIF
    272272      ierr = nf90_get_var(nid, nvarid, aire)
    … …  
    274274         write(lunout,*)"dynetat0: Lecture echouee pour <aire>"
    275275         write(lunout,*)trim(nf90_strerror(ierr))
    276          CALL abort
     276         CALL ABORT_gcm("dynetat0", "", 1)
    277277      ENDIF
    278278
    … …  
    281281         write(lunout,*)"dynetat0: Le champ <phisinit> est absent"
    282282         write(lunout,*)trim(nf90_strerror(ierr))
    283          CALL abort
     283         CALL ABORT_gcm("dynetat0", "", 1)
    284284      ENDIF
    285285      ierr = nf90_get_var(nid, nvarid, phis)
    … …  
    287287         write(lunout,*)"dynetat0: Lecture echouee pour <phisinit>"
    288288         write(lunout,*)trim(nf90_strerror(ierr))
    289          CALL abort
     289         CALL ABORT_gcm("dynetat0", "", 1)
    290290      ENDIF
    291291
    … …  
    299299           write(lunout,*)"dynetat0: Le champ <Time> est absent"
    300300           write(lunout,*)trim(nf90_strerror(ierr))
    301            CALL abort
     301           CALL ABORT_gcm("dynetat0", "", 1)
    302302        ENDIF
    303303        ! Get the length of the "Time" dimension
    … …  
    311311           write(lunout,*)"dynetat0: Lecture echouee <Time>"
    312312           write(lunout,*)trim(nf90_strerror(ierr))
    313            CALL abort
     313           CALL ABORT_gcm("dynetat0", "", 1)
    314314        ENDIF
    315315      ELSE   
    … …  
    324324           write(lunout,*)"dynetat0: Lecture echouee <temps>"
    325325           write(lunout,*)trim(nf90_strerror(ierr))
    326            CALL abort
     326           CALL ABORT_gcm("dynetat0", "", 1)
    327327        ENDIF
    328328      ENDIF
    … …  
    346346             PRINT*, time(i)
    347347          ENDDO
    348           CALL abort
     348          CALL ABORT_gcm("dynetat0", "", 1)
    349349        ENDIF
    350350      ENDIF
    … …  
    379379         write(lunout,*)"dynetat0: Le champ <vcov> est absent"
    380380         write(lunout,*)trim(nf90_strerror(ierr))
    381          CALL abort
     381         CALL ABORT_gcm("dynetat0", "", 1)
    382382      ENDIF
    383383      ierr=nf90_get_var(nid,nvarid,vcov,corner,edges)
    … …  
    385385         write(lunout,*)"dynetat0: Lecture echouee pour <vcov>"
    386386         write(lunout,*)trim(nf90_strerror(ierr))
    387          CALL abort
     387         CALL ABORT_gcm("dynetat0", "", 1)
    388388      ENDIF
    389389
    … …  
    401401         write(lunout,*)"dynetat0: Le champ <ucov> est absent"
    402402         write(lunout,*)trim(nf90_strerror(ierr))
    403          CALL abort
     403         CALL ABORT_gcm("dynetat0", "", 1)
    404404      ENDIF
    405405      ierr=nf90_get_var(nid,nvarid,ucov,corner,edges)
    … …  
    407407         write(lunout,*)"dynetat0: Lecture echouee pour <ucov>"
    408408         write(lunout,*)trim(nf90_strerror(ierr))
    409          CALL abort
     409         CALL ABORT_gcm("dynetat0", "", 1)
    410410      ENDIF
    411411 
    … …  
    415415         write(lunout,*)"dynetat0: Le champ <teta> est absent"
    416416         write(lunout,*)trim(nf90_strerror(ierr))
    417          CALL abort
     417         CALL ABORT_gcm("dynetat0", "", 1)
    418418      ENDIF
    419419      ierr=nf90_get_var(nid,nvarid,teta,corner,edges)
    … …  
    421421         write(lunout,*)"dynetat0: Lecture echouee pour <teta>"
    422422         write(lunout,*)trim(nf90_strerror(ierr))
    423          CALL abort
     423         CALL ABORT_gcm("dynetat0", "", 1)
    424424      ENDIF
    425425
    … …  
    439439     &                                //trim(tname(iq))
    440440            write(lunout,*)trim(nf90_strerror(ierr))
    441             CALL abort
     441            CALL ABORT_gcm("dynetat0", "", 1)
    442442          ENDIF
    443443        ENDIF
    … …  
    450450         write(lunout,*)"dynetat0: Le champ <masse> est absent"
    451451         write(lunout,*)trim(nf90_strerror(ierr))
    452          CALL abort
     452         CALL ABORT_gcm("dynetat0", "", 1)
    453453      ENDIF
    454454      ierr=nf90_get_var(nid,nvarid,masse,corner,edges)
    … …  
    456456         write(lunout,*)"dynetat0: Lecture echouee pour <masse>"
    457457         write(lunout,*)trim(nf90_strerror(ierr))
    458          CALL abort
     458         CALL ABORT_gcm("dynetat0", "", 1)
    459459      ENDIF
    460460
    … …  
    470470         write(lunout,*)"dynetat0: Le champ <ps> est absent"
    471471         write(lunout,*)trim(nf90_strerror(ierr))
    472          CALL abort
     472         CALL ABORT_gcm("dynetat0", "", 1)
    473473      ENDIF
    474474      ierr=nf90_get_var(nid,nvarid,ps,corner,edges)
    … …  
    476476         write(lunout,*)"dynetat0: Lecture echouee pour <ps>"
    477477         write(lunout,*)trim(nf90_strerror(ierr))
    478          CALL abort
     478         CALL ABORT_gcm("dynetat0", "", 1)
    479479      ENDIF
    480480

  • trunk/LMDZ.COMMON/libf/dyn3d_common/dynredem.F

    r1299 r1300  
    158158     &                  //trim(fichnom)
    159159         write(lunout,*)' ierr = ', ierr
    160          CALL ABORT
     160         CALL ABORT_GCM("DYNREDEM0", "", 1)
    161161      ENDIF
    162162c
    … …  
    636636      IF (ierr .NE. NF_NOERR) THEN
    637637         PRINT*, "dynredem1: Pb. d ouverture "//trim(fichnom)
    638          CALL abort
     638         CALL abort_gcm("dynredem1", "", 1)
    639639      ENDIF
    640640

  • trunk/LMDZ.COMMON/libf/dyn3d_common/fxhyp.F

    r1299 r1300  
    178178        WRITE(6,*)'Modifier les valeurs de  grossismx ,tau ou dzoomx ',
    179179     , ' et relancer ! ***  '
    180         CALL ABORT
     180        CALL ABORT_GCM("FXHYP", "", 1)
    181181       ENDIF
    182182c

  • trunk/LMDZ.COMMON/libf/dyn3d_common/fyhyp.F

    r1299 r1300  
    162162        WRITE(6,*)'Modifier les valeurs de  grossismy ,tauy ou dzoomy',
    163163     , ' et relancer ! ***  '
    164         CALL ABORT
     164        CALL ABORT_GCM("FYHYP", "", 1)
    165165
    166166       ENDIF

  • trunk/LMDZ.COMMON/libf/dyn3d_common/grid_atob.F

    r1299 r1300  
    5555      IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
    5656         PRINT*, 'imar ou jmar trop grand', imar, jmar
    57          CALL ABORT
     57         CALL ABORT_GCM("", "", 1)
    5858      ENDIF
    5959c
    … …  
    119119         ELSE
    120120         PRINT*, 'probleme,i,j=', i,j
    121 ccc         CALL ABORT
     121ccc         CALL ABORT_GCM("", "", 1)
    122122         CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans)
    123123#ifdef CRAY
    … …  
    193193      IF (imar.GT.400 .OR. jmar.GT.200) THEN
    194194         PRINT*, 'imar ou jmar trop grand', imar, jmar
    195          CALL ABORT
     195         CALL ABORT_GCM("", "", 1)
    196196      ENDIF
    197197c
    198198      IF (imdep.GT.400 .OR. jmdep.GT.200) THEN
    199199         PRINT*, 'imdep ou jmdep trop grand', imdep, jmdep
    200          CALL ABORT
     200         CALL ABORT_GCM("", "", 1)
    201201      ENDIF
    202202c
    … …  
    258258            PRINT*, 'Probleme grave,i,j,indx,indy=',
    259259     .              i,j,indx(i,j),indy(i,j)
    260             CALL abort
     260            CALL abort_gcm("", "", 1)
    261261         ENDIF
    262262      ENDDO
    … …  
    309309         ELSE
    310310         PRINT*, 'probleme,i,j=', i,j
    311          CALL ABORT
     311         CALL ABORT_GCM("", "", 1)
    312312         ENDIF
    313313      ENDDO
    … …  
    345345      IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
    346346         PRINT*, 'imar ou jmar trop grand', imar, jmar
    347          CALL ABORT
     347         CALL ABORT_GCM("", "", 1)
    348348      ENDIF
    349349c
    … …  
    409409         ELSE
    410410         PRINT*, 'probleme,i,j=', i,j
    411          CALL ABORT
     411         CALL ABORT_GCM("", "", 1)
    412412         ENDIF
    413413      ENDDO
    … …  
    452452      IF (imar.GT.400 .OR. jmar.GT.400) THEN
    453453         PRINT*, 'imar ou jmar trop grand', imar, jmar
    454          CALL ABORT
     454         CALL ABORT_GCM("", "", 1)
    455455      ENDIF
    456456c
    … …  
    513513         ELSE
    514514            PRINT*, 'probleme,i,j=', i,j
    515 ccc            CALL ABORT
     515ccc            CALL ABORT_GCM("", "", 1)
    516516         CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans)
    517517#ifdef CRAY
    … …  
    577577      IF (imar.GT.400 .OR. jmar.GT.400) THEN
    578578         PRINT*, 'imar ou jmar trop grand', imar, jmar
    579          CALL ABORT
     579         CALL ABORT_GCM("", "", 1)
    580580      ENDIF
    581581c
    … …  
    642642         ELSE
    643643           PRINT*, 'probleme,i,j=', i,j
    644 ccc           CALL ABORT
     644ccc           CALL ABORT_GCM("", "", 1)
    645645         CALL dist_sphe(x(i),y(j),xdata,ydata,imdep,jmdep,distans)
    646646#ifdef CRAY
    … …  
    712712      IF (immod.GT.2200 .OR. jmmod.GT.1100) THEN
    713713         PRINT*, 'immod ou jmmod trop grand', immod, jmmod
    714          CALL ABORT
     714         CALL ABORT_GCM("", "", 1)
    715715      ENDIF
    716716c
    … …  
    801801         ELSE
    802802         PRINT*, 'probleme,i,j=', i,j
    803          CALL ABORT
     803         CALL ABORT_GCM("", "", 1)
    804804         ENDIF
    805805      ENDDO
    … …  
    875875         ELSE
    876876         PRINT*, 'probleme,i,j=', i,j
    877 ccc         CALL ABORT
     877ccc         CALL ABORT_GCM("", "", 1)
    878878         CALL dist_sphe(xmod(i),ymod(j),xtmp,ytmp,imtmp,jmtmp,distans)
    879879#ifdef CRAY

  • trunk/LMDZ.COMMON/libf/dyn3d_common/inigrads.F

    r1299 r1300  
    1313      real xmin,xmax,ymin,ymax
    1414
    15       character file*10,titlel*40
     15      character(len=*),intent(in) :: file
     16      character(len=*),intent(in) :: titlel
    1617
    1718#include "gradsdef.h"

  • trunk/LMDZ.COMMON/libf/dyn3d_common/juldate.F

    r1299 r1300  
    99        implicit real (a-h,o-z)
    1010        frac=((os/60.+om)/60.+oh)/24.
    11         ojou=dfloat(ijou)+frac
    12             year=dfloat(ian)
    13             rmon=dfloat(imoi)
     11        ojou=dble(ijou)+frac
     12            year=dble(ian)
     13            rmon=dble(imoi)
    1414        if (imoi .le. 2) then
    1515            year=year-1.

  • trunk/LMDZ.COMMON/libf/dyn3d_common/pres2lev_mod.F90

    r1299 r1300  
    2727  INTEGER, INTENT(IN) :: lmn ! dimensions nouvelle couches
    2828 
     29  INTEGER, INTENT(IN) :: ni,nj ! nombre de point horizontal
    2930  REAL, INTENT(IN) :: po(ni*nj,lmo) ! niveau de pression ancienne grille
    3031  REAL, INTENT(IN) :: pn(ni*nj,lmn) ! niveau de pression nouvelle grille
    31 
    32   INTEGER, INTENT(IN) :: ni,nj ! nombre de point horizontal
    3332
    3433  REAL, INTENT(IN)  :: varo(ni*nj,lmo) ! var dans l'ancienne grille

  • trunk/LMDZ.COMMON/libf/dyn3d_common/ran1.F

    r1299 r1300  
    2828      IX3=MOD(IA3*IX3+IC3,M3)
    2929      J=1+(97*IX3)/M3
    30       IF(J.GT.97.OR.J.LT.1)PAUSE
     30      IF(J.GT.97.OR.J.LT.1) stop 1
    3131      RAN1=R(J)
    3232      R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
Note: See TracChangeset for help on using the changeset viewer.