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
Files:
4 deleted
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d/gcm.F

    r1107 r1300  
    348348          start_time = starttime
    349349        ELSE
    350           WRITE(lunout,*)'Je m''arrete'
    351           CALL abort
     350          call abort_gcm("gcm", "'Je m''arrete'", 1)
    352351        ENDIF
    353352      ENDIF

  • trunk/LMDZ.COMMON/libf/dyn3d/guide_mod.F90

    r979 r1300  
    8787    CALL getpar('guide_teta',.false.,guide_teta,'guidage de T par Teta')
    8888
    89     CALL getpar('guide_add',.false.,guide_add,'for�age constant?')
     89    CALL getpar('guide_add',.false.,guide_add,'for�age constant?')
    9090    CALL getpar('guide_zon',.false.,guide_zon,'guidage moy zonale')
    9191
    … …  
    104104    CALL getpar('guide_BL',.true.,guide_BL,'guidage dans C.Lim')
    105105   
    106 ! Sauvegarde du for�age
     106! Sauvegarde du for�age
    107107    CALL getpar('guide_sav',.false.,guide_sav,'sauvegarde guidage')
    108108    CALL getpar('iguide_sav',4,iguide_sav,'freq. sauvegarde guidage')
    … …  
    143143    ncidpl=-99
    144144    if (guide_modele) then
    145        if (ncidpl.eq.-99) rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
     145       if (ncidpl.eq.-99) then
     146          rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
     147          if (rcod.NE.NF_NOERR) THEN
     148             print *,'Guide: probleme -> pas de fichier apbp.nc'
     149             CALL abort_gcm(modname,abort_message,1)
     150          endif
     151       endif
    146152    else
    147153         if (guide_u) then
    148            if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
     154           if (ncidpl.eq.-99) then
     155               rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
     156               if (rcod.NE.NF_NOERR) THEN
     157                  print *,'Guide: probleme -> pas de fichier u.nc'
     158                  CALL abort_gcm(modname,abort_message,1)
     159               endif
     160           endif
    149161         elseif (guide_v) then
    150            if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
     162           if (ncidpl.eq.-99) then
     163               rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
     164               if (rcod.NE.NF_NOERR) THEN
     165                  print *,'Guide: probleme -> pas de fichier v.nc'
     166                  CALL abort_gcm(modname,abort_message,1)
     167               endif
     168           endif
    151169         elseif (guide_T) then
    152            if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
     170           if (ncidpl.eq.-99) then
     171               rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
     172               if (rcod.NE.NF_NOERR) THEN
     173                  print *,'Guide: probleme -> pas de fichier T.nc'
     174                  CALL abort_gcm(modname,abort_message,1)
     175               endif
     176           endif
    153177         elseif (guide_Q) then
    154            if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
     178           if (ncidpl.eq.-99) then
     179               rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
     180               if (rcod.NE.NF_NOERR) THEN
     181                  print *,'Guide: probleme -> pas de fichier hur.nc'
     182                  CALL abort_gcm(modname,abort_message,1)
     183               endif
     184           endif
    155185         endif
    156186    endif
    … …  
    9901020    INTEGER               :: status,rcode
    9911021
     1022    CHARACTER (len = 80)   :: abort_message
     1023    CHARACTER (len = 20)   :: modname = 'guide_read'
    9921024! -----------------------------------------------------------------
    9931025! Premier appel: initialisation de la lecture des fichiers
    … …  
    9981030! Niveaux de pression si non constants
    9991031         if (guide_modele) then
    1000              print *,'Lecture du guidage sur niveaux mod�le'
     1032             print *,'Lecture du guidage sur niveaux modele'
    10011033             rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
     1034             IF (rcode.NE.NF_NOERR) THEN
     1035              print *,'Guide: probleme -> pas de fichier apbp.nc'
     1036              CALL abort_gcm(modname,abort_message,1)
     1037             ENDIF
    10021038             rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
     1039             IF (rcode.NE.NF_NOERR) THEN
     1040              print *,'Guide: probleme -> pas de variable AP, fichier apbp.nc'
     1041              CALL abort_gcm(modname,abort_message,1)
     1042             ENDIF
    10031043             rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
     1044             IF (rcode.NE.NF_NOERR) THEN
     1045              print *,'Guide: probleme -> pas de variable BP, fichier apbp.nc'
     1046              CALL abort_gcm(modname,abort_message,1)
     1047             ENDIF
    10041048             print*,'ncidpl,varidap',ncidpl,varidap
    10051049         endif
    … …  
    10071051         if (guide_u) then
    10081052             rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
     1053             IF (rcode.NE.NF_NOERR) THEN
     1054              print *,'Guide: probleme -> pas de fichier u.nc'
     1055              CALL abort_gcm(modname,abort_message,1)
     1056             ENDIF
    10091057             rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
     1058             IF (rcode.NE.NF_NOERR) THEN
     1059              print *,'Guide: probleme -> pas de variable UWND, fichier u.nc'
     1060              CALL abort_gcm(modname,abort_message,1)
     1061             ENDIF
    10101062             print*,'ncidu,varidu',ncidu,varidu
    10111063             if (ncidpl.eq.-99) ncidpl=ncidu
    … …  
    10141066         if (guide_v) then
    10151067             rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
     1068             IF (rcode.NE.NF_NOERR) THEN
     1069              print *,'Guide: probleme -> pas de fichier v.nc'
     1070              CALL abort_gcm(modname,abort_message,1)
     1071             ENDIF
    10161072             rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
     1073             IF (rcode.NE.NF_NOERR) THEN
     1074              print *,'Guide: probleme -> pas de variable VWND, fichier v.nc'
     1075              CALL abort_gcm(modname,abort_message,1)
     1076             ENDIF
    10171077             print*,'ncidv,varidv',ncidv,varidv
    10181078             if (ncidpl.eq.-99) ncidpl=ncidv
    … …  
    10211081         if (guide_T) then
    10221082             rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
     1083             IF (rcode.NE.NF_NOERR) THEN
     1084              print *,'Guide: probleme -> pas de fichier T.nc'
     1085              CALL abort_gcm(modname,abort_message,1)
     1086             ENDIF
    10231087             rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
     1088             IF (rcode.NE.NF_NOERR) THEN
     1089              print *,'Guide: probleme -> pas de variable AIR, fichier T.nc'
     1090              CALL abort_gcm(modname,abort_message,1)
     1091             ENDIF
    10241092             print*,'ncidT,varidT',ncidt,varidt
    10251093             if (ncidpl.eq.-99) ncidpl=ncidt
    … …  
    10281096         if (guide_Q) then
    10291097             rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
     1098             IF (rcode.NE.NF_NOERR) THEN
     1099              print *,'Guide: probleme -> pas de fichier hur.nc'
     1100              CALL abort_gcm(modname,abort_message,1)
     1101             ENDIF
    10301102             rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
     1103             IF (rcode.NE.NF_NOERR) THEN
     1104              print *,'Guide: probleme -> pas de variable RH, fichier hur.nc'
     1105              CALL abort_gcm(modname,abort_message,1)
     1106             ENDIF
    10311107             print*,'ncidQ,varidQ',ncidQ,varidQ
    10321108             if (ncidpl.eq.-99) ncidpl=ncidQ
    … …  
    10351111         if ((guide_P).OR.(guide_modele)) then
    10361112             rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
     1113             IF (rcode.NE.NF_NOERR) THEN
     1114              print *,'Guide: probleme -> pas de fichier ps.nc'
     1115              CALL abort_gcm(modname,abort_message,1)
     1116             ENDIF
    10371117             rcode = nf90_inq_varid(ncidps, 'SP', varidps)
     1118             IF (rcode.NE.NF_NOERR) THEN
     1119              print *,'Guide: probleme -> pas de variable SP, fichier ps.nc'
     1120              CALL abort_gcm(modname,abort_message,1)
     1121             ENDIF
    10381122             print*,'ncidps,varidps',ncidps,varidps
    10391123         endif
    … …  
    11731257    INTEGER               :: i
    11741258
     1259    CHARACTER (len = 80)   :: abort_message
     1260    CHARACTER (len = 20)   :: modname = 'guide_read2D'
    11751261! -----------------------------------------------------------------
    11761262! Premier appel: initialisation de la lecture des fichiers
    … …  
    11811267! Niveaux de pression si non constants
    11821268         if (guide_modele) then
    1183              print *,'Lecture du guidage sur niveaux mod�le'
     1269             print *,'Lecture du guidage sur niveaux modele'
    11841270             rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
     1271             IF (rcode.NE.NF_NOERR) THEN
     1272              print *,'Guide: probleme -> pas de fichier apbp.nc'
     1273              CALL abort_gcm(modname,abort_message,1)
     1274             ENDIF
    11851275             rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
     1276             IF (rcode.NE.NF_NOERR) THEN
     1277              print *,'Guide: probleme -> pas de variable AP, fichier apbp.nc'
     1278              CALL abort_gcm(modname,abort_message,1)
     1279             ENDIF
    11861280             rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
     1281             IF (rcode.NE.NF_NOERR) THEN
     1282              print *,'Guide: probleme -> pas de variable BP, fichier apbp.nc'
     1283              CALL abort_gcm(modname,abort_message,1)
     1284             ENDIF
    11871285             print*,'ncidpl,varidap',ncidpl,varidap
    11881286         endif
    … …  
    11901288         if (guide_u) then
    11911289             rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
     1290             IF (rcode.NE.NF_NOERR) THEN
     1291              print *,'Guide: probleme -> pas de fichier u.nc'
     1292              CALL abort_gcm(modname,abort_message,1)
     1293             ENDIF
    11921294             rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
     1295             IF (rcode.NE.NF_NOERR) THEN
     1296              print *,'Guide: probleme -> pas de variable UWND, fichier u.nc'
     1297              CALL abort_gcm(modname,abort_message,1)
     1298             ENDIF
    11931299             print*,'ncidu,varidu',ncidu,varidu
    11941300             if (ncidpl.eq.-99) ncidpl=ncidu
    … …  
    11971303         if (guide_v) then
    11981304             rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
     1305             IF (rcode.NE.NF_NOERR) THEN
     1306              print *,'Guide: probleme -> pas de fichier v.nc'
     1307              CALL abort_gcm(modname,abort_message,1)
     1308             ENDIF
    11991309             rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
     1310             IF (rcode.NE.NF_NOERR) THEN
     1311              print *,'Guide: probleme -> pas de variable VWND, fichier v.nc'
     1312              CALL abort_gcm(modname,abort_message,1)
     1313             ENDIF
    12001314             print*,'ncidv,varidv',ncidv,varidv
    12011315             if (ncidpl.eq.-99) ncidpl=ncidv
    … …  
    12041318         if (guide_T) then
    12051319             rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
     1320             IF (rcode.NE.NF_NOERR) THEN
     1321              print *,'Guide: probleme -> pas de fichier T.nc'
     1322              CALL abort_gcm(modname,abort_message,1)
     1323             ENDIF
    12061324             rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
     1325             IF (rcode.NE.NF_NOERR) THEN
     1326              print *,'Guide: probleme -> pas de variable AIR, fichier T.nc'
     1327              CALL abort_gcm(modname,abort_message,1)
     1328             ENDIF
    12071329             print*,'ncidT,varidT',ncidt,varidt
    12081330             if (ncidpl.eq.-99) ncidpl=ncidt
    … …  
    12111333         if (guide_Q) then
    12121334             rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
     1335             IF (rcode.NE.NF_NOERR) THEN
     1336              print *,'Guide: probleme -> pas de fichier hur.nc'
     1337              CALL abort_gcm(modname,abort_message,1)
     1338             ENDIF
    12131339             rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
     1340             IF (rcode.NE.NF_NOERR) THEN
     1341              print *,'Guide: probleme -> pas de variable RH, fichier hur.nc'
     1342              CALL abort_gcm(modname,abort_message,1)
     1343             ENDIF
    12141344             print*,'ncidQ,varidQ',ncidQ,varidQ
    12151345             if (ncidpl.eq.-99) ncidpl=ncidQ
    … …  
    12181348         if ((guide_P).OR.(guide_modele)) then
    12191349             rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
     1350             IF (rcode.NE.NF_NOERR) THEN
     1351              print *,'Guide: probleme -> pas de fichier ps.nc'
     1352              CALL abort_gcm(modname,abort_message,1)
     1353             ENDIF
    12201354             rcode = nf90_inq_varid(ncidps, 'SP', varidps)
     1355             IF (rcode.NE.NF_NOERR) THEN
     1356              print *,'Guide: probleme -> pas de variable SP, fichier ps.nc'
     1357              CALL abort_gcm(modname,abort_message,1)
     1358             ENDIF
    12211359             print*,'ncidps,varidps',ncidps,varidps
    12221360         endif
    … …  
    14271565#endif
    14281566! --------------------------------------------------------------------
    1429 ! Cr�ation des variables sauvegard�es
     1567! Cr�ation des variables sauvegard�es
    14301568! --------------------------------------------------------------------
    14311569        ierr = NF_REDEF(nid)
    … …  
    15521690!===========================================================================
    15531691END MODULE guide_mod
     1692
Note: See TracChangeset for help on using the changeset viewer.