Ignore:
Timestamp:
Jul 22, 2024, 9:29:09 PM (2 months ago)
Author:
abarral
Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

Location:
LMDZ6/branches/Amaury_dev/libf/dyn3d
Files:
34 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/abort_gcm.F

    r5082 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44c
    55c

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/addfi.F

    r2598 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE addfi(pdt, leapf, forward,
    55     S          pucov, pvcov, pteta, pq   , pps ,

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F

    r5086 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE advect(ucov,vcov,teta,w,massebx,masseby,du,dv,dteta)
    55

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/advtrac.F90

    r5082 r5099  
    55SUBROUTINE advtrac(pbaru, pbarv, p, masse,q,iapptrac,teta, flxw, pk)
    66   !     Auteur :  F. Hourdin
    7    !
     7
    88   !     Modif. P. Le Van     (20/12/97)
    99   !            F. Codron     (10/99)
    1010   !            D. Le Croller (07/2001)
    1111   !            M.A Filiberti (04/2002)
    12    !
     12
    1313   USE infotrac,     ONLY: nqtot, tracers, isoCheck
    1414   USE control_mod,  ONLY: iapp_tracvl, day_step
    … …  
    1616
    1717   IMPLICIT NONE
    18    !
     18
    1919   include "dimensions.h"
    2020   include "paramet.h"
    … …  
    9696
    9797   !   ..  Modif P.Le Van  ( 20/12/97 )  ....
    98    !
     98
    9999   !   traitement des flux de masse avant advection.
    100100   !       1. calcul de w

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/bilan_dyn.F

    r5093 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE bilan_dyn (ntrac,dt_app,dt_cum,
    55     s  ps,masse,pk,flux_u,flux_v,teta,phi,ucov,vcov,trac)

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/caladvtrac.F

    r5082 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44c
    55c

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/caldyn.F

    r5082 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE caldyn
    55     $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
    … …  
    1212
    1313!=======================================================================
    14 !
     14
    1515!  Auteur :  P. Le Van
    16 !
     16
    1717!   Objet:
    1818!   ------
    19 !
     19
    2020!   Calcul des tendances dynamiques.
    21 !
     21
    2222! Modif 04/93 F.Forget
    2323!=======================================================================

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/covnat.F

    r1907 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
    55      IMPLICIT NONE

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dissip.F

    r2597 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE dissip( vcov,ucov,teta,p, dv,du,dh )
    55c

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dteta1.F

    r5086 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE dteta1 ( teta, pbaru, pbarv, dteta)
    55      IMPLICIT NONE

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv1.F

    r5086 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE dudv1 ( vorpot, pbaru, pbarv, du, dv )
    55      IMPLICIT NONE

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F

    r5086 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE dudv2 ( teta, pkf, bern, du, dv  )
    55

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dynetat0.F90

    r5091 r5099  
    11SUBROUTINE dynetat0(fichnom,vcov,ucov,teta,q,masse,ps,phis,time)
    2 !
     2
    33!-------------------------------------------------------------------------------
    44! Authors: P. Le Van , L.Fairhead
    … …  
    99  USE strings_mod, ONLY: maxlen, msg, strStack, real2str, int2str
    1010  USE netcdf,      ONLY: NF90_OPEN,  NF90_NOWRITE, NF90_INQ_VARID, &
    11                          NF90_CLOSE, NF90_GET_VAR, NF90_NoErr
     11                         NF90_CLOSE, nf90_get_var, nf90_noerr
    1212  USE readTracFiles_mod, ONLY: new2oldH2O, newHNO3, oldHNO3, getKey
    1313  USE control_mod, ONLY: planet_type
    … …  
    8282  pa         = tab_cntrl(idecal+13)
    8383  preff      = tab_cntrl(idecal+14)
    84 !
     84
    8585  clon       = tab_cntrl(idecal+15)
    8686  clat       = tab_cntrl(idecal+16)
    8787  grossismx  = tab_cntrl(idecal+17)
    8888  grossismy  = tab_cntrl(idecal+18)
    89 !
     89
    9090  IF ( tab_cntrl(idecal+19)==1. )  THEN
    9191    fxyhypb  = .TRUE.
    … …  
    116116  CALL get_var2("aire" ,aire)
    117117  var="temps"
    118   IF(NF90_INQ_VARID(fID,var,vID)/=NF90_NoErr) THEN
     118  IF(NF90_INQ_VARID(fID,var,vID)/=nf90_noerr) THEN
    119119    CALL msg('missing field <temps> ; trying with <Time>', modname)
    120120    var="Time"
    121121    CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
    122122  END IF
    123   CALL err(NF90_GET_VAR(fID,vID,time),"get",var)
     123  CALL err(nf90_get_var(fID,vID,time),"get",var)
    124124  CALL get_var2("phisinit",phis)
    125125  CALL get_var3("ucov",ucov)
    … …  
    132132  ll=.FALSE.
    133133#ifdef REPROBUS
    134   ll = NF90_INQ_VARID(fID, 'HNO3tot', vID) /= NF90_NoErr                                 !--- DETECT OLD REPRO start.nc FILE
     134  ll = NF90_INQ_VARID(fID, 'HNO3tot', vID) /= nf90_noerr                                 !--- DETECT OLD REPRO start.nc FILE
    135135#endif
    136136  DO iq=1,nqtot
    … …  
    145145    END IF
    146146    !--------------------------------------------------------------------------------------------------------------------------
    147     IF(NF90_INQ_VARID(fID, var, vID) == NF90_NoErr .AND. .NOT.lSkip) THEN                !=== REGULAR CASE: AVAILABLE VARIABLE
    148       CALL err(NF90_GET_VAR(fID,vID,q(:,:,:,iq)),"get",var)
    149     !--------------------------------------------------------------------------------------------------------------------------
    150     ELSE IF(NF90_INQ_VARID(fID, oldVar, vID) == NF90_NoErr) THEN                         !=== TRY WITH ALTERNATE NAME
     147    IF(NF90_INQ_VARID(fID, var, vID) == nf90_noerr .AND. .NOT.lSkip) THEN                !=== REGULAR CASE: AVAILABLE VARIABLE
     148      CALL err(nf90_get_var(fID,vID,q(:,:,:,iq)),"get",var)
     149    !--------------------------------------------------------------------------------------------------------------------------
     150    ELSE IF(NF90_INQ_VARID(fID, oldVar, vID) == nf90_noerr) THEN                         !=== TRY WITH ALTERNATE NAME
    151151      CALL msg('Tracer <'//TRIM(var)//'> is missing => initialized to <'//TRIM(oldVar)//'>', modname)
    152       CALL err(NF90_GET_VAR(fID,vID,q(:,:,:,iq)),"get",oldVar)
     152      CALL err(nf90_get_var(fID,vID,q(:,:,:,iq)),"get",oldVar)
    153153    !--------------------------------------------------------------------------------------------------------------------------
    154154    ELSE IF(tracers(iq)%iso_iGroup == iH2O .AND. niso > 0) THEN                          !=== WATER ISOTOPES
    … …  
    209209  REAL,             INTENT(OUT) :: v(:)
    210210  CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
    211   CALL err(NF90_GET_VAR(fID,vID,v),"get",var)
     211  CALL err(nf90_get_var(fID,vID,v),"get",var)
    212212END SUBROUTINE get_var1
    213213
    … …  
    217217  REAL,             INTENT(OUT) :: v(:,:)
    218218  CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
    219   CALL err(NF90_GET_VAR(fID,vID,v),"get",var)
     219  CALL err(nf90_get_var(fID,vID,v),"get",var)
    220220END SUBROUTINE get_var2
    221221
    … …  
    225225  REAL,             INTENT(OUT) :: v(:,:,:)
    226226  CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var)
    227   CALL err(NF90_GET_VAR(fID,vID,v),"get",var)
     227  CALL err(nf90_get_var(fID,vID,v),"get",var)
    228228END SUBROUTINE get_var3
    229229
    … …  
    233233  CHARACTER(LEN=*), INTENT(IN) :: typ    !--- TYPE OF OPERATION
    234234  CHARACTER(LEN=*), INTENT(IN) :: nam    !--- FIELD/FILE NAME
    235   IF(ierr==NF90_NoERR) RETURN
     235  IF(ierr==nf90_noerr) RETURN
    236236  SELECT CASE(typ)
    237237    CASE('inq');   mesg="Field <"//TRIM(nam)//"> is missing"

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dynredem.F90

    r5088 r5099  
    11SUBROUTINE dynredem0(fichnom,iday_end,phis)
    2 !
     2
    33!-------------------------------------------------------------------------------
    44! Write the NetCDF restart file (initialization).
    … …  
    8282  tab_cntrl(22) = grossismx
    8383  tab_cntrl(23) = grossismy
    84 !
     84
    8585  IF ( fxyhypb )   THEN
    8686    tab_cntrl(24) = 1.
    … …  
    155155
    156156END SUBROUTINE dynredem0
    157 !
    158 !-------------------------------------------------------------------------------
    159 
    160 
    161 !-------------------------------------------------------------------------------
    162 !
     157
     158!-------------------------------------------------------------------------------
     159
     160
     161!-------------------------------------------------------------------------------
     162
    163163SUBROUTINE dynredem1(fichnom,time,vcov,ucov,teta,q,masse,ps)
    164 !
     164
    165165!-------------------------------------------------------------------------------
    166166! Purpose: Write the NetCDF restart file (append).
    … …  
    169169  USE infotrac, ONLY: nqtot, tracers, type_trac
    170170  USE control_mod
    171   USE netcdf,   ONLY: NF90_OPEN,  NF90_NOWRITE, NF90_GET_VAR, NF90_INQ_VARID,  &
    172                       NF90_CLOSE, NF90_WRITE,   NF90_PUT_VAR, NF90_NoErr
     171  USE netcdf,   ONLY: NF90_OPEN,  NF90_NOWRITE, nf90_get_var, NF90_INQ_VARID,  &
     172                      NF90_CLOSE, NF90_WRITE,   NF90_PUT_VAR, nf90_noerr
    173173  USE dynredem_mod, ONLY: dynredem_write_u, dynredem_write_v, dynredem_read_u, &
    174174                          err, modname, fil, msg
    … …  
    214214  var="controle"
    215215  CALL err(NF90_INQ_VARID(nid,var,vID),"inq",var)
    216   CALL err(NF90_GET_VAR(nid,vID,tab_cntrl),"get",var)
     216  CALL err(nf90_get_var(nid,vID,tab_cntrl),"get",var)
    217217  tab_cntrl(31)=DBLE(itau_dyn + itaufin)
    218218  CALL err(NF90_INQ_VARID(nid,var,vID),"inq",var)
    … …  
    236236      fil="start_trac.nc"
    237237      ierr=NF90_INQ_VARID(nid_trac,var,vID_trac)
    238       dum='inq'; IF(ierr==NF90_NoErr) dum='fnd'
     238      dum='inq'; IF(ierr==nf90_noerr) dum='fnd'
    239239      WRITE(lunout,*)msg(dum,var)
    240240
    241241
    242       IF(ierr==NF90_NoErr) CALL dynredem_read_u(nid_trac,var,q(:,:,:,iq),llm)
     242      IF(ierr==nf90_noerr) CALL dynredem_read_u(nid_trac,var,q(:,:,:,iq),llm)
    243243    END IF
    244244    fil=fichnom

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dynredem_mod.F90

    r5090 r5099  
    1717
    1818!===============================================================================
    19 !
     19
    2020SUBROUTINE dynredem_write_u(ncid,id,var,ll)
    21 !
     21
    2222!===============================================================================
    2323! Arguments:
    … …  
    3535 
    3636END SUBROUTINE dynredem_write_u
    37 !
    38 !===============================================================================
    39 
    40 
    41 !===============================================================================
    42 !
     37
     38!===============================================================================
     39
     40
     41!===============================================================================
     42
    4343SUBROUTINE dynredem_write_v(ncid,id,var,ll)
    44 !
     44
    4545!===============================================================================
    4646! Arguments:
    … …  
    5858 
    5959END SUBROUTINE dynredem_write_v
    60 !
    61 !===============================================================================
    62 
    63 
    64 !===============================================================================
    65 !
     60
     61!===============================================================================
     62
     63
     64!===============================================================================
     65
    6666SUBROUTINE dynredem_read_u(ncid,id,var,ll)
    67 !
     67
    6868!===============================================================================
    6969! Arguments:
    … …  
    7878  start(:)=[1,1,1,1]; count(:)=[iip1,jjp1,ll,1]
    7979  CALL err(NF90_INQ_VARID(ncid,id,nvarid),"inq",id)
    80   CALL err(NF90_GET_VAR(ncid,nvarid,var,start,count),"get",id)
     80  CALL err(nf90_get_var(ncid,nvarid,var,start,count),"get",id)
    8181 
    8282END SUBROUTINE dynredem_read_u   
    83 !
    84 !===============================================================================
    85 
    86 
    87 !===============================================================================
    88 !
     83
     84!===============================================================================
     85
     86
     87!===============================================================================
     88
    8989SUBROUTINE cre_var(ncid,var,title,did,units)
    90 !
     90
    9191!===============================================================================
    9292! Arguments:
    … …  
    101101
    102102END SUBROUTINE cre_var
    103 !
    104 !===============================================================================
    105 
    106 
    107 !===============================================================================
    108 !
     103
     104!===============================================================================
     105
     106
     107!===============================================================================
     108
    109109SUBROUTINE put_var1(ncid,var,title,did,v,units)
    110 !
     110
    111111!===============================================================================
    112112! Arguments:
    … …  
    124124
    125125END SUBROUTINE put_var1
    126 !
    127 !===============================================================================
    128 
    129 
    130 !===============================================================================
    131 !
     126
     127!===============================================================================
     128
     129
     130!===============================================================================
     131
    132132SUBROUTINE put_var2(ncid,var,title,did,v,units)
    133 !
     133
    134134!===============================================================================
    135135! Arguments:
    … …  
    147147
    148148END SUBROUTINE put_var2
    149 !
    150 !===============================================================================
    151 
    152 
    153 !===============================================================================
    154 !
     149
     150!===============================================================================
     151
     152
     153!===============================================================================
     154
    155155FUNCTION msg(typ,nam)
    156 !
     156
    157157!===============================================================================
    158158! Arguments:
    … …  
    172172
    173173END FUNCTION msg
    174 !
    175 !===============================================================================
    176 
    177 
    178 !===============================================================================
    179 !
     174
     175!===============================================================================
     176
     177
     178!===============================================================================
     179
    180180SUBROUTINE err(ierr,typ,nam)
    181 !
     181
    182182!===============================================================================
    183183! Arguments:
    … …  
    186186  CHARACTER(LEN=*), OPTIONAL, INTENT(IN) :: nam    !--- FIELD NAME
    187187!===============================================================================
    188   IF(ierr==NF90_NoERR) RETURN
     188  IF(ierr==nf90_noerr) RETURN
    189189  IF(.NOT.PRESENT(typ)) THEN
    190190    CALL ABORT_gcm(modname,NF90_STRERROR(ierr),ierr)
    … …  
    194194
    195195END SUBROUTINE err
    196 !
     196
    197197!===============================================================================
    198198

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/fluxstokenc.F

    r5082 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE fluxstokenc(pbaru,pbarv,masse,teta,phi,phis,
    55     . time_step,itau )

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/friction.F

    r5082 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44c=======================================================================
    55      SUBROUTINE friction(ucov,vcov,pdt)
    … …  
    1616
    1717!=======================================================================
    18 !
     18
    1919!   Friction for the Newtonian case:
    2020!   --------------------------------

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/gcm.F90

    r5090 r5099  
    1 !
     1
    22! $Id$
    3 !
    4 !
    5 !
     3
     4
     5
    66PROGRAM gcm
    77
    … …  
    4444
    4545  !=======================================================================
    46   !
     46
    4747  !   Auteur:  P. Le Van /L. Fairhead/F.Hourdin
    4848  !   -------
    49   !
     49
    5050  !   Objet:
    5151  !   ------
    52   !
     52
    5353  !   GCM LMD nouvelle grille
    54   !
     54
    5555  !=======================================================================
    56   !
     56
    5757  !  ... Dans inigeom , nouveaux calculs pour les elongations  cu , cv
    5858  !      et possibilite d'appeler une fonction f(y)  a derivee tangente
    … …  
    6060  !  ... Possibilite de choisir le schema pour l'advection de
    6161  !        q  , en modifiant iadv dans traceur.def  (MAF,10/02) .
    62   !
     62
    6363  !      Pour Van-Leer + Vapeur d'eau saturee, iadv(1)=4. (F.Codron,10/99)
    6464  !      Pour Van-Leer iadv=10
    65   !
     65
    6666  !-----------------------------------------------------------------------
    6767  !   Declarations:
    … …  
    140140  !  lecture des fichiers gcm.def ou run.def
    141141  !  ---------------------------------------
    142   !
     142
    143143  CALL conf_gcm( 99, .TRUE.)
    144144
    … …  
    191191#endif
    192192  !-----------------------------------------------------------------------
    193   !
    194   !
     193
     194
    195195  !------------------------------------
    196196  !   Initialisation partie parallele
    197197  !------------------------------------
    198198
    199   !
    200   !
     199
    201200  !-----------------------------------------------------------------------
    202201  !   Initialisation des traceurs
    … …  
    269268  ENDIF
    270269
    271   !
    272270  ! on remet le calendrier \`a zero si demande
    273   !
     271
    274272  IF (start_time /= starttime) then
    275273     WRITE(lunout,*)' GCM: Attention l''heure de depart lue dans le' &
    … …  
    370368     CALL inifilr
    371369  endif ! of if (iflag_phys.eq.1)
    372   !
     370
    373371  !-----------------------------------------------------------------------
    374372  !   Initialisation de la dissipation :
    … …  
    453451  istphy=istdyn/iphysiq     
    454452
    455 
    456   !
    457453  !-----------------------------------------------------------------------
    458454  !   Integration temporelle du modele :

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/getparam.F90

    r2094 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44MODULE getparam
    55#ifdef CPP_IOIPSL
    … …  
    1919CONTAINS
    2020  SUBROUTINE ini_getparam(fichier)
    21     !
     21
    2222    IMPLICIT NONE
    23     !
     23
    2424    CHARACTER*(*) :: fichier
    2525    open(out_eff,file=fichier,status='unknown',form='formatted')
    … …  
    2727
    2828  SUBROUTINE fin_getparam
    29     !
     29
    3030    IMPLICIT NONE
    31     !
     31
    3232    close(out_eff)
    3333
    … …  
    3535
    3636  SUBROUTINE getparamr(TARGET,def_val,ret_val,comment)
    37     !
     37
    3838    IMPLICIT NONE
    39     !
     39
    4040    !   Get a real scalar. We first check if we find it
    4141    !   in the database and if not we get it from the run.def
    42     !
     42
    4343    !   getinr1d and getinr2d are written on the same pattern
    44     !
     44
    4545    CHARACTER*(*) :: TARGET
    4646    REAL :: def_val
    … …  
    5858
    5959  SUBROUTINE getparami(TARGET,def_val,ret_val,comment)
    60     !
     60
    6161    IMPLICIT NONE
    62     !
     62
    6363    !   Get a real scalar. We first check if we find it
    6464    !   in the database and if not we get it from the run.def
    65     !
     65
    6666    !   getinr1d and getinr2d are written on the same pattern
    67     !
     67
    6868    CHARACTER*(*) :: TARGET
    6969    INTEGER :: def_val
    … …  
    8282
    8383  SUBROUTINE getparaml(TARGET,def_val,ret_val,comment)
    84     !
     84
    8585    IMPLICIT NONE
    86     !
     86
    8787    !   Get a real scalar. We first check if we find it
    8888    !   in the database and if not we get it from the run.def
    89     !
     89
    9090    !   getinr1d and getinr2d are written on the same pattern
    91     !
     91
    9292    CHARACTER*(*) :: TARGET
    9393    LOGICAL :: def_val

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/groupe.F

    r5082 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      subroutine groupe(pext,pbaru,pbarv,pbarum,pbarvm,wm)
    55     

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/groupeun.F

    r2597 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE groupeun(jjmax,llmax,q)
    55     

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/guide_mod.F90

    r5093 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44MODULE guide_mod
    55
    … …  
    176176       if (ncidpl==-99) then
    177177          rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
    178           if (rcod/=NF90_NOERR) THEN
     178          if (rcod/=nf90_noerr) THEN
    179179             abort_message=' Nudging error -> no file apbp.nc'
    180180             CALL abort_gcm(modname,abort_message,1)
    … …  
    184184       if (ncidpl==-99) then
    185185          rcod=nf90_open('P.nc',Nf90_NOWRITe,ncidpl)
    186           if (rcod/=NF90_NOERR) THEN
     186          if (rcod/=nf90_noerr) THEN
    187187             abort_message=' Nudging error -> no file P.nc'
    188188             CALL abort_gcm(modname,abort_message,1)
    … …  
    193193           if (ncidpl==-99) then
    194194               rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
    195                if (rcod/=NF90_NOERR) THEN
     195               if (rcod/=nf90_noerr) THEN
    196196                  CALL abort_gcm(modname, &
    197197                       ' Nudging error -> no file u.nc',1)
    … …  
    202202           if (ncidpl==-99) then
    203203               rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
    204                if (rcod/=NF90_NOERR) THEN
     204               if (rcod/=nf90_noerr) THEN
    205205                  CALL abort_gcm(modname, &
    206206                       ' Nudging error -> no file v.nc',1)
    … …  
    210210           if (ncidpl==-99) then
    211211               rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
    212                if (rcod/=NF90_NOERR) THEN
     212               if (rcod/=nf90_noerr) THEN
    213213                  CALL abort_gcm(modname, &
    214214                       ' Nudging error -> no file T.nc',1)
    … …  
    218218           if (ncidpl==-99) then
    219219               rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
    220                if (rcod/=NF90_NOERR) THEN
     220               if (rcod/=nf90_noerr) THEN
    221221                  CALL abort_gcm(modname, &
    222222                       ' Nudging error -> no file hur.nc',1)
    … …  
    227227    endif
    228228    error=nf90_inq_dimid(ncidpl,'LEVEL',rid)
    229     IF (error/=NF90_NOERR) error=nf90_inq_dimid(ncidpl,'PRESSURE',rid)
    230     IF (error/=NF90_NOERR) THEN
     229    IF (error/=nf90_noerr) error=nf90_inq_dimid(ncidpl,'PRESSURE',rid)
     230    IF (error/=nf90_noerr) THEN
    231231        CALL abort_gcm(modname,'Nudging: error reading pressure levels',1)
    232232    ENDIF
    … …  
    11201120             write(*,*) trim(modname)//' Reading nudging on model levels'
    11211121             rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
    1122              IF (rcode/=NF90_NOERR) THEN
     1122             IF (rcode/=nf90_noerr) THEN
    11231123              abort_message='Nudging: error -> no file apbp.nc'
    11241124              CALL abort_gcm(modname,abort_message,1)
    11251125             ENDIF
    11261126             rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
    1127              IF (rcode/=NF90_NOERR) THEN
     1127             IF (rcode/=nf90_noerr) THEN
    11281128              abort_message='Nudging: error -> no AP variable in file apbp.nc'
    11291129              CALL abort_gcm(modname,abort_message,1)
    11301130             ENDIF
    11311131             rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
    1132              IF (rcode/=NF90_NOERR) THEN
     1132             IF (rcode/=nf90_noerr) THEN
    11331133              abort_message='Nudging: error -> no BP variable in file apbp.nc'
    11341134              CALL abort_gcm(modname,abort_message,1)
    … …  
    11401140         if (guide_plevs==2) then
    11411141             rcode = nf90_open('P.nc', nf90_nowrite, ncidp)
    1142              IF (rcode/=NF90_NOERR) THEN
     1142             IF (rcode/=nf90_noerr) THEN
    11431143              abort_message='Nudging: error -> no file P.nc'
    11441144              CALL abort_gcm(modname,abort_message,1)
    11451145             ENDIF
    11461146             rcode = nf90_inq_varid(ncidp, 'PRES', varidp)
    1147              IF (rcode/=NF90_NOERR) THEN
     1147             IF (rcode/=nf90_noerr) THEN
    11481148              abort_message='Nudging: error -> no PRES variable in file P.nc'
    11491149              CALL abort_gcm(modname,abort_message,1)
    … …  
    11561156         if (guide_u) then
    11571157             rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
    1158              IF (rcode/=NF90_NOERR) THEN
     1158             IF (rcode/=nf90_noerr) THEN
    11591159              abort_message='Nudging: error -> no file u.nc'
    11601160              CALL abort_gcm(modname,abort_message,1)
    11611161             ENDIF
    11621162             rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
    1163              IF (rcode/=NF90_NOERR) THEN
     1163             IF (rcode/=nf90_noerr) THEN
    11641164              abort_message='Nudging: error -> no UWND variable in file u.nc'
    11651165              CALL abort_gcm(modname,abort_message,1)
    … …  
    11691169
    11701170             status=nf90_inq_dimid(ncidu, "LONU", dimid)
    1171              status=NF90_INQUIRE_DIMENSION(ncidu,dimid,namedim,lendim)
     1171             status=nf90_inquire_dimension(ncidu,dimid,namedim,lendim)
    11721172             IF (lendim /= iip1) THEN
    11731173                abort_message='dimension LONU different from iip1 in u.nc'
    … …  
    11761176
    11771177             status=nf90_inq_dimid(ncidu, "LATU", dimid)
    1178              status=NF90_INQUIRE_DIMENSION(ncidu,dimid,namedim,lendim)
     1178             status=nf90_inquire_dimension(ncidu,dimid,namedim,lendim)
    11791179             IF (lendim /= jjp1) THEN
    11801180                abort_message='dimension LATU different from jjp1 in u.nc'
    … …  
    11871187         if (guide_v) then
    11881188             rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
    1189              IF (rcode/=NF90_NOERR) THEN
     1189             IF (rcode/=nf90_noerr) THEN
    11901190              abort_message='Nudging: error -> no file v.nc'
    11911191              CALL abort_gcm(modname,abort_message,1)
    11921192             ENDIF
    11931193             rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
    1194              IF (rcode/=NF90_NOERR) THEN
     1194             IF (rcode/=nf90_noerr) THEN
    11951195              abort_message='Nudging: error -> no VWND variable in file v.nc'
    11961196              CALL abort_gcm(modname,abort_message,1)
    … …  
    12001200             
    12011201             status=nf90_inq_dimid(ncidv, "LONV", dimid)
    1202              status=NF90_INQUIRE_DIMENSION(ncidv,dimid,namedim,lendim)
     1202             status=nf90_inquire_dimension(ncidv,dimid,namedim,lendim)
    12031203             
    12041204                IF (lendim /= iip1) THEN
    … …  
    12091209
    12101210             status=nf90_inq_dimid(ncidv, "LATV", dimid)
    1211              status=NF90_INQUIRE_DIMENSION(ncidv,dimid,namedim,lendim)
     1211             status=nf90_inquire_dimension(ncidv,dimid,namedim,lendim)
    12121212             IF (lendim /= jjm) THEN
    12131213                abort_message='dimension LATV different from jjm in v.nc'
    … …  
    12201220         if (guide_T) then
    12211221             rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
    1222              IF (rcode/=NF90_NOERR) THEN
     1222             IF (rcode/=nf90_noerr) THEN
    12231223              abort_message='Nudging: error -> no file T.nc'
    12241224              CALL abort_gcm(modname,abort_message,1)
    12251225             ENDIF
    12261226             rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
    1227              IF (rcode/=NF90_NOERR) THEN
     1227             IF (rcode/=nf90_noerr) THEN
    12281228              abort_message='Nudging: error -> no AIR variable in file T.nc'
    12291229              CALL abort_gcm(modname,abort_message,1)
    … …  
    12331233
    12341234             status=nf90_inq_dimid(ncidt, "LONV", dimid)
    1235              status=NF90_INQUIRE_DIMENSION(ncidt,dimid,namedim,lendim)
     1235             status=nf90_inquire_dimension(ncidt,dimid,namedim,lendim)
    12361236             IF (lendim /= iip1) THEN
    12371237                abort_message='dimension LONV different from iip1 in T.nc'
    … …  
    12401240
    12411241             status=nf90_inq_dimid(ncidt, "LATU", dimid)
    1242              status=NF90_INQUIRE_DIMENSION(ncidt,dimid,namedim,lendim)
     1242             status=nf90_inquire_dimension(ncidt,dimid,namedim,lendim)
    12431243             IF (lendim /= jjp1) THEN
    12441244                abort_message='dimension LATU different from jjp1 in T.nc'
    … …  
    12511251         if (guide_Q) then
    12521252             rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
    1253              IF (rcode/=NF90_NOERR) THEN
     1253             IF (rcode/=nf90_noerr) THEN
    12541254              abort_message='Nudging: error -> no file hur.nc'
    12551255              CALL abort_gcm(modname,abort_message,1)
    12561256             ENDIF
    12571257             rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
    1258              IF (rcode/=NF90_NOERR) THEN
     1258             IF (rcode/=nf90_noerr) THEN
    12591259              abort_message='Nudging: error -> no RH variable in file hur.nc'
    12601260              CALL abort_gcm(modname,abort_message,1)
    … …  
    12641264
    12651265             status=nf90_inq_dimid(ncidQ, "LONV", dimid)
    1266              status=NF90_INQUIRE_DIMENSION(ncidQ,dimid,namedim,lendim)
     1266             status=nf90_inquire_dimension(ncidQ,dimid,namedim,lendim)
    12671267             IF (lendim /= iip1) THEN
    12681268                abort_message='dimension LONV different from iip1 in hur.nc'
    … …  
    12711271
    12721272             status=nf90_inq_dimid(ncidQ, "LATU", dimid)
    1273              status=NF90_INQUIRE_DIMENSION(ncidQ,dimid,namedim,lendim)
     1273             status=nf90_inquire_dimension(ncidQ,dimid,namedim,lendim)
    12741274             IF (lendim /= jjp1) THEN
    12751275                abort_message='dimension LATU different from jjp1 in hur.nc'
    … …  
    12821282         if ((guide_P).OR.(guide_modele)) then
    12831283             rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
    1284              IF (rcode/=NF90_NOERR) THEN
     1284             IF (rcode/=nf90_noerr) THEN
    12851285              abort_message='Nudging: error -> no file ps.nc'
    12861286              CALL abort_gcm(modname,abort_message,1)
    12871287             ENDIF
    12881288             rcode = nf90_inq_varid(ncidps, 'SP', varidps)
    1289              IF (rcode/=NF90_NOERR) THEN
     1289             IF (rcode/=nf90_noerr) THEN
    12901290              abort_message='Nudging: error -> no SP variable in file ps.nc'
    12911291              CALL abort_gcm(modname,abort_message,1)
    … …  
    13011301! Coefs ap, bp pour calcul de la pression aux differents niveaux
    13021302         if (guide_plevs==1) then
    1303              status=NF90_GET_VAR(ncidpl,varidap,apnc,[1],[nlevnc])
    1304              status=NF90_GET_VAR(ncidpl,varidbp,bpnc,[1],[nlevnc])
     1303             status=nf90_get_var(ncidpl,varidap,apnc,[1],[nlevnc])
     1304             status=nf90_get_var(ncidpl,varidbp,bpnc,[1],[nlevnc])
    13051305         ELSEIF (guide_plevs==0) THEN
    1306              status=NF90_GET_VAR(ncidpl,varidpl,apnc,[1],[nlevnc])
     1306             status=nf90_get_var(ncidpl,varidpl,apnc,[1],[nlevnc])
    13071307!FC Pour les corrections la pression est deja en Pascals on commente la ligne ci-dessous
    13081308             IF(convert_Pa) apnc=apnc*100.! conversion en Pascals
    … …  
    13291329! Pression
    13301330     if (guide_plevs==2) then
    1331          status=NF90_GET_VAR(ncidp,varidp,pnat2,start,count)
     1331         status=nf90_get_var(ncidp,varidp,pnat2,start,count)
    13321332         IF (invert_y) THEN
    13331333!           PRINT*,"Invertion impossible actuellement"
    … …  
    13391339!  Vent zonal
    13401340     if (guide_u) then
    1341          status=NF90_GET_VAR(ncidu,varidu,unat2,start,count)
     1341         status=nf90_get_var(ncidu,varidu,unat2,start,count)
    13421342         IF (invert_y) THEN
    13431343           CALL invert_lat(iip1,jjp1,nlevnc,unat2)
    … …  
    13471347!  Temperature
    13481348     if (guide_T) then
    1349          status=NF90_GET_VAR(ncidt,varidt,tnat2,start,count)
     1349         status=nf90_get_var(ncidt,varidt,tnat2,start,count)
    13501350         IF (invert_y) THEN
    13511351           CALL invert_lat(iip1,jjp1,nlevnc,tnat2)
    … …  
    13551355!  Humidite
    13561356     if (guide_Q) then
    1357          status=NF90_GET_VAR(ncidQ,varidQ,qnat2,start,count)
     1357         status=nf90_get_var(ncidQ,varidQ,qnat2,start,count)
    13581358         IF (invert_y) THEN
    13591359           CALL invert_lat(iip1,jjp1,nlevnc,qnat2)
    … …  
    13651365     if (guide_v) then
    13661366         count(2)=jjm
    1367          status=NF90_GET_VAR(ncidv,varidv,vnat2,start,count)
     1367         status=nf90_get_var(ncidv,varidv,vnat2,start,count)
    13681368         IF (invert_y) THEN
    13691369           CALL invert_lat(iip1,jjm,nlevnc,vnat2)
    … …  
    13781378         count(3)=1
    13791379         count(4)=0
    1380          status=NF90_GET_VAR(ncidps,varidps,psnat2,start,count)
     1380         status=nf90_get_var(ncidps,varidps,psnat2,start,count)
    13811381         IF (invert_y) THEN
    13821382           CALL invert_lat(iip1,jjp1,1,psnat2)
    … …  
    14201420           write(*,*)trim(modname)//' Reading nudging on model levels'
    14211421           rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
    1422            IF (rcode/=NF90_NOERR) THEN
     1422           IF (rcode/=nf90_noerr) THEN
    14231423             abort_message='Nudging: error -> no file apbp.nc'
    14241424           CALL abort_gcm(modname,abort_message,1)
    14251425           ENDIF
    14261426           rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
    1427            IF (rcode/=NF90_NOERR) THEN
     1427           IF (rcode/=nf90_noerr) THEN
    14281428             abort_message='Nudging: error -> no AP variable in file apbp.nc'
    14291429           CALL abort_gcm(modname,abort_message,1)
    14301430           ENDIF
    14311431           rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
    1432            IF (rcode/=NF90_NOERR) THEN
     1432           IF (rcode/=nf90_noerr) THEN
    14331433             abort_message='Nudging: error -> no BP variable in file apbp.nc'
    14341434             CALL abort_gcm(modname,abort_message,1)
    … …  
    14391439         if (guide_plevs==2) then
    14401440           rcode = nf90_open('P.nc', nf90_nowrite, ncidp)
    1441            IF (rcode/=NF90_NOERR) THEN
     1441           IF (rcode/=nf90_noerr) THEN
    14421442             abort_message='Nudging: error -> no file P.nc'
    14431443             CALL abort_gcm(modname,abort_message,1)
    14441444           ENDIF
    14451445           rcode = nf90_inq_varid(ncidp, 'PRES', varidp)
    1446            IF (rcode/=NF90_NOERR) THEN
     1446           IF (rcode/=nf90_noerr) THEN
    14471447             abort_message='Nudging: error -> no PRES variable in file P.nc'
    14481448             CALL abort_gcm(modname,abort_message,1)
    … …  
    14541454         if (guide_u) then
    14551455           rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
    1456            IF (rcode/=NF90_NOERR) THEN
     1456           IF (rcode/=nf90_noerr) THEN
    14571457             abort_message='Nudging: error -> no file u.nc'
    14581458             CALL abort_gcm(modname,abort_message,1)
    14591459           ENDIF
    14601460           rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
    1461            IF (rcode/=NF90_NOERR) THEN
     1461           IF (rcode/=nf90_noerr) THEN
    14621462             abort_message='Nudging: error -> no UWND variable in file u.nc'
    14631463             CALL abort_gcm(modname,abort_message,1)
    … …  
    14691469         if (guide_v) then
    14701470           rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
    1471            IF (rcode/=NF90_NOERR) THEN
     1471           IF (rcode/=nf90_noerr) THEN
    14721472             abort_message='Nudging: error -> no file v.nc'
    14731473             CALL abort_gcm(modname,abort_message,1)
    14741474           ENDIF
    14751475           rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
    1476            IF (rcode/=NF90_NOERR) THEN
     1476           IF (rcode/=nf90_noerr) THEN
    14771477             abort_message='Nudging: error -> no VWND variable in file v.nc'
    14781478             CALL abort_gcm(modname,abort_message,1)
    … …  
    14841484         if (guide_T) then
    14851485           rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
    1486            IF (rcode/=NF90_NOERR) THEN
     1486           IF (rcode/=nf90_noerr) THEN
    14871487             abort_message='Nudging: error -> no file T.nc'
    14881488             CALL abort_gcm(modname,abort_message,1)
    14891489           ENDIF
    14901490           rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
    1491            IF (rcode/=NF90_NOERR) THEN
     1491           IF (rcode/=nf90_noerr) THEN
    14921492             abort_message='Nudging: error -> no AIR variable in file T.nc'
    14931493             CALL abort_gcm(modname,abort_message,1)
    … …  
    14991499         if (guide_Q) then
    15001500           rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
    1501            IF (rcode/=NF90_NOERR) THEN
     1501           IF (rcode/=nf90_noerr) THEN
    15021502             abort_message='Nudging: error -> no file hur.nc'
    15031503             CALL abort_gcm(modname,abort_message,1)
    15041504           ENDIF
    15051505           rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
    1506            IF (rcode/=NF90_NOERR) THEN
     1506           IF (rcode/=nf90_noerr) THEN
    15071507             abort_message='Nudging: error -> no RH,variable in file hur.nc'
    15081508             CALL abort_gcm(modname,abort_message,1)
    … …  
    15141514         if ((guide_P).OR.(guide_modele)) then
    15151515           rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
    1516            IF (rcode/=NF90_NOERR) THEN
     1516           IF (rcode/=nf90_noerr) THEN
    15171517             abort_message='Nudging: error -> no file ps.nc'
    15181518             CALL abort_gcm(modname,abort_message,1)
    15191519           ENDIF
    15201520           rcode = nf90_inq_varid(ncidps, 'SP', varidps)
    1521            IF (rcode/=NF90_NOERR) THEN
     1521           IF (rcode/=nf90_noerr) THEN
    15221522             abort_message='Nudging: error -> no SP variable in file ps.nc'
    15231523             CALL abort_gcm(modname,abort_message,1)
    … …  
    15331533! Coefs ap, bp pour calcul de la pression aux differents niveaux
    15341534         if (guide_plevs==1) then
    1535              status=NF90_GET_VAR(ncidpl,varidap,apnc,[1],[nlevnc])
    1536              status=NF90_GET_VAR(ncidpl,varidbp,bpnc,[1],[nlevnc])
     1535             status=nf90_get_var(ncidpl,varidap,apnc,[1],[nlevnc])
     1536             status=nf90_get_var(ncidpl,varidbp,bpnc,[1],[nlevnc])
    15371537         elseif (guide_plevs==0) THEN
    1538              status=NF90_GET_VAR(ncidpl,varidpl,apnc,[1],[nlevnc])
     1538             status=nf90_get_var(ncidpl,varidpl,apnc,[1],[nlevnc])
    15391539             apnc=apnc*100.! conversion en Pascals
    15401540             bpnc(:)=0.
    … …  
    15601560!  Pression
    15611561     if (guide_plevs==2) then
    1562          status=NF90_GET_VAR(ncidp,varidp,zu,start,count)
     1562         status=nf90_get_var(ncidp,varidp,zu,start,count)
    15631563         DO i=1,iip1
    15641564             pnat2(i,:,:)=zu(:,:)
    … …  
    15731573!  Vent zonal
    15741574     if (guide_u) then
    1575          status=NF90_GET_VAR(ncidu,varidu,zu,start,count)
     1575         status=nf90_get_var(ncidu,varidu,zu,start,count)
    15761576         DO i=1,iip1
    15771577             unat2(i,:,:)=zu(:,:)
    … …  
    15861586!  Temperature
    15871587     if (guide_T) then
    1588          status=NF90_GET_VAR(ncidt,varidt,zu,start,count)
     1588         status=nf90_get_var(ncidt,varidt,zu,start,count)
    15891589         DO i=1,iip1
    15901590             tnat2(i,:,:)=zu(:,:)
    … …  
    15991599!  Humidite
    16001600     if (guide_Q) then
    1601          status=NF90_GET_VAR(ncidQ,varidQ,zu,start,count)
     1601         status=nf90_get_var(ncidQ,varidQ,zu,start,count)
    16021602         DO i=1,iip1
    16031603             qnat2(i,:,:)=zu(:,:)
    … …  
    16131613     if (guide_v) then
    16141614         count(2)=jjm
    1615          status=NF90_GET_VAR(ncidv,varidv,zv,start,count)
     1615         status=nf90_get_var(ncidv,varidv,zv,start,count)
    16161616         DO i=1,iip1
    16171617             vnat2(i,:,:)=zv(:,:)
    … …  
    16311631         count(3)=1
    16321632         count(4)=0
    1633          status=NF90_GET_VAR(ncidps,varidps,zu(:,1),start,count)
     1633         status=nf90_get_var(ncidps,varidps,zu(:,1),start,count)
    16341634         DO i=1,iip1
    16351635             psnat2(i,:)=zu(:,1)

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/iniacademic.F90

    r5088 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44SUBROUTINE iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
    55
    … …  
    2222  USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
    2323  USE readTracFiles_mod, ONLY: addPhase
    24   use netcdf, only : NF90_NOWRITE,NF90_OPEN,NF90_NOERR,NF90_INQ_VARID,NF90_CLOSE,NF90_GET_VAR
     24  use netcdf, only : NF90_NOWRITE,NF90_OPEN,nf90_noerr,NF90_INQ_VARID,NF90_CLOSE,nf90_get_var
    2525
    2626  !   Author:    Frederic Hourdin      original: 15/01/93
    … …  
    9595  ! 1. Initializations for Earth-like case
    9696  ! --------------------------------------
    97   !
     97
    9898  ! initialize planet radius, rotation rate,...
    9999  call conf_planete
    … …  
    141141     relief=0.
    142142     ierr = NF90_OPEN ('relief_in.nc', NF90_NOWRITE,nid_relief)
    143      if (ierr==NF90_NOERR) THEN
     143     if (ierr==nf90_noerr) THEN
    144144         ierr=NF90_INQ_VARID(nid_relief,'RELIEF',varid)
    145          if (ierr==NF90_NOERR) THEN
    146               ierr=NF90_GET_VAR(nid_relief,varid,relief(1:iim,1:jjp1))
     145         if (ierr==nf90_noerr) THEN
     146              ierr=nf90_get_var(nid_relief,varid,relief(1:iim,1:jjp1))
    147147              relief(iip1,:)=relief(1,:)
    148148         else

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/integrd.F

    r5091 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE integrd
    55     $  (  nq,vcovm1,ucovm1,tetam1,psm1,massem1,

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/leapfrog.F

    r5091 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44c
    55c

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/logic_mod.F90

    r4996 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44MODULE logic_mod
    55

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/paramet.h

    r1907 r5099  
    1 !
     1
    22! $Header$
    3 !
    4 !
     3
     4
    55!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
    66!                 veillez  n'utiliser que des ! pour les commentaires
    77!                 et  bien positionner les & des lignes de continuation
    88!                 (les placer en colonne 6 et en colonne 73)
    9 !
    10 !
     9
     10
    1111!-----------------------------------------------------------------------
    1212!   INCLUDE 'paramet.h'

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/qminimum.F

    r5082 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      SUBROUTINE qminimum( q,nqtot,deltap )
    55

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/replay3d.F90

    r5090 r5099  
    1 !
     1
    22! $Id: gcm.F90 3579 2019-10-09 13:11:07Z fairhead $
    3 !
    4 !
    5 !
     3
     4
     5
    66PROGRAM replay3d
    77
    … …  
    3030
    3131  !=======================================================================
    32   !
     32
    3333  !   Auteur:  P. Le Van /L. Fairhead/F.Hourdin
    3434  !   -------
    35   !
     35
    3636  !   Objet:
    3737  !   ------
    38   !
     38
    3939  !   GCM LMD nouvelle grille
    40   !
     40
    4141  !=======================================================================
    42   !
     42
    4343  !  ... Dans inigeom , nouveaux calculs pour les elongations  cu , cv
    4444  !      et possibilite d'appeler une fonction f(y)  a derivee tangente
    … …  
    4646  !  ... Possibilite de choisir le schema pour l'advection de
    4747  !        q  , en modifiant iadv dans traceur.def  (MAF,10/02) .
    48   !
     48
    4949  !      Pour Van-Leer + Vapeur d'eau saturee, iadv(1)=4. (F.Codron,10/99)
    5050  !      Pour Van-Leer iadv=10
    51   !
     51
    5252  !-----------------------------------------------------------------------
    5353  !   Declarations:
    … …  
    8585  !  lecture des fichiers gcm.def ou run.def
    8686  !  ---------------------------------------
    87   !
     87
    8888preff=101325.
    8989 pa=50000.
    … …  
    128128  zdtvr    = daysec/REAL(day_step)
    129129
    130   !
    131130  ! on remet le calendrier \`a zero si demande
    132   !
     131
    133132     annee_ref = anneeref
    134133     day_ref = dayref

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/sw_case_williamson91_6.F

    r2600 r5099  
    1 !
     1
    22! $Id $
    3 !
     3
    44      SUBROUTINE sw_case_williamson91_6(vcov,ucov,teta,masse,ps)
    55

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/temps_mod.F90

    r3540 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44MODULE temps_mod
    55

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/top_bound.F

    r5093 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44      SUBROUTINE top_bound(vcov,ucov,teta,masse,dt)
    55     

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/wrgrads.F

    r5082 r5099  
    1 !
     1
    22! $Header$
    3 !
     3
    44      subroutine wrgrads(if,nl,field,name,titlevar)
    55      implicit none

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/write_paramLMDZ_dyn.h

    r3103 r5099  
    136136!     CALL histwrite(nid_ctesGCM, "cycle_diurne", itau_w,
    137137!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    138 !
     138
    139139!     IF(soil_model) THEN
    140140!      zx_tmp_2d(1:iip1,1:jjp1)=1.
    … …  
    144144!     CALL histwrite(nid_ctesGCM, "soil_model", itau_w,
    145145!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    146 !
     146
    147147!     IF(new_oliq) THEN
    148148!      zx_tmp_2d(1:iip1,1:jjp1)=1.
    … …  
    152152!     CALL histwrite(nid_ctesGCM, "new_oliq", itau_w,
    153153!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    154 !
     154
    155155!     IF(ok_orodr) THEN
    156156!      zx_tmp_2d(1:iip1,1:jjp1)=1.
    … …  
    160160!     CALL histwrite(nid_ctesGCM, "ok_orodr", itau_w,
    161161!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    162 !
     162
    163163!     IF(ok_orolf) THEN
    164164!      zx_tmp_2d(1:iip1,1:jjp1)=1.
    … …  
    168168!     CALL histwrite(nid_ctesGCM, "ok_orolf", itau_w,
    169169!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    170 !
     170
    171171!     IF(ok_limitvrai) THEN
    172172!      zx_tmp_2d(1:iip1,1:jjp1)=1.
    … …  
    176176!     CALL histwrite(nid_ctesGCM, "ok_limitvrai", itau_w,
    177177!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    178 !
     178
    179179!     zx_tmp_2d(1:iip1,1:jjp1)=nbapp_rad
    180180!     CALL histwrite(nid_ctesGCM, "nbapp_rad", itau_w,
    181181!    .               zx_tmp_2d,iip1*jjp1,ndex2d)
    182 !
     182
    183183!     zx_tmp_2d(1:iip1,1:jjp1)=iflag_con
    184184!     CALL histwrite(nid_ctesGCM, "iflag_con", itau_w,
Note: See TracChangeset for help on using the changeset viewer.