Ignore:
Timestamp:
Jan 17, 2014, 5:45:09 PM (10 years ago)
Author:
lguez
Message:

abort, dfloat and pause are not in the Fortran standard. Replaced
abort by abort_gcm and dfloat by dble. Note: I modified dyn3dpar files
that were identical to dyn3d modified files.

Location:
LMDZ5/trunk/libf/dyn3d
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/dyn3d/defrun.F

    r1907 r1930  
    364364        WRITE(tapeout,*)' La valeur de dzoomx passee par run.def est dif
    365365     *ferente de celle lue sur le fichier  start '
    366         CALL ABORT
     366        CALL ABORT_gcm("defrun", "", 1)
    367367       ENDIF
    368368
    … …  
    370370        WRITE(tapeout,*)' La valeur de dzoomy 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(6,*)' La valeur de taux passee par run.def est differente
    377377     *  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 tauy 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

  • LMDZ5/trunk/libf/dyn3d/dynetat0.F

    r1907 r1930  
    6767        write(lunout,*)'dynetat0: Pb d''ouverture du fichier start.nc'
    6868        write(lunout,*)' ierr = ', ierr
    69         CALL ABORT
     69        CALL ABORT_gcm("dynetat0", "", 1)
    7070      ENDIF
    7171
    … …  
    7474      IF (ierr .NE. NF_NOERR) THEN
    7575         write(lunout,*)"dynetat0: Le champ <controle> est absent"
    76          CALL abort
     76         CALL ABORT_gcm("dynetat0", "", 1)
    7777      ENDIF
    7878      ierr = nf90_get_var(nid, nvarid, tab_cntrl)
    7979      IF (ierr .NE. NF_NOERR) THEN
    8080         write(lunout,*)"dynetat0: Lecture echoue pour <controle>"
    81          CALL abort
     81         CALL ABORT_gcm("dynetat0", "", 1)
    8282      ENDIF
    8383
    … …  
    154154      IF (ierr .NE. NF_NOERR) THEN
    155155         write(lunout,*)"dynetat0: Le champ <rlonu> est absent"
    156          CALL abort
     156         CALL ABORT_gcm("dynetat0", "", 1)
    157157      ENDIF
    158158      ierr = nf90_get_var(nid, nvarid, rlonu)
    159159      IF (ierr .NE. NF_NOERR) THEN
    160160         write(lunout,*)"dynetat0: Lecture echouee pour <rlonu>"
    161          CALL abort
     161         CALL ABORT_gcm("dynetat0", "", 1)
    162162      ENDIF
    163163
    … …  
    165165      IF (ierr .NE. NF_NOERR) THEN
    166166         write(lunout,*)"dynetat0: Le champ <rlatu> est absent"
    167          CALL abort
     167         CALL ABORT_gcm("dynetat0", "", 1)
    168168      ENDIF
    169169      ierr = nf90_get_var(nid, nvarid, rlatu)
    170170      IF (ierr .NE. NF_NOERR) THEN
    171171         write(lunout,*)"dynetat0: Lecture echouee pour <rlatu>"
    172          CALL abort
     172         CALL ABORT_gcm("dynetat0", "", 1)
    173173      ENDIF
    174174
    … …  
    176176      IF (ierr .NE. NF_NOERR) THEN
    177177         write(lunout,*)"dynetat0: Le champ <rlonv> est absent"
    178          CALL abort
     178         CALL ABORT_gcm("dynetat0", "", 1)
    179179      ENDIF
    180180      ierr = nf90_get_var(nid, nvarid, rlonv)
    181181      IF (ierr .NE. NF_NOERR) THEN
    182182         write(lunout,*)"dynetat0: Lecture echouee pour <rlonv>"
    183          CALL abort
     183         CALL ABORT_gcm("dynetat0", "", 1)
    184184      ENDIF
    185185
    … …  
    187187      IF (ierr .NE. NF_NOERR) THEN
    188188         write(lunout,*)"dynetat0: Le champ <rlatv> est absent"
    189          CALL abort
     189         CALL ABORT_gcm("dynetat0", "", 1)
    190190      ENDIF
    191191      ierr = nf90_get_var(nid, nvarid, rlatv)
    192192      IF (ierr .NE. NF_NOERR) THEN
    193193         write(lunout,*)"dynetat0: Lecture echouee pour rlatv"
    194          CALL abort
     194         CALL ABORT_gcm("dynetat0", "", 1)
    195195      ENDIF
    196196
    … …  
    198198      IF (ierr .NE. NF_NOERR) THEN
    199199         write(lunout,*)"dynetat0: Le champ <cu> est absent"
    200          CALL abort
     200         CALL ABORT_gcm("dynetat0", "", 1)
    201201      ENDIF
    202202      ierr = nf90_get_var(nid, nvarid, cu)
    203203      IF (ierr .NE. NF_NOERR) THEN
    204204         write(lunout,*)"dynetat0: Lecture echouee pour <cu>"
    205          CALL abort
     205         CALL ABORT_gcm("dynetat0", "", 1)
    206206      ENDIF
    207207
    … …  
    209209      IF (ierr .NE. NF_NOERR) THEN
    210210         write(lunout,*)"dynetat0: Le champ <cv> est absent"
    211          CALL abort
     211         CALL ABORT_gcm("dynetat0", "", 1)
    212212      ENDIF
    213213      ierr = nf90_get_var(nid, nvarid, cv)
    214214      IF (ierr .NE. NF_NOERR) THEN
    215215         write(lunout,*)"dynetat0: Lecture echouee pour <cv>"
    216          CALL abort
     216         CALL ABORT_gcm("dynetat0", "", 1)
    217217      ENDIF
    218218
    … …  
    220220      IF (ierr .NE. NF_NOERR) THEN
    221221         write(lunout,*)"dynetat0: Le champ <aire> est absent"
    222          CALL abort
     222         CALL ABORT_gcm("dynetat0", "", 1)
    223223      ENDIF
    224224      ierr = nf90_get_var(nid, nvarid, aire)
    225225      IF (ierr .NE. NF_NOERR) THEN
    226226         write(lunout,*)"dynetat0: Lecture echouee pour <aire>"
    227          CALL abort
     227         CALL ABORT_gcm("dynetat0", "", 1)
    228228      ENDIF
    229229
    … …  
    231231      IF (ierr .NE. NF_NOERR) THEN
    232232         write(lunout,*)"dynetat0: Le champ <phisinit> est absent"
    233          CALL abort
     233         CALL ABORT_gcm("dynetat0", "", 1)
    234234      ENDIF
    235235      ierr = nf90_get_var(nid, nvarid, phis)
    236236      IF (ierr .NE. NF_NOERR) THEN
    237237         write(lunout,*)"dynetat0: Lecture echouee pour <phisinit>"
    238          CALL abort
     238         CALL ABORT_gcm("dynetat0", "", 1)
    239239      ENDIF
    240240
    … …  
    246246         IF (ierr .NE. NF_NOERR) THEN
    247247            write(lunout,*)"dynetat0: Le champ <Time> est absent"
    248             CALL abort
     248            CALL ABORT_gcm("dynetat0", "", 1)
    249249         ENDIF
    250250      ENDIF
    … …  
    252252      IF (ierr .NE. NF_NOERR) THEN
    253253         write(lunout,*)"dynetat0: Lecture echouee <temps>"
    254          CALL abort
     254         CALL ABORT_gcm("dynetat0", "", 1)
    255255      ENDIF
    256256
    … …  
    258258      IF (ierr .NE. NF_NOERR) THEN
    259259         write(lunout,*)"dynetat0: Le champ <ucov> est absent"
    260          CALL abort
     260         CALL ABORT_gcm("dynetat0", "", 1)
    261261      ENDIF
    262262      ierr = nf90_get_var(nid, nvarid, ucov)
    263263      IF (ierr .NE. NF_NOERR) THEN
    264264         write(lunout,*)"dynetat0: Lecture echouee pour <ucov>"
    265          CALL abort
     265         CALL ABORT_gcm("dynetat0", "", 1)
    266266      ENDIF
    267267 
    … …  
    269269      IF (ierr .NE. NF_NOERR) THEN
    270270         write(lunout,*)"dynetat0: Le champ <vcov> est absent"
    271          CALL abort
     271         CALL ABORT_gcm("dynetat0", "", 1)
    272272      ENDIF
    273273      ierr = nf90_get_var(nid, nvarid, vcov)
    274274      IF (ierr .NE. NF_NOERR) THEN
    275275         write(lunout,*)"dynetat0: Lecture echouee pour <vcov>"
    276          CALL abort
     276         CALL ABORT_gcm("dynetat0", "", 1)
    277277      ENDIF
    278278
    … …  
    280280      IF (ierr .NE. NF_NOERR) THEN
    281281         write(lunout,*)"dynetat0: Le champ <teta> est absent"
    282          CALL abort
     282         CALL ABORT_gcm("dynetat0", "", 1)
    283283      ENDIF
    284284      ierr = nf90_get_var(nid, nvarid, teta)
    285285      IF (ierr .NE. NF_NOERR) THEN
    286286         write(lunout,*)"dynetat0: Lecture echouee pour <teta>"
    287          CALL abort
     287         CALL ABORT_gcm("dynetat0", "", 1)
    288288      ENDIF
    289289
    … …  
    301301          IF (ierr .NE. NF_NOERR) THEN
    302302            write(lunout,*)"dynetat0: Lecture echouee pour "//tname(iq)
    303             CALL abort
     303            CALL ABORT_gcm("dynetat0", "", 1)
    304304          ENDIF
    305305        ENDIF
    … …  
    310310      IF (ierr .NE. NF_NOERR) THEN
    311311         write(lunout,*)"dynetat0: Le champ <masse> est absent"
    312          CALL abort
     312         CALL ABORT_gcm("dynetat0", "", 1)
    313313      ENDIF
    314314      ierr = nf90_get_var(nid, nvarid, masse)
    315315      IF (ierr .NE. NF_NOERR) THEN
    316316         write(lunout,*)"dynetat0: Lecture echouee pour <masse>"
    317          CALL abort
     317         CALL ABORT_gcm("dynetat0", "", 1)
    318318      ENDIF
    319319
    … …  
    321321      IF (ierr .NE. NF_NOERR) THEN
    322322         write(lunout,*)"dynetat0: Le champ <ps> est absent"
    323          CALL abort
     323         CALL ABORT_gcm("dynetat0", "", 1)
    324324      ENDIF
    325325      ierr = nf90_get_var(nid, nvarid, ps)
    326326      IF (ierr .NE. NF_NOERR) THEN
    327327         write(lunout,*)"dynetat0: Lecture echouee pour <ps>"
    328          CALL abort
     328         CALL ABORT_gcm("dynetat0", "", 1)
    329329      ENDIF
    330330

  • LMDZ5/trunk/libf/dyn3d/dynredem.F

    r1907 r1930  
    133133     &                  //trim(fichnom)
    134134         write(lunout,*)' ierr = ', ierr
    135          CALL ABORT
     135         CALL ABORT_GCM("DYNREDEM0", "", 1)
    136136      ENDIF
    137137c
    … …  
    512512      IF (ierr .NE. NF_NOERR) THEN
    513513         write(lunout,*)"dynredem1: Pb. d ouverture "//trim(fichnom)
    514          CALL abort
     514         call abort_gcm("dynredem1", "", 1)
    515515      ENDIF
    516516

  • LMDZ5/trunk/libf/dyn3d/fxhyp.F

    r1907 r1930  
    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

  • LMDZ5/trunk/libf/dyn3d/fyhyp.F

    r1907 r1930  
    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

  • LMDZ5/trunk/libf/dyn3d/gcm.F

    r1907 r1930  
    327327          start_time = starttime
    328328        ELSE
    329           WRITE(lunout,*)'Je m''arrete'
    330           CALL abort
     329          call abort_gcm("gcm", "'Je m''arrete'", 1)
    331330        ENDIF
    332331      ENDIF

  • LMDZ5/trunk/libf/dyn3d/grid_atob.F

    r1907 r1930  
    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

  • LMDZ5/trunk/libf/dyn3d/grid_noro.F

    r1907 r1930  
    111111      IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
    112112         PRINT*, 'imar or jmar too big', imar, jmar
    113          CALL ABORT
     113         CALL ABORT_GCM("GRID_NORO", "", 1)
    114114      ENDIF
    115115
    116116      IF(imdep.ne.iusn.or.jmdep.ne.jusn)then
    117117         print *,' imdep or jmdep bad dimensions:',imdep,jmdep
    118          call abort
     118         call abort_gcm("grid_noro", "", 1)
    119119      ENDIF
    120120
    121121      IF(imar+1.ne.iim+1.or.jmar.ne.jjm+1)THEN
    122122        print *,' imar or jmar bad dimensions:',imar,jmar
    123         call abort
     123        call abort_gcm("grid_noro", "", 1)
    124124      ENDIF
    125125

  • LMDZ5/trunk/libf/dyn3d/juldate.F

    r1907 r1930  
    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.

  • LMDZ5/trunk/libf/dyn3d/ran1.F

    r1907 r1930  
    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.