Changeset 1299 for LMDZ4/branches


Ignore:
Timestamp:
Jan 20, 2010, 3:27:21 PM (15 years ago)
Author:
Laurent Fairhead
Message:

Nettoyage general pour se rapprocher des normes et éviter des erreurs a la
compilation:

  • tous les FLOAT() sont remplacés par des REAL()
  • tous les STOP dans phylmd sont remplacés par des appels à abort_gcm
  • le common control défini dans le fichier control.h est remplacé par le module control_mod pour éviter des messages sur l'alignement des variables dans les déclarations
  • des $Header$ remplacés par des $Id$ pour svn

Quelques remplacements à faire ont pu m'échapper

General cleanup of the code to try and adhere to norms and to prevent some
compilation errors:

  • all FLOAT() instructions have been replaced by REAL() instructions
  • all STOP instructions in phylmd have been replaced by calls to abort_gcm
  • the common block control defined in the control.h file has been replaced by the control_mod to prevent compilation warnings on the alignement of declared variables
  • $Header$ replaced by $Id$ for svn

Some changes which should have been made might have escaped me

Location:
LMDZ4/branches/LMDZ4V5.0-dev/libf
Files:
2 added
147 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/adaptdt.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      subroutine adaptdt(nadv,dtbon,n,pbaru,
    55     c                   masse)
    66
     7      USE control_mod
    78      IMPLICIT NONE
    89
    … …  
    1617#include "logic.h"
    1718#include "temps.h"
    18 #include "control.h"
    1919#include "ener.h"
    2020#include "description.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/advtrac.F

    r1279 r1299  
    1616c
    1717      USE infotrac
     18      USE control_mod
     19 
    1820
    1921      IMPLICIT NONE
    … …  
    2729#include "logic.h"
    2830#include "temps.h"
    29 #include "control.h"
    3031#include "ener.h"
    3132#include "description.h"
    … …  
    121122
    122123      ! ... Flux de masse diaganostiques traceurs
    123       flxw = wg / FLOAT(iapp_tracvl)
     124      flxw = wg / REAL(iapp_tracvl)
    124125
    125126c  test sur l'eventuelle creation de valeurs negatives de la masse

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/bilan_dyn.F

    r1279 r1299  
    423423         Q_cum(:,:,:,iQ)=Q_cum(:,:,:,iQ)/masse_cum(:,:,:)
    424424      enddo
    425       zz=1./float(ncum)
     425      zz=1./REAL(ncum)
    426426      ps_cum=ps_cum*zz
    427427      masse_cum=masse_cum*zz

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/caladvtrac.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    99c
    1010      USE infotrac
     11      USE control_mod
     12 
    1113      IMPLICIT NONE
    1214c
    … …  
    2426#include "paramet.h"
    2527#include "comconst.h"
    26 #include "control.h"
    2728
    2829c   Arguments:

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/calfis.F

    r1279 r1299  
    3131c   .........
    3232      USE infotrac
     33      USE control_mod
     34 
    3335
    3436      IMPLICIT NONE
    … …  
    9698#include "comvert.h"
    9799#include "comgeom2.h"
    98 #include "control.h"
    99100
    100101c    Arguments :

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/conf_gcm.F

    r1279 r1299  
    66      SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
    77c
     8      USE control_mod
    89#ifdef CPP_IOIPSL
    910      use IOIPSL
    … …  
    1112! if not using IOIPSL, we still need to use (a local version of) getin
    1213      use ioipsl_getincom
     14 
    1315#endif
    1416      IMPLICIT NONE
    … …  
    3436#include "dimensions.h"
    3537#include "paramet.h"
    36 #include "control.h"
    3738#include "logic.h"
    3839#include "serre.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/create_etat0_limit.F

    r1279 r1299  
    66! This prog. is designed to work for Earth
    77       USE dimphy
     8       USE control_mod
    89       USE comgeomphy
    910       USE infotrac
    … …  
    3435#include "paramet.h"
    3536#include "indicesol.h"
    36 #include  "control.h"
    3737      REAL :: masque(iip1,jjp1)
    3838!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/defrun.F

    r956 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    66      SUBROUTINE defrun( tapedef, etatinit, clesphy0 )
    77c
     8      USE control_mod
     9 
    810      IMPLICIT NONE
    911c-----------------------------------------------------------------------
    … …  
    2830#include "dimensions.h"
    2931#include "paramet.h"
    30 #include "control.h"
    3132#include "logic.h"
    3233#include "serre.h"
    … …  
    239240       clesphy0(i) = 0.
    240241      ENDDO
    241                           clesphy0(1) = FLOAT( iflag_con )
    242                           clesphy0(2) = FLOAT( nbapp_rad )
     242                          clesphy0(1) = REAL( iflag_con )
     243                          clesphy0(2) = REAL( nbapp_rad )
    243244
    244245       IF( cycle_diurne  ) clesphy0(3) =  1.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/disvert.F

    r1279 r1299  
    111111      snorm  = 0.
    112112      DO l = 1, llm
    113          x = 2.*asin(1.) * (FLOAT(l)-0.5) / float(llm+1)
     113         x = 2.*asin(1.) * (REAL(l)-0.5) / REAL(llm+1)
    114114
    115115         IF (ok_strato) THEN
    … …  
    135135
    136136      DO l=1,llm
    137         nivsigs(l) = FLOAT(l)
     137        nivsigs(l) = REAL(l)
    138138      ENDDO
    139139
    140140      DO l=1,llmp1
    141         nivsig(l)= FLOAT(l)
     141        nivsig(l)= REAL(l)
    142142      ENDDO
    143143

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/dynredem.F

    r1279 r1299  
    88#endif
    99      USE infotrac
     10 
    1011      IMPLICIT NONE
    1112c=======================================================================
    … …  
    7273       tab_cntrl(l) = 0.
    7374      ENDDO
    74        tab_cntrl(1)  = FLOAT(iim)
    75        tab_cntrl(2)  = FLOAT(jjm)
    76        tab_cntrl(3)  = FLOAT(llm)
    77        tab_cntrl(4)  = FLOAT(day_ref)
    78        tab_cntrl(5)  = FLOAT(annee_ref)
     75       tab_cntrl(1)  = REAL(iim)
     76       tab_cntrl(2)  = REAL(jjm)
     77       tab_cntrl(3)  = REAL(llm)
     78       tab_cntrl(4)  = REAL(day_ref)
     79       tab_cntrl(5)  = REAL(annee_ref)
    7980       tab_cntrl(6)  = rad
    8081       tab_cntrl(7)  = omeg
    … …  
    116117      ENDIF
    117118
    118        tab_cntrl(30) = FLOAT(iday_end)
    119        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     119       tab_cntrl(30) = REAL(iday_end)
     120       tab_cntrl(31) = REAL(itau_dyn + itaufin)
    120121c
    121122c    .........................................................
    … …  
    517518     .                     vcov,ucov,teta,q,masse,ps)
    518519      USE infotrac
     520      USE control_mod
     521 
    519522      IMPLICIT NONE
    520523c=================================================================
    … …  
    528531#include "comgeom.h"
    529532#include "temps.h"
    530 #include "control.h"
    531533
    532534      INTEGER l
    … …  
    589591      ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
    590592#endif
    591        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     593       tab_cntrl(31) = REAL(itau_dyn + itaufin)
    592594#ifdef NC_DOUBLE
    593595      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/etat0_netcdf.F

    r1293 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    99      USE ioipsl
    1010      USE dimphy
     11      USE control_mod
    1112      USE infotrac
    1213      USE fonte_neige_mod
    … …  
    8990
    9091#include "comdissnew.h"
    91 #include "control.h"
    9292#include "serre.h"
    9393#include "clesphys.h"
    … …  
    189189      co2_ppm0 = co2_ppm
    190190
    191       dtvr   = daysec/FLOAT(day_step)
     191      dtvr   = daysec/REAL(day_step)
    192192      print*,'dtvr',dtvr
    193193
    … …  
    691691C
    692692      write(*,*)'phystep ',dtvr,iphysiq,nbapp_rad
    693       phystep   = dtvr * FLOAT(iphysiq)
    694       radpas    = NINT (86400./phystep/ FLOAT(nbapp_rad) )
     693      phystep   = dtvr * REAL(iphysiq)
     694      radpas    = NINT (86400./phystep/ REAL(nbapp_rad) )
    695695      write(*,*)'phystep =', phystep, radpas
    696696cIM : lecture de co2_ppm & solaire ds physiq.def

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/extrapol.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44C
    … …  
    158158               jlat = jy(k)
    159159               pwork(i,j) = pwork(i,j)
    160      $                      + pfild(ilon,jlat) * zmask(k)/FLOAT(inbor)
     160     $                      + pfild(ilon,jlat) * zmask(k)/REAL(inbor)
    161161            ENDDO
    162162         ENDIF

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/fluxstokenc.F

    r1279 r1299  
    114114      DO l=1,llm
    115115         DO ij = 1,ip1jmp1
    116             pbaruc(ij,l) = pbaruc(ij,l)/float(istdyn)
    117             tetac(ij,l) = tetac(ij,l)/float(istdyn)
    118             phic(ij,l) = phic(ij,l)/float(istdyn)
     116            pbaruc(ij,l) = pbaruc(ij,l)/REAL(istdyn)
     117            tetac(ij,l) = tetac(ij,l)/REAL(istdyn)
     118            phic(ij,l) = phic(ij,l)/REAL(istdyn)
    119119         ENDDO
    120120         DO ij = 1,ip1jm
    121             pbarvc(ij,l) = pbarvc(ij,l)/float(istdyn)
     121            pbarvc(ij,l) = pbarvc(ij,l)/REAL(istdyn)
    122122         ENDDO
    123123      ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/friction.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c=======================================================================
    55      SUBROUTINE friction(ucov,vcov,pdt)
     6
     7      USE control_mod
     8 
    69      IMPLICIT NONE
    710
    … …  
    2124#include "paramet.h"
    2225#include "comgeom2.h"
    23 #include "control.h"
    2426#include "comconst.h"
    2527

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/fxhyp.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    8989
    9090       DO i = 0, nmax2
    91         xtild(i) = - pi + FLOAT(i) * depi /nmax2
     91        xtild(i) = - pi + REAL(i) * depi /nmax2
    9292       ENDDO
    9393
    … …  
    235235      DO 1500 i = ii1, ii2
    236236
    237       xlon2 = - pi + (FLOAT(i) + xuv - decalx) * depi / FLOAT(iim)
     237      xlon2 = - pi + (REAL(i) + xuv - decalx) * depi / REAL(iim)
    238238
    239239      Xfi    = xlon2
    … …  
    280280550   CONTINUE
    281281
    282        xxprim(i) = depi/ ( FLOAT(iim) * Xprimin )
     282       xxprim(i) = depi/ ( REAL(iim) * Xprimin )
    283283       xvrai(i)  =  xi + xzoom
    284284

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/fxy.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE fxy (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
    … …  
    3232c
    3333       DO j = 1, jjm + 1
    34           rlatu(j) = fy    ( FLOAT( j )        )
    35          yprimu(j) = fyprim( FLOAT( j )        )
     34          rlatu(j) = fy    ( REAL( j )        )
     35         yprimu(j) = fyprim( REAL( j )        )
    3636       ENDDO
    3737
    … …  
    3939       DO j = 1, jjm
    4040
    41          rlatv(j)  = fy    ( FLOAT( j ) + 0.5  )
    42          rlatu1(j) = fy    ( FLOAT( j ) + 0.25 )
    43          rlatu2(j) = fy    ( FLOAT( j ) + 0.75 )
     41         rlatv(j)  = fy    ( REAL( j ) + 0.5  )
     42         rlatu1(j) = fy    ( REAL( j ) + 0.25 )
     43         rlatu2(j) = fy    ( REAL( j ) + 0.75 )
    4444
    45         yprimv(j)  = fyprim( FLOAT( j ) + 0.5  )
    46         yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
    47         yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )
     45        yprimv(j)  = fyprim( REAL( j ) + 0.5  )
     46        yprimu1(j) = fyprim( REAL( j ) + 0.25 )
     47        yprimu2(j) = fyprim( REAL( j ) + 0.75 )
    4848
    4949       ENDDO
    … …  
    5353c
    5454       DO i = 1, iim + 1
    55            rlonv(i)     = fx    (   FLOAT( i )          )
    56            rlonu(i)     = fx    (   FLOAT( i ) + 0.5    )
    57         rlonm025(i)     = fx    (   FLOAT( i ) - 0.25  )
    58         rlonp025(i)     = fx    (   FLOAT( i ) + 0.25  )
     55           rlonv(i)     = fx    (   REAL( i )          )
     56           rlonu(i)     = fx    (   REAL( i ) + 0.5    )
     57        rlonm025(i)     = fx    (   REAL( i ) - 0.25  )
     58        rlonp025(i)     = fx    (   REAL( i ) + 0.25  )
    5959
    60          xprimv  (i)    = fxprim (  FLOAT( i )          )
    61          xprimu  (i)    = fxprim (  FLOAT( i ) + 0.5    )
    62         xprimm025(i)    = fxprim (  FLOAT( i ) - 0.25   )
    63         xprimp025(i)    = fxprim (  FLOAT( i ) + 0.25   )
     60         xprimv  (i)    = fxprim (  REAL( i )          )
     61         xprimu  (i)    = fxprim (  REAL( i ) + 0.5    )
     62        xprimm025(i)    = fxprim (  REAL( i ) - 0.25   )
     63        xprimp025(i)    = fxprim (  REAL( i ) + 0.25   )
    6464       ENDDO
    6565

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/fxysinus.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE fxysinus (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
    … …  
    3232c
    3333       DO j = 1, jjm + 1
    34           rlatu(j) = fy    ( FLOAT( j )        )
    35          yprimu(j) = fyprim( FLOAT( j )        )
     34          rlatu(j) = fy    ( REAL( j )        )
     35         yprimu(j) = fyprim( REAL( j )        )
    3636       ENDDO
    3737
    … …  
    3939       DO j = 1, jjm
    4040
    41          rlatv(j)  = fy    ( FLOAT( j ) + 0.5  )
    42          rlatu1(j) = fy    ( FLOAT( j ) + 0.25 )
    43          rlatu2(j) = fy    ( FLOAT( j ) + 0.75 )
     41         rlatv(j)  = fy    ( REAL( j ) + 0.5  )
     42         rlatu1(j) = fy    ( REAL( j ) + 0.25 )
     43         rlatu2(j) = fy    ( REAL( j ) + 0.75 )
    4444
    45         yprimv(j)  = fyprim( FLOAT( j ) + 0.5  )
    46         yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
    47         yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )
     45        yprimv(j)  = fyprim( REAL( j ) + 0.5  )
     46        yprimu1(j) = fyprim( REAL( j ) + 0.25 )
     47        yprimu2(j) = fyprim( REAL( j ) + 0.75 )
    4848
    4949       ENDDO
    … …  
    5353c
    5454       DO i = 1, iim + 1
    55            rlonv(i)     = fx    (   FLOAT( i )          )
    56            rlonu(i)     = fx    (   FLOAT( i ) + 0.5    )
    57         rlonm025(i)     = fx    (   FLOAT( i ) - 0.25  )
    58         rlonp025(i)     = fx    (   FLOAT( i ) + 0.25  )
     55           rlonv(i)     = fx    (   REAL( i )          )
     56           rlonu(i)     = fx    (   REAL( i ) + 0.5    )
     57        rlonm025(i)     = fx    (   REAL( i ) - 0.25  )
     58        rlonp025(i)     = fx    (   REAL( i ) + 0.25  )
    5959
    60          xprimv  (i)    = fxprim (  FLOAT( i )          )
    61          xprimu  (i)    = fxprim (  FLOAT( i ) + 0.5    )
    62         xprimm025(i)    = fxprim (  FLOAT( i ) - 0.25   )
    63         xprimp025(i)    = fxprim (  FLOAT( i ) + 0.25   )
     60         xprimv  (i)    = fxprim (  REAL( i )          )
     61         xprimu  (i)    = fxprim (  REAL( i ) + 0.5    )
     62        xprimm025(i)    = fxprim (  REAL( i ) - 0.25   )
     63        xprimp025(i)    = fxprim (  REAL( i ) + 0.25   )
    6464       ENDDO
    6565

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/fyhyp.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    7575       depi     = 2. * pi
    7676       pis2     = pi/2.
    77        pisjm    = pi/ FLOAT(jjm)
     77       pisjm    = pi/ REAL(jjm)
    7878       epsilon  = 1.e-3
    7979       y0       =  yzoomdeg * pi/180.
    … …  
    9494
    9595       DO i = 0, nmax2
    96         yt(i) = - pis2  + FLOAT(i)* pi /nmax2
     96        yt(i) = - pis2  + REAL(i)* pi /nmax2
    9797       ENDDO
    9898
    … …  
    210210       DO 1500 j =  1,jlat
    211211        yo1   = 0.
    212         ylon2 =  - pis2 + pisjm * ( FLOAT(j)  + yuv  -1.) 
     212        ylon2 =  - pis2 + pisjm * ( REAL(j)  + yuv  -1.) 
    213213        yfi    = ylon2
    214214c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/gcm.F

    r1279 r1299  
    1515      USE filtreg_mod
    1616      USE infotrac
     17      USE control_mod
    1718
    1819!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    6869#include "logic.h"
    6970#include "temps.h"
    70 #include "control.h"
     71!!!!!!!!!!!#include "control.h"
    7172#include "ener.h"
    7273#include "description.h"
    … …  
    292293      ENDIF
    293294
    294       zdtvr    = daysec/FLOAT(day_step)
     295      zdtvr    = daysec/REAL(day_step)
    295296        IF(dtvr.NE.zdtvr) THEN
    296297         WRITE(lunout,*)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/grid_atob.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE grille_m(imdep, jmdep, xdata, ydata, entree,
    … …  
    717717c Calculs intermediares:
    718718c
    719       xtmp(1) = -180.0 + 360.0/FLOAT(imtmp) / 2.0
     719      xtmp(1) = -180.0 + 360.0/REAL(imtmp) / 2.0
    720720      DO i = 2, imtmp
    721          xtmp(i) = xtmp(i-1) + 360.0/FLOAT(imtmp)
     721         xtmp(i) = xtmp(i-1) + 360.0/REAL(imtmp)
    722722      ENDDO
    723723      DO i = 1, imtmp
    724724         xtmp(i) = xtmp(i) /180.0 * 4.0*ATAN(1.0)
    725725      ENDDO
    726       ytmp(1) = -90.0 + 180.0/FLOAT(jmtmp) / 2.0
     726      ytmp(1) = -90.0 + 180.0/REAL(jmtmp) / 2.0
    727727      DO j = 2, jmtmp
    728          ytmp(j) = ytmp(j-1) + 180.0/FLOAT(jmtmp)
     728         ytmp(j) = ytmp(j-1) + 180.0/REAL(jmtmp)
    729729      ENDDO
    730730      DO j = 1, jmtmp

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/grid_noro.F

    r773 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    9393      xpi=acos(-1.)
    9494      rad    = 6 371 229.
    95       zdeltay=2.*xpi/float(jusn)*rad
     95      zdeltay=2.*xpi/REAL(jusn)*rad
    9696c
    9797c utilise-t'on un masque lu?
    … …  
    215215c  SUMMATION OVER GRIDPOINT AREA
    216216c
    217       zleny=xpi/float(jusn)*rad
    218       xincr=xpi/2./float(jusn)
     217      zleny=xpi/REAL(jusn)*rad
     218      xincr=xpi/2./REAL(jusn)
    219219       DO ii = 1, imar+1
    220220       DO jj = 1, jmar
    … …  
    468468      DO IS=-1,1
    469469        DO JS=-1,1
    470           WEIGHTpb(IS,JS)=1./FLOAT((1+IS**2)*(1+JS**2))
     470          WEIGHTpb(IS,JS)=1./REAL((1+IS**2)*(1+JS**2))
    471471          SUM=SUM+WEIGHTpb(IS,JS)
    472472        ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/grilles_gcm_netcdf.F

    r636 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    218218      open (20,file='grille.dat',form='unformatted',access='direct'
    219219     s      ,recl=4*ip1jmp1)
    220       write(20,rec=1) ((float(mod(i,2)+mod(j,2)),i=1,iip1),j=1,jjp1)
    221       write(20,rec=2) ((float(mod(i,2)*mod(j,2)),i=1,iip1),j=1,jjp1)
     220      write(20,rec=1) ((REAL(mod(i,2)+mod(j,2)),i=1,iip1),j=1,jjp1)
     221      write(20,rec=2) ((REAL(mod(i,2)*mod(j,2)),i=1,iip1),j=1,jjp1)
    222222      do j=2,jjm
    223223         dlat1(j)=180.*(rlatv(j)-rlatv(j-1))/pi
    224 c        dlat2(j)=180.*fyprim(float(j))/pi
     224c        dlat2(j)=180.*fyprim(REAL(j))/pi
    225225      enddo
    226226      do i=2,iip1
    227227         dlon1(i)=180.*(rlonu(i)-rlonu(i-1))/pi
    228 c        dlon2(i)=180.*fxprim(float(i))/pi
     228c        dlon2(i)=180.*fxprim(REAL(i))/pi
    229229      enddo
    230230      do j=2,jjm

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/guide_mod.F90

    r1279 r1299  
    6262  SUBROUTINE guide_init
    6363
     64    USE control_mod
     65
    6466    IMPLICIT NONE
    6567 
    … …  
    6769    INCLUDE "paramet.h"
    6870    INCLUDE "netcdf.inc"
    69     INCLUDE "control.h"
    7071
    7172    INTEGER                :: error,ncidpl,rid,rcod
    … …  
    269270!=======================================================================
    270271  SUBROUTINE guide_main(itau,ucov,vcov,teta,q,masse,ps)
     272
     273    USE control_mod
    271274 
    272275    IMPLICIT NONE
    … …  
    274277    INCLUDE "dimensions.h"
    275278    INCLUDE "paramet.h"
    276     INCLUDE "control.h"
    277279    INCLUDE "comconst.h"
    278280    INCLUDE "comvert.h"
    … …  
    354356      dday_step=real(day_step)
    355357      IF (iguide_read.LT.0) THEN
    356           tau=ditau/dday_step/FLOAT(iguide_read)
     358          tau=ditau/dday_step/REAL(iguide_read)
    357359      ELSE
    358           tau=FLOAT(iguide_read)*ditau/dday_step
     360          tau=REAL(iguide_read)*ditau/dday_step
    359361      ENDIF
    360362      reste=tau-AINT(tau)
    … …  
    541543            ENDDO
    542544        ENDDO
    543         fieldm(:,l)=fieldm(:,l)/FLOAT(imax(typ)-imin(typ)+1)
     545        fieldm(:,l)=fieldm(:,l)/REAL(imax(typ)-imin(typ)+1)
    544546    ! Compute forcing
    545547        DO j=1,hsize

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/infotrac.F90

    r1279 r1299  
    3131
    3232  SUBROUTINE infotrac_init
     33
     34    USE control_mod
     35 
    3336    IMPLICIT NONE
    3437!=======================================================================
    … …  
    4952
    5053    INCLUDE "dimensions.h"
    51     INCLUDE "control.h"
    5254    INCLUDE "iniprint.h"
    5355
    … …  
    217219          new_iq=new_iq+10 ! 9 tracers added
    218220       ELSE
    219           WRITE(lunout,*) 'This choice of advection schema is not available'
     221          WRITE(lunout,*) 'This choice of advection schema is not available',iq,hadv(iq),vadv(iq)
    220222          CALL abort_gcm('infotrac_init','Bad choice of advection schema - 1',1)
    221223       END IF
    … …  
    258260          iadv(new_iq)=11
    259261       ELSE
    260           WRITE(lunout,*)'This choice of advection schema is not available'
     262          WRITE(lunout,*)'This choice of advection schema is not available',iq,hadv(iq),vadv(iq)
    261263          CALL abort_gcm('infotrac_init','Bad choice of advection schema - 2',1)
    262264       END IF

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/iniacademic.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    88      USE filtreg_mod
    99      USE infotrac, ONLY : nqtot
     10      USE control_mod
     11 
    1012
    1113c%W%    %G%
    … …  
    4446#include "ener.h"
    4547#include "temps.h"
    46 #include "control.h"
    4748#include "iniprint.h"
    4849
    … …  
    9394        g      = 9.8
    9495        daysec = 86400.
    95         dtvr    = daysec/FLOAT(day_step)
     96        dtvr    = daysec/REAL(day_step)
    9697        zdtvr=dtvr
    9798        kappa  = 0.2857143

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/iniconst.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE iniconst
     5
     6      USE control_mod
    57
    68      IMPLICIT NONE
    … …  
    1618#include "comconst.h"
    1719#include "temps.h"
    18 #include "control.h"
    1920#include "comvert.h"
    2021

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/inidissip.F

    r1279 r1299  
    1111c   -------------
    1212
     13      USE control_mod
     14
    1315      IMPLICIT NONE
    1416#include "dimensions.h"
    … …  
    1719#include "comconst.h"
    1820#include "comvert.h"
    19 #include "control.h"
    2021#include "logic.h"
    2122
    … …  
    165166
    166167c     IF(.NOT.lstardis) THEN
    167          fact    = rad*24./float(jjm)
     168         fact    = rad*24./REAL(jjm)
    168169         fact    = fact*fact
    169170         PRINT*,'coef u ', fact/cdivu, 1./cdivu

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/inigeom.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    168168c
    169169      IF( nitergdiv.NE.2 ) THEN
    170         gamdi_gdiv = coefdis/ ( float(nitergdiv) -2. )
     170        gamdi_gdiv = coefdis/ ( REAL(nitergdiv) -2. )
    171171      ELSE
    172172        gamdi_gdiv = 0.
    173173      ENDIF
    174174      IF( nitergrot.NE.2 ) THEN
    175         gamdi_grot = coefdis/ ( float(nitergrot) -2. )
     175        gamdi_grot = coefdis/ ( REAL(nitergrot) -2. )
    176176      ELSE
    177177        gamdi_grot = 0.
    178178      ENDIF
    179179      IF( niterh.NE.2 ) THEN
    180         gamdi_h = coefdis/ ( float(niterh) -2. )
     180        gamdi_h = coefdis/ ( REAL(niterh) -2. )
    181181      ELSE
    182182        gamdi_h = 0.
    … …  
    381381       yprp               = yprimu2(j-1)
    382382       rlatp              = rlatu2 (j-1)
    383 ccc       yprp             = fyprim( FLOAT(j) - 0.25 )
    384 ccc       rlatp            = fy    ( FLOAT(j) - 0.25 )
     383ccc       yprp             = fyprim( REAL(j) - 0.25 )
     384ccc       rlatp            = fy    ( REAL(j) - 0.25 )
    385385c
    386386      coslatp             = COS( rlatp )
    … …  
    416416        rlatm    = rlatu1 (  j  )
    417417        yprm     = yprimu1(  j  )
    418 cc         rlatp    = fy    ( FLOAT(j) - 0.25 )
    419 cc         yprp     = fyprim( FLOAT(j) - 0.25 )
    420 cc         rlatm    = fy    ( FLOAT(j) + 0.25 )
    421 cc         yprm     = fyprim( FLOAT(j) + 0.25 )
     418cc         rlatp    = fy    ( REAL(j) - 0.25 )
     419cc         yprp     = fyprim( REAL(j) - 0.25 )
     420cc         rlatm    = fy    ( REAL(j) + 0.25 )
     421cc         yprm     = fyprim( REAL(j) + 0.25 )
    422422
    423423         coslatm  = COS( rlatm )

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/integrd.F

    r1279 r1299  
    55     $  (  nq,vcovm1,ucovm1,tetam1,psm1,massem1,
    66     $     dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps,masse,phis,finvmaold )
     7
     8      USE control_mod
    79
    810      IMPLICIT NONE
    … …  
    3234#include "temps.h"
    3335#include "serre.h"
    34 #include "control.h"
    3536
    3637c   Arguments:

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/interp_horiz.F

    r616 r1299  
    11c
    2 c $Header$
     2c $Id$
    33c
    44      subroutine interp_horiz (varo,varn,imo,jmo,imn,jmn,lm,
    … …  
    101101           end do
    102102           do ii =1, imn+1
    103              varn(ii,1,l) = totn/float(imn+1)
    104              varn(ii,jmn+1,l) = tots/float(imn+1)
     103             varn(ii,1,l) = totn/REAL(imn+1)
     104             varn(ii,jmn+1,l) = tots/REAL(imn+1)
    105105           end do
    106106       end do

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/interpre.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44       subroutine interpre(q,qppm,w,fluxwppm,masse,
    55     s            apppm,bpppm,massebx,masseby,pbaru,pbarv,
    66     s            unatppm,vnatppm,psppm)
     7
     8      USE control_mod
    79
    810       implicit none
    … …  
    1719#include "logic.h"
    1820#include "temps.h"
    19 #include "control.h"
    2021#include "ener.h"
    2122#include "description.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/juldate.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44        subroutine juldate(ian,imoi,ijou,oh,om,os,tjd,tjdsec)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/leapfrog.F

    r1286 r1299  
    1515      USE guide_mod, ONLY : guide_main
    1616      USE write_field
     17      USE control_mod
    1718      IMPLICIT NONE
    1819
    … …  
    5657#include "logic.h"
    5758#include "temps.h"
    58 #include "control.h"
    5959#include "ener.h"
    6060#include "description.h"
    … …  
    198198      itau = 0
    199199c$$$      iday = day_ini+itau/day_step
    200 c$$$      time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
     200c$$$      time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
    201201c$$$         IF(time.GT.1.) THEN
    202202c$$$          time = time-1.
    … …  
    523523              itau= itau + 1
    524524c$$$              iday= day_ini+itau/day_step
    525 c$$$              time= FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
     525c$$$              time= REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
    526526c$$$                IF(time.GT.1.) THEN
    527527c$$$                  time = time-1.
    … …  
    646646             itau =  itau + 1
    647647c$$$             iday = day_ini+itau/day_step
    648 c$$$             time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
     648c$$$             time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
    649649c$$$
    650650c$$$                  IF(time.GT.1.) THEN

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/limit_netcdf.F

    r1293 r1299  
    1010      use phys_state_var_mod , ONLY : pctsrf
    1111      use inter_barxy_m, only: inter_barxy
     12      USE control_mod
    1213
    1314      IMPLICIT none
    … …  
    2930#include "dimensions.h"
    3031#include "paramet.h"
    31 #include "control.h"
    3232#include "logic.h"
    3333#include "netcdf.inc"
    … …  
    148148      kappa  = 0.2857143
    149149      cpp    = 1004.70885
    150       dtvr    = daysec/FLOAT(day_step)
     150      dtvr    = daysec/REAL(day_step)
    151151      CALL inigeom
    152152c
    … …  
    418418          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    419419          DO k = 1, 360
    420             time = FLOAT(k-1)
     420            time = REAL(k-1)
    421421            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    422422            champan(i,j,k) = by
    … …  
    632632          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    633633          DO k = 1, 360
    634             time = FLOAT(k-1)
     634            time = REAL(k-1)
    635635            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    636636            champan(i,j,k) = by
    … …  
    951951          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    952952          DO k = 1, 360
    953             time = FLOAT(k-1)
     953            time = REAL(k-1)
    954954            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    955955            champan(i,j,k) = by
    … …  
    11461146          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    11471147          DO k = 1, 360
    1148             time = FLOAT(k-1)
     1148            time = REAL(k-1)
    11491149            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    11501150            champan(i,j,k) = by
    … …  
    13001300      ierr = NF_PUT_VARA_DOUBLE (nid,id_RUG,debut,epais,phy_rug(1,k))
    13011301#else
    1302       ierr = NF_PUT_VAR1_REAL (nid,id_tim,k,FLOAT(k))
     1302      ierr = NF_PUT_VAR1_REAL (nid,id_tim,k,REAL(k))
    13031303      IF (newlmt ) THEN
    13041304          ierr = NF_PUT_VARA_REAL (nid,id_FOCE,debut,epais

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/ppm3d.F

    r695 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44
    … …  
    345345C
    346346      PI = 4. * ATAN(1.)
    347       DL = 2.*PI / float(IMR)
    348       DP =    PI / float(JMR)
     347      DL = 2.*PI / REAL(IMR)
     348      DP =    PI / REAL(JMR)
    349349C
    350350      if(IGD.eq.0) then
    … …  
    388388      ZTC  = acos(CR1) * (180./PI)
    389389C
    390       JS0 = float(JMR)*(90.-ZTC)/180. + 2
     390      JS0 = REAL(JMR)*(90.-ZTC)/180. + 2
    391391      JS0 = max(JS0, J1+1)
    392392      IML = min(6*JS0/(J1-1)+2, 4*IMR/5)
    … …  
    628628C Contribution from the N-S advection
    629629      do i=1,imr*(j2-j1+1)
    630       JT = float(J1) - VA(i,j1)
     630      JT = REAL(J1) - VA(i,j1)
    631631      wk1(i,j1,2) = VA(i,j1) * (q(i,jt,k,IC) - q(i,jt+1,k,IC))
    632632      enddo
    … …  
    949949      IF(IORD.eq.1 .or. j.eq.j1. or. j.eq.j2) THEN
    950950      DO 1406 i=1,IMR
    951       iu = float(i) - uc(i,j)
     951      iu = REAL(i) - uc(i,j)
    9529521406  fx1(i) = qtmp(iu)
    953953      ELSE
    … …  
    957957      if(IORD.eq.2 .or. j.le.j1vl .or. j.ge.j2vl) then
    958958      DO 1408 i=1,IMR
    959       iu = float(i) - uc(i,j)
     959      iu = REAL(i) - uc(i,j)
    9609601408  fx1(i) = qtmp(iu) + DC(iu)*(sign(1.,uc(i,j))-uc(i,j))
    961961      else
    … …  
    11111111      if(JORD.eq.1) then
    11121112      DO 1000 i=1,len
    1113       JT = float(J1) - VC(i,J1)
     1113      JT = REAL(J1) - VC(i,J1)
    111411141000  fx(i,j1) = p(i,JT)
    11151115      else
    … …  
    11231123      else
    11241124      DO 1200 i=1,len
    1125       JT = float(J1) - VC(i,J1)
     1125      JT = REAL(J1) - VC(i,J1)
    112611261200  fx(i,j1) = p(i,JT) + (sign(1.,VC(i,j1))-VC(i,j1))*DC2(i,JT)
    11271127      endif
    … …  
    13581358        do j=j1-1,j2+1
    13591359      do i=1,imr
    1360       JP = float(j)-VA(i,j)
     1360      JP = REAL(j)-VA(i,j)
    13611361      ady(i,j) = VA(i,j)*(wk(i,jp)-wk(i,jp+1))
    13621362      enddo
    … …  
    15821582      JMR = JNP-1
    15831583      do 55 j=2,JNP
    1584         ph5  =  -0.5*PI + (FLOAT(J-1)-0.5)*DP
     1584        ph5  =  -0.5*PI + (REAL(J-1)-0.5)*DP
    1585158555      cose(j) = cos(ph5)
    15861586C
    … …  
    18341834C
    18351835c      if(first) then
    1836       DP = 4.*ATAN(1.)/float(JNP-1)
     1836      DP = 4.*ATAN(1.)/REAL(JNP-1)
    18371837      CAP1 = IMR*(1.-COS((j1-1.5)*DP))/DP
    18381838c      first = .false.
    … …  
    18891889C Check Poles.
    18901890      if(q(1,1).lt.0.) then
    1891       dq = q(1,1)*cap1/float(IMR)*acosp(j1)
     1891      dq = q(1,1)*cap1/REAL(IMR)*acosp(j1)
    18921892      do i=1,imr
    18931893      q(i,1) = 0.
    … …  
    18981898C
    18991899      if(q(1,JNP).lt.0.) then
    1900       dq = q(1,JNP)*cap1/float(IMR)*acosp(j2)
     1900      dq = q(1,JNP)*cap1/REAL(IMR)*acosp(j2)
    19011901      do i=1,imr
    19021902      q(i,JNP) = 0.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/ran1.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      FUNCTION RAN1(IDUM)
    … …  
    2020          IX1=MOD(IA1*IX1+IC1,M1)
    2121          IX2=MOD(IA2*IX2+IC2,M2)
    22           R(J)=(FLOAT(IX1)+FLOAT(IX2)*RM2)*RM1
     22          R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
    232311      CONTINUE
    2424        IDUM=1
    … …  
    3030      IF(J.GT.97.OR.J.LT.1)PAUSE
    3131      RAN1=R(J)
    32       R(J)=(FLOAT(IX1)+FLOAT(IX2)*RM2)*RM1
     32      R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
    3333      RETURN
    3434      END

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/sortvarc.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE sortvarc
    … …  
    5959
    6060       dtvrs1j   = dtvr/daysec
    61        rjour     = FLOAT( INT( itau * dtvrs1j ))
     61       rjour     = REAL( INT( itau * dtvrs1j ))
    6262       heure     = ( itau*dtvrs1j-rjour ) * 24.
    6363       imjmp1    = iim * jjp1
    … …  
    129129      ang   = SSUM(     llm,  angl, 1 )
    130130
    131 c      rday = FLOAT(INT ( day_ini + time ))
     131c      rday = REAL(INT ( day_ini + time ))
    132132c
    133        rday = FLOAT(INT(time-jD_ref-jH_ref))
     133       rday = REAL(INT(time-jD_ref-jH_ref))
    134134      IF(ptot0.eq.0.)  THEN
    135135         PRINT 3500, itau, rday, heure,time

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/sortvarc0.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE sortvarc0
    … …  
    6060
    6161       dtvrs1j   = dtvr/daysec
    62        rjour     = FLOAT( INT( itau * dtvrs1j ))
     62       rjour     = REAL( INT( itau * dtvrs1j ))
    6363       heure     = ( itau*dtvrs1j-rjour ) * 24.
    6464       imjmp1    = iim * jjp1
    … …  
    130130      ang0   = SSUM(     llm,  angl, 1 )
    131131
    132       rday = FLOAT(INT (time ))
     132      rday = REAL(INT (time ))
    133133c
    134134      PRINT 3500, itau, rday, heure, time

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/startvar.F

    r1293 r1299  
    352352      phis = phis * 9.81
    353353    !
    354 !      masque(:,:) = FLOAT(tmp_int(:,:))
     354!      masque(:,:) = REAL(tmp_int(:,:))
    355355    !
    356356    !  Compute surface roughness

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/tourabs.F

    r644 r1299  
    5757        ELSE
    5858         rot( ij,l ) = (vcov(ij+1,l)/cv(ij+1)-vcov(ij,l)/cv(ij))/
    59      $                 (2.*pi*RAD*cos(rlatv(j)))*float(iim)
     59     $                 (2.*pi*RAD*cos(rlatv(j)))*REAL(iim)
    6060     $                +(ucov(ij+iip1,l)/cu(ij+iip1)-ucov(ij,l)/cu(ij))/
    61      $                 (pi*RAD)*(float(jjm)-1.)
     61     $                 (pi*RAD)*(REAL(jjm)-1.)
    6262c
    6363        ENDIF

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/traceurpole.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44          subroutine traceurpole(q,masse)
     5
     6      USE control_mod
    57
    68          implicit none
    … …  
    1517#include "logic.h"
    1618#include "temps.h"
    17 #include "control.h"
    1819#include "ener.h"
    1920#include "description.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/ugeostr.F

    r1279 r1299  
    4040            DO i=1,iim
    4141               u(i,j,l)=fact*(phi(i,j+1,l)-phi(i,j,l))
    42                um(j,l)=um(j,l)+u(i,j,l)/float(iim)
     42               um(j,l)=um(j,l)+u(i,j,l)/REAL(iim)
    4343            ENDDO
    4444         ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/write_paramLMDZ_dyn.h

    r1279 r1299  
    77      itau_w=itau_dyn+itau
    88c
    9       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(prt_level)
     9      zx_tmp_2d(1:iip1,1:jjp1)=REAL(prt_level)
    1010      CALL histwrite(nid_ctesGCM, "prt_level", itau_w,
    1111     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    1212c
    13       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(dayref)
     13      zx_tmp_2d(1:iip1,1:jjp1)=REAL(dayref)
    1414      CALL histwrite(nid_ctesGCM, "dayref", itau_w,
    1515     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    1616c
    17       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(anneeref)
     17      zx_tmp_2d(1:iip1,1:jjp1)=REAL(anneeref)
    1818      CALL histwrite(nid_ctesGCM, "anneeref", itau_w,
    1919     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    2020c
    21       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(raz_date)
     21      zx_tmp_2d(1:iip1,1:jjp1)=REAL(raz_date)
    2222      CALL histwrite(nid_ctesGCM, "raz_date", itau_w,
    2323     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    2424c
    25       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nday)
     25      zx_tmp_2d(1:iip1,1:jjp1)=REAL(nday)
    2626      CALL histwrite(nid_ctesGCM, "nday", itau_w,
    2727     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    2828c
    29       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(day_step)
     29      zx_tmp_2d(1:iip1,1:jjp1)=REAL(day_step)
    3030      CALL histwrite(nid_ctesGCM, "day_step", itau_w,
    3131     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    3232c
    33       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iperiod)
     33      zx_tmp_2d(1:iip1,1:jjp1)=REAL(iperiod)
    3434      CALL histwrite(nid_ctesGCM, "iperiod", itau_w,
    3535     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    3636c
    37       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iapp_tracvl)
     37      zx_tmp_2d(1:iip1,1:jjp1)=REAL(iapp_tracvl)
    3838      CALL histwrite(nid_ctesGCM, "iapp_tracvl", itau_w,
    3939     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    4040c
    41       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iconser)
     41      zx_tmp_2d(1:iip1,1:jjp1)=REAL(iconser)
    4242      CALL histwrite(nid_ctesGCM, "iconser", itau_w,
    4343     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    4444c
    45       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iecri)
     45      zx_tmp_2d(1:iip1,1:jjp1)=REAL(iecri)
    4646      CALL histwrite(nid_ctesGCM, "iecri", itau_w,
    4747     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    … …  
    5151     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    5252c
    53       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(idissip)
     53      zx_tmp_2d(1:iip1,1:jjp1)=REAL(idissip)
    5454      CALL histwrite(nid_ctesGCM, "idissip", itau_w,
    5555     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    … …  
    6363     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    6464c
    65       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergdiv)
     65      zx_tmp_2d(1:iip1,1:jjp1)=REAL(nitergdiv)
    6666      CALL histwrite(nid_ctesGCM, "nitergdiv", itau_w,
    6767     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    6868c
    69       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(nitergrot)
     69      zx_tmp_2d(1:iip1,1:jjp1)=REAL(nitergrot)
    7070      CALL histwrite(nid_ctesGCM, "nitergrot", itau_w,
    7171     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    7272c
    73       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(niterh)
     73      zx_tmp_2d(1:iip1,1:jjp1)=REAL(niterh)
    7474      CALL histwrite(nid_ctesGCM, "niterh", itau_w,
    7575     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    … …  
    118118     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    119119c
    120       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iflag_phys)
     120      zx_tmp_2d(1:iip1,1:jjp1)=REAL(iflag_phys)
    121121      CALL histwrite(nid_ctesGCM, "iflag_phys", itau_w,
    122122     .               zx_tmp_2d,iip1*jjp1,ndex2d)
    123123c
    124       zx_tmp_2d(1:iip1,1:jjp1)=FLOAT(iphysiq)
     124      zx_tmp_2d(1:iip1,1:jjp1)=REAL(iphysiq)
    125125      CALL histwrite(nid_ctesGCM, "iphysiq", itau_w,
    126126     .               zx_tmp_2d,iip1*jjp1,ndex2d)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/adaptdt.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      subroutine adaptdt(nadv,dtbon,n,pbaru,
    55     c                   masse)
     6
     7      USE control_mod
    68
    79      IMPLICIT NONE
    … …  
    1618#include "logic.h"
    1719#include "temps.h"
    18 #include "control.h"
    1920#include "ener.h"
    2021#include "description.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/advtrac_p.F

    r1146 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    2323      USE times
    2424      USE infotrac
     25      USE control_mod
    2526      IMPLICIT NONE
    2627c
    … …  
    3334#include "logic.h"
    3435#include "temps.h"
    35 #include "control.h"
    3636#include "ener.h"
    3737#include "description.h"
    … …  
    215215         ijb=ij_begin
    216216         ije=ij_end
    217          flxw(ijb:ije,1:llm)=wg(ijb:ije,1:llm)/FLOAT(iapp_tracvl)
     217         flxw(ijb:ije,1:llm)=wg(ijb:ije,1:llm)/REAL(iapp_tracvl)
    218218
    219219c  test sur l'eventuelle creation de valeurs negatives de la masse

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/bilan_dyn_p.F

    r1279 r1299  
    511511     .                        /masse_cum(:,jjb:jje,:)
    512512      enddo
    513       zz=1./float(ncum)
     513      zz=1./REAL(ncum)
    514514
    515515      jjb=jj_begin

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/caladvtrac_p.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    99      USE parallel
    1010      USE infotrac
     11      USE control_mod
    1112c
    1213      IMPLICIT NONE
    … …  
    2526#include "paramet.h"
    2627#include "comconst.h"
    27 #include "control.h"
    2828
    2929c   Arguments:

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/calfis_p.F

    r1279 r1299  
    4141      USE IOPHY
    4242      USE infotrac
     43      USE control_mod
    4344
    4445      IMPLICIT NONE
    … …  
    107108#include "comvert.h"
    108109#include "comgeom2.h"
    109 #include "control.h"
    110110#ifdef CPP_MPI
    111111      include 'mpif.h'

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/conf_gcm.F

    r1279 r1299  
    1616      use mod_hallo, ONLY : use_mpi_alloc
    1717      use parallel, ONLY : omp_chunk
     18      USE control_mod
    1819      IMPLICIT NONE
    1920c-----------------------------------------------------------------------
    … …  
    3839#include "dimensions.h"
    3940#include "paramet.h"
    40 #include "control.h"
    4141#include "logic.h"
    4242#include "serre.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/create_etat0_limit.F

    r1279 r1299  
    1010       USE mod_const_mpi
    1111       USE infotrac
     12      USE control_mod
    1213#ifdef CPP_IOIPSL
    1314       use ioipsl, only: ioconf_calendar
    … …  
    3637#include "paramet.h"
    3738#include "indicesol.h"
    38 #include  "control.h"
    3939      REAL :: masque(iip1,jjp1)
    4040!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/defrun.F

    r985 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    66      SUBROUTINE defrun( tapedef, etatinit, clesphy0 )
    77c
     8      USE control_mod
    89      IMPLICIT NONE
    910c-----------------------------------------------------------------------
    … …  
    2829#include "dimensions.h"
    2930#include "paramet.h"
    30 #include "control.h"
    3131#include "logic.h"
    3232#include "serre.h"
    … …  
    241241       clesphy0(i) = 0.
    242242      ENDDO
    243                           clesphy0(1) = FLOAT( iflag_con )
    244                           clesphy0(2) = FLOAT( nbapp_rad )
     243                          clesphy0(1) = REAL( iflag_con )
     244                          clesphy0(2) = REAL( nbapp_rad )
    245245
    246246       IF( cycle_diurne  ) clesphy0(3) =  1.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/disvert.F

    r1279 r1299  
    111111      snorm  = 0.
    112112      DO l = 1, llm
    113          x = 2.*asin(1.) * (FLOAT(l)-0.5) / float(llm+1)
     113         x = 2.*asin(1.) * (REAL(l)-0.5) / REAL(llm+1)
    114114
    115115         IF (ok_strato) THEN
    … …  
    135135
    136136      DO l=1,llm
    137         nivsigs(l) = FLOAT(l)
     137        nivsigs(l) = REAL(l)
    138138      ENDDO
    139139
    140140      DO l=1,llmp1
    141         nivsig(l)= FLOAT(l)
     141        nivsig(l)= REAL(l)
    142142      ENDDO
    143143

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/dynredem.F

    r1279 r1299  
    7272       tab_cntrl(l) = 0.
    7373      ENDDO
    74        tab_cntrl(1)  = FLOAT(iim)
    75        tab_cntrl(2)  = FLOAT(jjm)
    76        tab_cntrl(3)  = FLOAT(llm)
    77        tab_cntrl(4)  = FLOAT(day_ref)
    78        tab_cntrl(5)  = FLOAT(annee_ref)
     74       tab_cntrl(1)  =  REAL(iim)
     75       tab_cntrl(2)  =  REAL(jjm)
     76       tab_cntrl(3)  =  REAL(llm)
     77       tab_cntrl(4)  =  REAL(day_ref)
     78       tab_cntrl(5)  =  REAL(annee_ref)
    7979       tab_cntrl(6)  = rad
    8080       tab_cntrl(7)  = omeg
    … …  
    116116      ENDIF
    117117
    118        tab_cntrl(30) = FLOAT(iday_end)
    119        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     118       tab_cntrl(30) =  REAL(iday_end)
     119       tab_cntrl(31) =  REAL(itau_dyn + itaufin)
    120120c
    121121c    .........................................................
    … …  
    517517     .                     vcov,ucov,teta,q,masse,ps)
    518518      USE infotrac
     519      USE control_mod
    519520      IMPLICIT NONE
    520521c=================================================================
    … …  
    528529#include "comgeom.h"
    529530#include "temps.h"
    530 #include "control.h"
    531531
    532532      INTEGER l
    … …  
    589589      ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
    590590#endif
    591        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     591       tab_cntrl(31) =  REAL(itau_dyn + itaufin)
    592592#ifdef NC_DOUBLE
    593593      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/dynredem_p.F

    r1279 r1299  
    7474       tab_cntrl(l) = 0.
    7575      ENDDO
    76        tab_cntrl(1)  = FLOAT(iim)
    77        tab_cntrl(2)  = FLOAT(jjm)
    78        tab_cntrl(3)  = FLOAT(llm)
    79        tab_cntrl(4)  = FLOAT(day_ref)
    80        tab_cntrl(5)  = FLOAT(annee_ref)
     76       tab_cntrl(1)  =  REAL(iim)
     77       tab_cntrl(2)  =  REAL(jjm)
     78       tab_cntrl(3)  =  REAL(llm)
     79       tab_cntrl(4)  =  REAL(day_ref)
     80       tab_cntrl(5)  =  REAL(annee_ref)
    8181       tab_cntrl(6)  = rad
    8282       tab_cntrl(7)  = omeg
    … …  
    118118      ENDIF
    119119
    120        tab_cntrl(30) = FLOAT(iday_end)
    121        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     120       tab_cntrl(30) =  REAL(iday_end)
     121       tab_cntrl(31) =  REAL(itau_dyn + itaufin)
    122122c
    123123c    .........................................................
    … …  
    521521      USE parallel
    522522      USE infotrac
     523      USE control_mod
    523524      IMPLICIT NONE
    524525c=================================================================
    … …  
    532533#include "comgeom.h"
    533534#include "temps.h"
    534 #include "control.h"
    535535
    536536      INTEGER l
    … …  
    608608      ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
    609609#endif
    610        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
     610       tab_cntrl(31) =  REAL(itau_dyn + itaufin)
    611611#ifdef NC_DOUBLE
    612612      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tab_cntrl)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/etat0_netcdf.F

    r1293 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    1414      USE phys_state_var_mod
    1515      USE filtreg_mod
     16      USE control_mod
    1617      use regr_lat_time_climoz_m, only: regr_lat_time_climoz
    1718      use conf_phys_m, only: conf_phys
    … …  
    8990
    9091#include "comdissnew.h"
    91 #include "control.h"
    9292#include "serre.h"
    9393#include "clesphys.h"
    … …  
    189189      co2_ppm0 = co2_ppm
    190190
    191       dtvr   = daysec/FLOAT(day_step)
     191      dtvr   = daysec/REAL(day_step)
    192192      print*,'dtvr',dtvr
    193193
    … …  
    691691C
    692692      write(*,*)'phystep ',dtvr,iphysiq,nbapp_rad
    693       phystep   = dtvr * FLOAT(iphysiq)
    694       radpas    = NINT (86400./phystep/ FLOAT(nbapp_rad) )
     693      phystep   = dtvr * REAL(iphysiq)
     694      radpas    = NINT (86400./phystep/ REAL(nbapp_rad) )
    695695      write(*,*)'phystep =', phystep, radpas
    696696cIM : lecture de co2_ppm & solaire ds physiq.def

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/extrapol.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44C
    … …  
    158158               jlat = jy(k)
    159159               pwork(i,j) = pwork(i,j)
    160      $                      + pfild(ilon,jlat) * zmask(k)/FLOAT(inbor)
     160     $                      + pfild(ilon,jlat) * zmask(k)/ REAL(inbor)
    161161            ENDDO
    162162         ENDIF

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/fluxstokenc_p.F

    r1279 r1299  
    153153      DO l=1,llm
    154154         DO ij = ijb,ije
    155             pbaruc(ij,l) = pbaruc(ij,l)/float(istdyn)
    156             tetac(ij,l) = tetac(ij,l)/float(istdyn)
    157             phic(ij,l) = phic(ij,l)/float(istdyn)
     155            pbaruc(ij,l) = pbaruc(ij,l)/REAL(istdyn)
     156            tetac(ij,l) = tetac(ij,l)/REAL(istdyn)
     157            phic(ij,l) = phic(ij,l)/REAL(istdyn)
    158158         ENDDO
    159159      ENDDO
    … …  
    165165      DO l=1,llm
    166166          DO ij = ijb,ije
    167             pbarvc(ij,l) = pbarvc(ij,l)/float(istdyn)
     167            pbarvc(ij,l) = pbarvc(ij,l)/REAL(istdyn)
    168168         ENDDO
    169169      ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/friction_p.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c=======================================================================
    55      SUBROUTINE friction_p(ucov,vcov,pdt)
    66      USE parallel
     7      USE control_mod
    78      IMPLICIT NONE
    89
    … …  
    2223#include "paramet.h"
    2324#include "comgeom2.h"
    24 #include "control.h"
    2525#include "comconst.h"
    2626

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/fxhyp.F

    r764 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    8989
    9090       DO i = 0, nmax2
    91         xtild(i) = - pi + FLOAT(i) * depi /nmax2
     91        xtild(i) = - pi + REAL(i) * depi /nmax2
    9292       ENDDO
    9393
    … …  
    235235      DO 1500 i = ii1, ii2
    236236
    237       xlon2 = - pi + (FLOAT(i) + xuv - decalx) * depi / FLOAT(iim)
     237      xlon2 = - pi + (REAL(i) + xuv - decalx) * depi / REAL(iim)
    238238
    239239      Xfi    = xlon2
    … …  
    280280550   CONTINUE
    281281
    282        xxprim(i) = depi/ ( FLOAT(iim) * Xprimin )
     282       xxprim(i) = depi/ ( REAL(iim) * Xprimin )
    283283       xvrai(i)  =  xi + xzoom
    284284

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/fxy.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE fxy (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
    … …  
    3232c
    3333       DO j = 1, jjm + 1
    34           rlatu(j) = fy    ( FLOAT( j )        )
    35          yprimu(j) = fyprim( FLOAT( j )        )
     34          rlatu(j) = fy    ( REAL( j )        )
     35         yprimu(j) = fyprim( REAL( j )        )
    3636       ENDDO
    3737
    … …  
    3939       DO j = 1, jjm
    4040
    41          rlatv(j)  = fy    ( FLOAT( j ) + 0.5  )
    42          rlatu1(j) = fy    ( FLOAT( j ) + 0.25 )
    43          rlatu2(j) = fy    ( FLOAT( j ) + 0.75 )
     41         rlatv(j)  = fy    ( REAL( j ) + 0.5  )
     42         rlatu1(j) = fy    ( REAL( j ) + 0.25 )
     43         rlatu2(j) = fy    ( REAL( j ) + 0.75 )
    4444
    45         yprimv(j)  = fyprim( FLOAT( j ) + 0.5  )
    46         yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
    47         yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )
     45        yprimv(j)  = fyprim( REAL( j ) + 0.5  )
     46        yprimu1(j) = fyprim( REAL( j ) + 0.25 )
     47        yprimu2(j) = fyprim( REAL( j ) + 0.75 )
    4848
    4949       ENDDO
    … …  
    5353c
    5454       DO i = 1, iim + 1
    55            rlonv(i)     = fx    (   FLOAT( i )          )
    56            rlonu(i)     = fx    (   FLOAT( i ) + 0.5    )
    57         rlonm025(i)     = fx    (   FLOAT( i ) - 0.25  )
    58         rlonp025(i)     = fx    (   FLOAT( i ) + 0.25  )
     55           rlonv(i)     = fx    (   REAL( i )          )
     56           rlonu(i)     = fx    (   REAL( i ) + 0.5    )
     57        rlonm025(i)     = fx    (   REAL( i ) - 0.25  )
     58        rlonp025(i)     = fx    (   REAL( i ) + 0.25  )
    5959
    60          xprimv  (i)    = fxprim (  FLOAT( i )          )
    61          xprimu  (i)    = fxprim (  FLOAT( i ) + 0.5    )
    62         xprimm025(i)    = fxprim (  FLOAT( i ) - 0.25   )
    63         xprimp025(i)    = fxprim (  FLOAT( i ) + 0.25   )
     60         xprimv  (i)    = fxprim (  REAL( i )          )
     61         xprimu  (i)    = fxprim (  REAL( i ) + 0.5    )
     62        xprimm025(i)    = fxprim (  REAL( i ) - 0.25   )
     63        xprimp025(i)    = fxprim (  REAL( i ) + 0.25   )
    6464       ENDDO
    6565

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/fxysinus.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE fxysinus (rlatu,yprimu,rlatv,yprimv,rlatu1,yprimu1,
    … …  
    3232c
    3333       DO j = 1, jjm + 1
    34           rlatu(j) = fy    ( FLOAT( j )        )
    35          yprimu(j) = fyprim( FLOAT( j )        )
     34          rlatu(j) = fy    ( REAL( j )        )
     35         yprimu(j) = fyprim( REAL( j )        )
    3636       ENDDO
    3737
    … …  
    3939       DO j = 1, jjm
    4040
    41          rlatv(j)  = fy    ( FLOAT( j ) + 0.5  )
    42          rlatu1(j) = fy    ( FLOAT( j ) + 0.25 )
    43          rlatu2(j) = fy    ( FLOAT( j ) + 0.75 )
     41         rlatv(j)  = fy    ( REAL( j ) + 0.5  )
     42         rlatu1(j) = fy    ( REAL( j ) + 0.25 )
     43         rlatu2(j) = fy    ( REAL( j ) + 0.75 )
    4444
    45         yprimv(j)  = fyprim( FLOAT( j ) + 0.5  )
    46         yprimu1(j) = fyprim( FLOAT( j ) + 0.25 )
    47         yprimu2(j) = fyprim( FLOAT( j ) + 0.75 )
     45        yprimv(j)  = fyprim( REAL( j ) + 0.5  )
     46        yprimu1(j) = fyprim( REAL( j ) + 0.25 )
     47        yprimu2(j) = fyprim( REAL( j ) + 0.75 )
    4848
    4949       ENDDO
    … …  
    5353c
    5454       DO i = 1, iim + 1
    55            rlonv(i)     = fx    (   FLOAT( i )          )
    56            rlonu(i)     = fx    (   FLOAT( i ) + 0.5    )
    57         rlonm025(i)     = fx    (   FLOAT( i ) - 0.25  )
    58         rlonp025(i)     = fx    (   FLOAT( i ) + 0.25  )
     55           rlonv(i)     = fx    (   REAL( i )          )
     56           rlonu(i)     = fx    (   REAL( i ) + 0.5    )
     57        rlonm025(i)     = fx    (   REAL( i ) - 0.25  )
     58        rlonp025(i)     = fx    (   REAL( i ) + 0.25  )
    5959
    60          xprimv  (i)    = fxprim (  FLOAT( i )          )
    61          xprimu  (i)    = fxprim (  FLOAT( i ) + 0.5    )
    62         xprimm025(i)    = fxprim (  FLOAT( i ) - 0.25   )
    63         xprimp025(i)    = fxprim (  FLOAT( i ) + 0.25   )
     60         xprimv  (i)    = fxprim (  REAL( i )          )
     61         xprimu  (i)    = fxprim (  REAL( i ) + 0.5    )
     62        xprimm025(i)    = fxprim (  REAL( i ) - 0.25   )
     63        xprimp025(i)    = fxprim (  REAL( i ) + 0.25   )
    6464       ENDDO
    6565

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/fyhyp.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    7575       depi     = 2. * pi
    7676       pis2     = pi/2.
    77        pisjm    = pi/ FLOAT(jjm)
     77       pisjm    = pi/ REAL(jjm)
    7878       epsilon  = 1.e-3
    7979       y0       =  yzoomdeg * pi/180.
    … …  
    9494
    9595       DO i = 0, nmax2
    96         yt(i) = - pis2  + FLOAT(i)* pi /nmax2
     96        yt(i) = - pis2  + REAL(i)* pi /nmax2
    9797       ENDDO
    9898
    … …  
    210210       DO 1500 j =  1,jlat
    211211        yo1   = 0.
    212         ylon2 =  - pis2 + pisjm * ( FLOAT(j)  + yuv  -1.) 
     212        ylon2 =  - pis2 + pisjm * ( REAL(j)  + yuv  -1.) 
    213213        yfi    = ylon2
    214214c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/gcm.F

    r1279 r1299  
    1818      USE getparam
    1919      USE filtreg_mod
     20      USE control_mod
    2021
    2122! Ehouarn: for now these only apply to Earth:
    … …  
    6667#include "logic.h"
    6768#include "temps.h"
    68 #include "control.h"
    6969#include "ener.h"
    7070#include "description.h"
    … …  
    308308      ENDIF
    309309
    310       zdtvr    = daysec/FLOAT(day_step)
     310      zdtvr    = daysec/REAL(day_step)
    311311        IF(dtvr.NE.zdtvr) THEN
    312312         WRITE(lunout,*)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/grid_atob.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE grille_m(imdep, jmdep, xdata, ydata, entree,
    … …  
    717717c Calculs intermediares:
    718718c
    719       xtmp(1) = -180.0 + 360.0/FLOAT(imtmp) / 2.0
     719      xtmp(1) = -180.0 + 360.0/REAL(imtmp) / 2.0
    720720      DO i = 2, imtmp
    721          xtmp(i) = xtmp(i-1) + 360.0/FLOAT(imtmp)
     721         xtmp(i) = xtmp(i-1) + 360.0/REAL(imtmp)
    722722      ENDDO
    723723      DO i = 1, imtmp
    724724         xtmp(i) = xtmp(i) /180.0 * 4.0*ATAN(1.0)
    725725      ENDDO
    726       ytmp(1) = -90.0 + 180.0/FLOAT(jmtmp) / 2.0
     726      ytmp(1) = -90.0 + 180.0/REAL(jmtmp) / 2.0
    727727      DO j = 2, jmtmp
    728          ytmp(j) = ytmp(j-1) + 180.0/FLOAT(jmtmp)
     728         ytmp(j) = ytmp(j-1) + 180.0/REAL(jmtmp)
    729729      ENDDO
    730730      DO j = 1, jmtmp

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/grid_noro.F

    r764 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    9393      xpi=acos(-1.)
    9494      rad    = 6 371 229.
    95       zdeltay=2.*xpi/float(jusn)*rad
     95      zdeltay=2.*xpi/REAL(jusn)*rad
    9696c
    9797c utilise-t'on un masque lu?
    … …  
    215215c  SUMMATION OVER GRIDPOINT AREA
    216216c
    217       zleny=xpi/float(jusn)*rad
    218       xincr=xpi/2./float(jusn)
     217      zleny=xpi/REAL(jusn)*rad
     218      xincr=xpi/2./REAL(jusn)
    219219       DO ii = 1, imar+1
    220220       DO jj = 1, jmar
    … …  
    468468      DO IS=-1,1
    469469        DO JS=-1,1
    470           WEIGHTpb(IS,JS)=1./FLOAT((1+IS**2)*(1+JS**2))
     470          WEIGHTpb(IS,JS)=1./REAL((1+IS**2)*(1+JS**2))
    471471          SUM=SUM+WEIGHTpb(IS,JS)
    472472        ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/grilles_gcm_netcdf.F

    r764 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    218218      open (20,file='grille.dat',form='unformatted',access='direct'
    219219     s      ,recl=4*ip1jmp1)
    220       write(20,rec=1) ((float(mod(i,2)+mod(j,2)),i=1,iip1),j=1,jjp1)
    221       write(20,rec=2) ((float(mod(i,2)*mod(j,2)),i=1,iip1),j=1,jjp1)
     220      write(20,rec=1) (( REAL(mod(i,2)+mod(j,2)),i=1,iip1),j=1,jjp1)
     221      write(20,rec=2) (( REAL(mod(i,2)*mod(j,2)),i=1,iip1),j=1,jjp1)
    222222      do j=2,jjm
    223223         dlat1(j)=180.*(rlatv(j)-rlatv(j-1))/pi
    224 c        dlat2(j)=180.*fyprim(float(j))/pi
     224c        dlat2(j)=180.*fyprim( REAL(j))/pi
    225225      enddo
    226226      do i=2,iip1
    227227         dlon1(i)=180.*(rlonu(i)-rlonu(i-1))/pi
    228 c        dlon2(i)=180.*fxprim(float(i))/pi
     228c        dlon2(i)=180.*fxprim( REAL(i))/pi
    229229      enddo
    230230      do j=2,jjm

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/guide_p_mod.F90

    r1279 r1299  
    11!
    2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/guide.F,v 1.3.4.1 2006/11/06 15:51:16 fairhead Exp $
     2! $Id$
    33!
    44MODULE guide_p_mod
    … …  
    6666  SUBROUTINE guide_init
    6767
     68    USE control_mod
    6869    IMPLICIT NONE
    6970 
    … …  
    7172    INCLUDE "paramet.h"
    7273    INCLUDE "netcdf.inc"
    73     INCLUDE "control.h"
    7474
    7575    INTEGER                :: error,ncidpl,rid,rcod
    … …  
    274274  SUBROUTINE guide_main(itau,ucov,vcov,teta,q,masse,ps)
    275275    use parallel
     276    USE control_mod
    276277   
    277278    IMPLICIT NONE
    … …  
    279280    INCLUDE "dimensions.h"
    280281    INCLUDE "paramet.h"
    281     INCLUDE "control.h"
    282282    INCLUDE "comconst.h"
    283283    INCLUDE "comvert.h"
    … …  
    380380      dday_step=real(day_step)
    381381      IF (iguide_read.LT.0) THEN
    382           tau=ditau/dday_step/FLOAT(iguide_read)
     382          tau=ditau/dday_step/ REAL(iguide_read)
    383383      ELSE
    384           tau=FLOAT(iguide_read)*ditau/dday_step
     384          tau= REAL(iguide_read)*ditau/dday_step
    385385      ENDIF
    386386      reste=tau-AINT(tau)
    … …  
    580580              ENDDO
    581581          ENDDO
    582           fieldm(:,l)=fieldm(:,l)/FLOAT(imax(typ)-imin(typ)+1)
     582          fieldm(:,l)=fieldm(:,l)/ REAL(imax(typ)-imin(typ)+1)
    583583    ! Compute forcing
    584584          DO j=jjb_v,jje_v
    … …  
    598598              ENDDO
    599599          ENDDO
    600           fieldm(:,l)=fieldm(:,l)/FLOAT(imax(typ)-imin(typ)+1)
     600          fieldm(:,l)=fieldm(:,l)/ REAL(imax(typ)-imin(typ)+1)
    601601    ! Compute forcing
    602602          DO j=jjb_u,jje_u

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/infotrac.F90

    r1279 r1299  
    3131
    3232  SUBROUTINE infotrac_init
     33    USE control_mod
    3334    IMPLICIT NONE
    3435!=======================================================================
    … …  
    4950
    5051    INCLUDE "dimensions.h"
    51     INCLUDE "control.h"
    5252    INCLUDE "iniprint.h"
    5353

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/iniacademic.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    3131c
    3232c=======================================================================
     33      USE control_mod
    3334      IMPLICIT NONE
    3435c-----------------------------------------------------------------------
    … …  
    4445#include "ener.h"
    4546#include "temps.h"
    46 #include "control.h"
    4747#include "iniprint.h"
    4848
    … …  
    9393        g      = 9.8
    9494        daysec = 86400.
    95         dtvr    = daysec/FLOAT(day_step)
     95        dtvr    = daysec/REAL(day_step)
    9696        zdtvr=dtvr
    9797        kappa  = 0.2857143

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/iniconst.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE iniconst
    55
     6      USE control_mod
    67      IMPLICIT NONE
    78c
    … …  
    1617#include "comconst.h"
    1718#include "temps.h"
    18 #include "control.h"
    1919#include "comvert.h"
    2020

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/inidissip.F

    r1279 r1299  
    1111c   -------------
    1212
     13      USE control_mod
     14
    1315      IMPLICIT NONE
    1416#include "dimensions.h"
    … …  
    1719#include "comconst.h"
    1820#include "comvert.h"
    19 #include "control.h"
    2021#include "logic.h"
    2122
    … …  
    165166
    166167c     IF(.NOT.lstardis) THEN
    167          fact    = rad*24./float(jjm)
     168         fact    = rad*24./REAL(jjm)
    168169         fact    = fact*fact
    169170         PRINT*,'coef u ', fact/cdivu, 1./cdivu

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/inigeom.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    168168c
    169169      IF( nitergdiv.NE.2 ) THEN
    170         gamdi_gdiv = coefdis/ ( float(nitergdiv) -2. )
     170        gamdi_gdiv = coefdis/ ( REAL(nitergdiv) -2. )
    171171      ELSE
    172172        gamdi_gdiv = 0.
    173173      ENDIF
    174174      IF( nitergrot.NE.2 ) THEN
    175         gamdi_grot = coefdis/ ( float(nitergrot) -2. )
     175        gamdi_grot = coefdis/ ( REAL(nitergrot) -2. )
    176176      ELSE
    177177        gamdi_grot = 0.
    178178      ENDIF
    179179      IF( niterh.NE.2 ) THEN
    180         gamdi_h = coefdis/ ( float(niterh) -2. )
     180        gamdi_h = coefdis/ ( REAL(niterh) -2. )
    181181      ELSE
    182182        gamdi_h = 0.
    … …  
    381381       yprp               = yprimu2(j-1)
    382382       rlatp              = rlatu2 (j-1)
    383 ccc       yprp             = fyprim( FLOAT(j) - 0.25 )
    384 ccc       rlatp            = fy    ( FLOAT(j) - 0.25 )
     383ccc       yprp             = fyprim( REAL(j) - 0.25 )
     384ccc       rlatp            = fy    ( REAL(j) - 0.25 )
    385385c
    386386      coslatp             = COS( rlatp )
    … …  
    416416        rlatm    = rlatu1 (  j  )
    417417        yprm     = yprimu1(  j  )
    418 cc         rlatp    = fy    ( FLOAT(j) - 0.25 )
    419 cc         yprp     = fyprim( FLOAT(j) - 0.25 )
    420 cc         rlatm    = fy    ( FLOAT(j) + 0.25 )
    421 cc         yprm     = fyprim( FLOAT(j) + 0.25 )
     418cc         rlatp    = fy    ( REAL(j) - 0.25 )
     419cc         yprp     = fyprim( REAL(j) - 0.25 )
     420cc         rlatm    = fy    ( REAL(j) + 0.25 )
     421cc         yprm     = fyprim( REAL(j) + 0.25 )
    422422
    423423         coslatm  = COS( rlatm )

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/integrd_p.F

    r1279 r1299  
    66     $     dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps0,masse,phis,finvmaold)
    77      USE parallel
     8      USE control_mod
    89      IMPLICIT NONE
    910
    … …  
    3233#include "temps.h"
    3334#include "serre.h"
    34 #include "control.h"
    3535
    3636c   Arguments:

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/interpre.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44       subroutine interpre(q,qppm,w,fluxwppm,masse,
    … …  
    66     s            unatppm,vnatppm,psppm)
    77
    8        implicit none
     8      USE control_mod
     9      implicit none
    910
    1011#include "dimensions.h"
    … …  
    1718#include "logic.h"
    1819#include "temps.h"
    19 #include "control.h"
    2020#include "ener.h"
    2121#include "description.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/leapfrog_p.F

    r1286 r1299  
    2020       USE guide_p_mod, ONLY : guide_main
    2121       USE getparam
     22       USE control_mod
    2223
    2324      IMPLICIT NONE
    … …  
    6263#include "logic.h"
    6364#include "temps.h"
    64 #include "control.h"
    6565#include "ener.h"
    6666#include "description.h"
    … …  
    212212      itau = 0
    213213!      iday = day_ini+itau/day_step
    214 !      time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
     214!      time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
    215215!         IF(time.GT.1.) THEN
    216216!          time = time-1.
    … …  
    12801280              itau= itau + 1
    12811281!              iday= day_ini+itau/day_step
    1282 !              time= FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
     1282!              time= REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
    12831283!                IF(time.GT.1.) THEN
    12841284!                  time = time-1.
    … …  
    14581458             itau =  itau + 1
    14591459!             iday = day_ini+itau/day_step
    1460 !             time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
     1460!             time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
    14611461!
    14621462!                  IF(time.GT.1.) THEN

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/limit_netcdf.F

    r1293 r1299  
    88! This routine is designed to work for Earth
    99      USE dimphy
     10      USE control_mod
    1011      use phys_state_var_mod , ONLY : pctsrf
    1112      use inter_barxy_m, only: inter_barxy
    … …  
    2930#include "dimensions.h"
    3031#include "paramet.h"
    31 #include "control.h"
    3232#include "logic.h"
    3333#include "netcdf.inc"
    … …  
    148148      kappa  = 0.2857143
    149149      cpp    = 1004.70885
    150       dtvr    = daysec/FLOAT(day_step)
     150      dtvr    = daysec/ REAL(day_step)
    151151      CALL inigeom
    152152c
    … …  
    418418          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    419419          DO k = 1, 360
    420             time = FLOAT(k-1)
     420            time =  REAL(k-1)
    421421            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    422422            champan(i,j,k) = by
    … …  
    632632          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    633633          DO k = 1, 360
    634             time = FLOAT(k-1)
     634            time =  REAL(k-1)
    635635            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    636636            champan(i,j,k) = by
    … …  
    951951          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    952952          DO k = 1, 360
    953             time = FLOAT(k-1)
     953            time =  REAL(k-1)
    954954            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    955955            champan(i,j,k) = by
    … …  
    11461146          CALL SPLINE(ax,ay,lmdep,1.e30,1.e30,yder)
    11471147          DO k = 1, 360
    1148             time = FLOAT(k-1)
     1148            time =  REAL(k-1)
    11491149            CALL SPLINT(ax,ay,yder,lmdep,time,by)
    11501150            champan(i,j,k) = by
    … …  
    13001300      ierr = NF_PUT_VARA_DOUBLE (nid,id_RUG,debut,epais,phy_rug(1,k))
    13011301#else
    1302       ierr = NF_PUT_VAR1_REAL (nid,id_tim,k,FLOAT(k))
     1302      ierr = NF_PUT_VAR1_REAL (nid,id_tim,k, REAL(k))
    13031303      IF (newlmt ) THEN
    13041304          ierr = NF_PUT_VARA_REAL (nid,id_FOCE,debut,epais

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/ppm3d.F

    r764 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44
    … …  
    345345C
    346346      PI = 4. * ATAN(1.)
    347       DL = 2.*PI / float(IMR)
    348       DP =    PI / float(JMR)
     347      DL = 2.*PI / REAL(IMR)
     348      DP =    PI / REAL(JMR)
    349349C
    350350      if(IGD.eq.0) then
    … …  
    388388      ZTC  = acos(CR1) * (180./PI)
    389389C
    390       JS0 = float(JMR)*(90.-ZTC)/180. + 2
     390      JS0 = REAL(JMR)*(90.-ZTC)/180. + 2
    391391      JS0 = max(JS0, J1+1)
    392392      IML = min(6*JS0/(J1-1)+2, 4*IMR/5)
    … …  
    628628C Contribution from the N-S advection
    629629      do i=1,imr*(j2-j1+1)
    630       JT = float(J1) - VA(i,j1)
     630      JT = REAL(J1) - VA(i,j1)
    631631      wk1(i,j1,2) = VA(i,j1) * (q(i,jt,k,IC) - q(i,jt+1,k,IC))
    632632      enddo
    … …  
    949949      IF(IORD.eq.1 .or. j.eq.j1. or. j.eq.j2) THEN
    950950      DO 1406 i=1,IMR
    951       iu = float(i) - uc(i,j)
     951      iu = REAL(i) - uc(i,j)
    9529521406  fx1(i) = qtmp(iu)
    953953      ELSE
    … …  
    957957      if(IORD.eq.2 .or. j.le.j1vl .or. j.ge.j2vl) then
    958958      DO 1408 i=1,IMR
    959       iu = float(i) - uc(i,j)
     959      iu = REAL(i) - uc(i,j)
    9609601408  fx1(i) = qtmp(iu) + DC(iu)*(sign(1.,uc(i,j))-uc(i,j))
    961961      else
    … …  
    11111111      if(JORD.eq.1) then
    11121112      DO 1000 i=1,len
    1113       JT = float(J1) - VC(i,J1)
     1113      JT = REAL(J1) - VC(i,J1)
    111411141000  fx(i,j1) = p(i,JT)
    11151115      else
    … …  
    11231123      else
    11241124      DO 1200 i=1,len
    1125       JT = float(J1) - VC(i,J1)
     1125      JT = REAL(J1) - VC(i,J1)
    112611261200  fx(i,j1) = p(i,JT) + (sign(1.,VC(i,j1))-VC(i,j1))*DC2(i,JT)
    11271127      endif
    … …  
    13581358        do j=j1-1,j2+1
    13591359      do i=1,imr
    1360       JP = float(j)-VA(i,j)
     1360      JP = REAL(j)-VA(i,j)
    13611361      ady(i,j) = VA(i,j)*(wk(i,jp)-wk(i,jp+1))
    13621362      enddo
    … …  
    15821582      JMR = JNP-1
    15831583      do 55 j=2,JNP
    1584         ph5  =  -0.5*PI + (FLOAT(J-1)-0.5)*DP
     1584        ph5  =  -0.5*PI + (REAL(J-1)-0.5)*DP
    1585158555      cose(j) = cos(ph5)
    15861586C
    … …  
    18341834C
    18351835c      if(first) then
    1836       DP = 4.*ATAN(1.)/float(JNP-1)
     1836      DP = 4.*ATAN(1.)/REAL(JNP-1)
    18371837      CAP1 = IMR*(1.-COS((j1-1.5)*DP))/DP
    18381838c      first = .false.
    … …  
    18891889C Check Poles.
    18901890      if(q(1,1).lt.0.) then
    1891       dq = q(1,1)*cap1/float(IMR)*acosp(j1)
     1891      dq = q(1,1)*cap1/REAL(IMR)*acosp(j1)
    18921892      do i=1,imr
    18931893      q(i,1) = 0.
    … …  
    18981898C
    18991899      if(q(1,JNP).lt.0.) then
    1900       dq = q(1,JNP)*cap1/float(IMR)*acosp(j2)
     1900      dq = q(1,JNP)*cap1/REAL(IMR)*acosp(j2)
    19011901      do i=1,imr
    19021902      q(i,JNP) = 0.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/ran1.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      FUNCTION RAN1(IDUM)
    … …  
    2020          IX1=MOD(IA1*IX1+IC1,M1)
    2121          IX2=MOD(IA2*IX2+IC2,M2)
    22           R(J)=(FLOAT(IX1)+FLOAT(IX2)*RM2)*RM1
     22          R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
    232311      CONTINUE
    2424        IDUM=1
    … …  
    3030      IF(J.GT.97.OR.J.LT.1)PAUSE
    3131      RAN1=R(J)
    32       R(J)=(FLOAT(IX1)+FLOAT(IX2)*RM2)*RM1
     32      R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
    3333      RETURN
    3434      END

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/sortvarc.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE sortvarc
    … …  
    5959
    6060       dtvrs1j   = dtvr/daysec
    61        rjour     = FLOAT( INT( itau * dtvrs1j ))
     61       rjour     = REAL( INT( itau * dtvrs1j ))
    6262       heure     = ( itau*dtvrs1j-rjour ) * 24.
    6363       imjmp1    = iim * jjp1
    … …  
    129129      ang   = SSUM(     llm,  angl, 1 )
    130130
    131 c      rday = FLOAT(INT ( day_ini + time ))
     131c      rday = REAL(INT ( day_ini + time ))
    132132c
    133        rday = FLOAT(INT(time-jD_ref-jH_ref))
     133       rday = REAL(INT(time-jD_ref-jH_ref))
    134134      IF(ptot0.eq.0.)  THEN
    135135         PRINT 3500, itau, rday, heure,time

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/sortvarc0.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE sortvarc0
    … …  
    6060
    6161       dtvrs1j   = dtvr/daysec
    62        rjour     = FLOAT( INT( itau * dtvrs1j ))
     62       rjour     = REAL( INT( itau * dtvrs1j ))
    6363       heure     = ( itau*dtvrs1j-rjour ) * 24.
    6464       imjmp1    = iim * jjp1
    … …  
    130130      ang0   = SSUM(     llm,  angl, 1 )
    131131
    132       rday = FLOAT(INT (time ))
     132      rday = REAL(INT (time ))
    133133c
    134134      PRINT 3500, itau, rday, heure, time

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/startvar.F

    r1293 r1299  
    352352      phis = phis * 9.81
    353353    !
    354 !      masque(:,:) = FLOAT(tmp_int(:,:))
     354!      masque(:,:) = REAL(tmp_int(:,:))
    355355    !
    356356    !  Compute surface roughness

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/tourabs.F

    r763 r1299  
    5757        ELSE
    5858         rot( ij,l ) = (vcov(ij+1,l)/cv(ij+1)-vcov(ij,l)/cv(ij))/
    59      $                 (2.*pi*RAD*cos(rlatv(j)))*float(iim)
     59     $                 (2.*pi*RAD*cos(rlatv(j)))*REAL(iim)
    6060     $                +(ucov(ij+iip1,l)/cu(ij+iip1)-ucov(ij,l)/cu(ij))/
    61      $                 (pi*RAD)*(float(jjm)-1.)
     61     $                 (pi*RAD)*(REAL(jjm)-1.)
    6262c
    6363        ENDIF

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/traceurpole.F

    r774 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44          subroutine traceurpole(q,masse)
     5
     6      USE control_mod
    57
    68          implicit none
    … …  
    1517#include "logic.h"
    1618#include "temps.h"
    17 #include "control.h"
    1819#include "ener.h"
    1920#include "description.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/ugeostr.F

    r1279 r1299  
    4040            DO i=1,iim
    4141               u(i,j,l)=fact*(phi(i,j+1,l)-phi(i,j,l))
    42                um(j,l)=um(j,l)+u(i,j,l)/float(iim)
     42               um(j,l)=um(j,l)+u(i,j,l)/REAL(iim)
    4343            ENDDO
    4444         ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/grid/fxy_new.h

    r524 r1299  
    88c....stretching in x...
    99c
    10         ripx(  ri )= (ri-1.0) *2.*pi/FLOAT(iim)
     10        ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim)
    1111        fx  (  ri )= ripx(ri) + transx  +
    1212     *         alphax * SIN( ripx(ri)+transx-pxo ) - pi
    13         fxprim(ri) = 2.*pi/FLOAT(iim)  *
     13        fxprim(ri) = 2.*pi/REAL(iim)  *
    1414     *        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
    1515
    1616c....stretching in y...
    1717c
    18         bigy(rj)   = 2.* (FLOAT(jjp1)-rj ) *pi/jjm
     18        bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
    1919        fy(rj)     =  ( bigy(rj) + transy  +
    2020     *        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/grid/fxy_reg.h

    r524 r1299  
    1313c
    1414c
    15       fy    ( rj ) =    pi/FLOAT(jjm) * ( 0.5 * FLOAT(jjm) +  1. - rj  )
    16       fyprim( rj ) =    pi/FLOAT(jjm)
     15      fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
     16      fyprim( rj ) =    pi/REAL(jjm)
    1717
    18 c     fy(rj)=ASIN(1.+2.*((1.-rj)/FLOAT(jjm)))
     18c     fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
    1919c     fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
    2020
    21       fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5*  FLOAT(iim) - 1. )
    22 c     fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5* ( FLOAT(iim) + 1.) )
    23       fxprim( ri ) = 2.*pi/FLOAT(iim)
     21      fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
     22c     fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
     23      fxprim( ri ) = 2.*pi/REAL(iim)
    2424c
    2525c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/grid/fxy_sin.h

    r524 r1299  
    1313c
    1414c
    15       fy(rj)=ASIN(1.+2.*((1.-rj)/FLOAT(jjm)))
     15      fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
    1616      fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
    1717
    18       fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5*  FLOAT(iim) - 1. )
    19 c     fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5* ( FLOAT(iim) + 1.) )
    20       fxprim( ri ) = 2.*pi/FLOAT(iim)
     18      fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
     19c     fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
     20      fxprim( ri ) = 2.*pi/REAL(iim)
    2121c
    2222c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/grid/fxyprim.h

    r524 r1299  
    1313c
    1414c
    15       fy    ( rj ) =    pi/FLOAT(jjm) * ( 0.5 * FLOAT(jjm) +  1. - rj  )
    16       fyprim( rj ) =    pi/FLOAT(jjm)
     15      fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
     16      fyprim( rj ) =    pi/REAL(jjm)
    1717
    18 c     fy(rj)=ASIN(1.+2.*((1.-rj)/FLOAT(jjm)))
     18c     fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
    1919c     fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
    2020
    21       fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5*  FLOAT(iim) - 1. )
    22 c     fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5* ( FLOAT(iim) + 1.) )
    23       fxprim( ri ) = 2.*pi/FLOAT(iim)
     21      fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
     22c     fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
     23      fxprim( ri ) = 2.*pi/REAL(iim)
    2424c
    2525c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/aaam_bud.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      subroutine aaam_bud (iam,nlon,nlev,rjour,rsec,
    … …  
    117117      REAL BLSU(801,401),BLSV(801,401)
    118118      REAL ZLON(801),ZLAT(401)
     119
     120      CHARACTER (LEN=20) :: modname='aaam_bud'
     121      CHARACTER (LEN=80) :: abort_message
     122
     123
    119124C
    120125C  PUT AAM QUANTITIES AT ZERO:
    121126C
    122127      if(iim+1.gt.801.or.jjm+1.gt.401)then
    123       print *,' Pb de dimension dans aaam_bud'
    124       stop
     128        abort_message = 'Pb de dimension dans aaam_bud'
     129        CALL abort_gcm (modname,abort_message,1)
    125130      endif
    126131
    … …  
    128133      hadley=1.e18
    129134      hadday=1.e18*24.*3600.
    130       dlat=xpi/float(jjm)
    131       dlon=2.*xpi/float(iim)
     135      dlat=xpi/REAL(jjm)
     136      dlon=2.*xpi/REAL(iim)
    132137     
    133138      do iax=1,3

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/aeropt.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE aeropt(pplay, paprs, t_seri, msulfate, RHcl,
    … …  
    3939      REAL alpha_aer_sulfate(nbre_RH,5)   !--unit m2/g SO4
    4040      REAL alphasulfate     
     41
     42      CHARACTER (LEN=20) :: modname='aeropt'
     43      CHARACTER (LEN=80) :: abort_message
     44
    4145c
    4246c Proprietes optiques
    … …  
    8589        rh=MIN(RHcl(i,k)*100.,RH_MAX)
    8690        RH_num = INT( rh/10. + 1.)
    87         IF (rh.LT.0.) STOP 'aeropt: RH < 0 not possible'
     91        IF (rh.LT.0.) THEN
     92          abort_message = 'aeropt: RH < 0 not possible'
     93          CALL abort_gcm (modname,abort_message,1)
     94        ENDIF
    8895        IF (rh.gt.85.) RH_num=10
    8996        IF (rh.gt.90.) RH_num=11

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/aeropt_5wv.F90

    r1279 r1299  
    621621  DO k=1, klev
    622622    DO i=1, klon
    623 !      IF (t_seri(i,k).EQ.0) stop 'stop aeropt_5wv T '
    624 !      IF (pplay(i,k).EQ.0) stop  'stop aeropt_5wv p '
    625623      zrho=pplay(i,k)/t_seri(i,k)/RD                  ! kg/m3
    626624!CDIR UNROLL=naero_spc

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/albedo.F

    r900 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    6767c prend en compte l'autre moitie de la journee):
    6868      DO k = 1, npts
    69          rmu = aa + bb * COS(FLOAT(k)/FLOAT(npts)*zpi)
     69         rmu = aa + bb * COS(REAL(k)/REAL(npts)*zpi)
    7070         rmu = MAX(0.0, rmu)
    7171         fauxo = (1.47-ACOS(rmu))/.15
    … …  
    110110c prend en compte l'autre moitie de la journee):
    111111      DO k = 1, npts
    112          rmu = aa + bb * COS(FLOAT(k)/FLOAT(npts)*zpi)
     112         rmu = aa + bb * COS(REAL(k)/REAL(npts)*zpi)
    113113         rmu = MAX(0.0, rmu)
    114114cIM cf. PB      alb = 0.058/(rmu + 0.30)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/calcul_simulISCCP.h

    r1279 r1299  
    11c
    2 c $Header$
     2c $Id$
    33c
    44c on appelle le simulateur ISCCP toutes les 3h
    … …  
    1818       sunlit(i)=1
    1919       IF(rmu0(i).EQ.0.) sunlit(i)=0
    20        nbsunlit(1,i,n)=FLOAT(sunlit(i))
     20       nbsunlit(1,i,n)=REAL(sunlit(i))
    2121      ENDDO
    2222c
    … …  
    8888           print*,'seed=0 i paprs aa seed_re',
    8989     .     i,paprs(i,2),aa,seed_re(i,n)
    90            STOP
     90           abort_message = ''
     91           CALL abort_gcm (modname,abort_message,1)
    9192          ELSE IF(seed(i,n).LT.0) THEN
    9293           print*,'seed < 0, i seed itap paprs',i,
    9394     .     seed(i,n),itap,paprs(i,2)
    94            STOP
     95           abort_message = ''
     96           CALL abort_gcm (modname,abort_message,1)
    9597          ENDIF
    9698c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/calltherm.F90

    r1295 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      subroutine calltherm(dtime  &
    … …  
    8383!      save zentr_therm,zfm_therm
    8484
     85      character (len=20) :: modname='calltherm'
     86      character (len=80) :: abort_message
     87
    8588      integer i,k
    8689      logical, save :: first=.true.
    … …  
    137140         if(nbptspb.GT.0) print*,'Number of points with q_seri(i,k)<=0 ',nbptspb   
    138141
    139          zdt=dtime/float(nsplit_thermals)
     142         zdt=dtime/REAL(nsplit_thermals)
    140143         do isplit=1,nsplit_thermals
    141144
    … …  
    173176     &      ,tau_thermals,3)
    174177          else if (iflag_thermals.eq.11) then
    175             stop 'cas non prevu dans calltherm'
     178              abort_message = 'cas non prevu dans calltherm'
     179              CALL abort_gcm (modname,abort_message,1)
     180
    176181!           CALL thermcell_pluie(klon,klev,zdt  &
    177182!   &      ,pplay,paprs,pphi,zlev  &
    … …  
    210215           if (prt_level.gt.10) write(lunout,*)'Apres thermcell_main OK'
    211216         else
    212             STOP'Cas des thermiques non prevu'
     217           abort_message = 'Cas des thermiques non prevu'
     218           CALL abort_gcm (modname,abort_message,1)
    213219         endif
    214220
    … …  
    218224      DO i=1,klon
    219225       logexpr1(i)=flag_bidouille_stratocu.or.weak_inversion(i).gt.0.5
    220        IF(logexpr1(i)) fact(i)=1./float(nsplit_thermals)
     226       IF(logexpr1(i)) fact(i)=1./REAL(nsplit_thermals)
    221227      ENDDO
    222228
    … …  
    255261        if(prt_level.GE.10) print*,'calltherm i Alp_bl Alp Ale_bl Ale',i,Alp_bl(i),Alp(i),Ale_bl(i),Ale(i)
    256262            fm_therm(i,klev+1)=0.
    257             Ale_bl(i)=Ale_bl(i)+Ale(i)/float(nsplit_thermals)
     263            Ale_bl(i)=Ale_bl(i)+Ale(i)/REAL(nsplit_thermals)
    258264!            write(22,*)'ALE CALLTHERM',Ale_bl(i),Ale(i)
    259             Alp_bl(i)=Alp_bl(i)+Alp(i)/float(nsplit_thermals)
     265            Alp_bl(i)=Alp_bl(i)+Alp(i)/REAL(nsplit_thermals)
    260266!            write(23,*)'ALP CALLTHERM',Alp_bl(i),Alp(i)
    261267       ENDDO
    … …  
    276282!    &         'fm=',zfm_therm(i,k),'entr=',entr_therm(i,k)
    277283                 endif
    278 !       stop
    279284            ENDDO
    280285            ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/conema3.F

    r1146 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE conema3 (dtime,paprs,pplay,t,q,u,v,tra,ntra,
    … …  
    360360      cape(i) = em_CAPE
    361361      wd(i) = em_wd
    362       rflag(i) = float(iflag)
     362      rflag(i) = REAL(iflag)
    363363c SB      kbas(i) = em_bas
    364364c SB      ktop(i) = em_top

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/conf_phys.F90

    r1286 r1299  
    2727   USE surface_data
    2828   USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl
     29   use control_mod
    2930
    3031 include "conema3.h"
    … …  
    3637include "clesphys.h"
    3738include "compbl.h"
    38 include "control.h"
    3939include "comsoil.h"
    4040!

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/convect2.F

    r766 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      subroutine convect2(ncum,idcum,len,nd,ndp1,nl,minorig,

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cpl_mod.F90

    r1279 r1299  
    2424  USE oasis
    2525  USE write_field_phy
     26  USE control_mod
     27
    2628 
    2729! Global attributes
    … …  
    101103    INCLUDE "dimensions.h"
    102104    INCLUDE "indicesol.h"
    103     INCLUDE "control.h"
    104105    INCLUDE "temps.h"
    105106    INCLUDE "iniprint.h"
    … …  
    583584    DO ig = 1, knon
    584585       cpl_sols(ig,cpl_index) = cpl_sols(ig,cpl_index) + &
    585             swdown(ig)      / FLOAT(nexca)
     586            swdown(ig)      / REAL(nexca)
    586587       cpl_nsol(ig,cpl_index) = cpl_nsol(ig,cpl_index) + &
    587             (lwdown(ig) + fluxlat(ig) +fluxsens(ig)) / FLOAT(nexca)
     588            (lwdown(ig) + fluxlat(ig) +fluxsens(ig)) / REAL(nexca)
    588589       cpl_rain(ig,cpl_index) = cpl_rain(ig,cpl_index) + &
    589             precip_rain(ig) / FLOAT(nexca)
     590            precip_rain(ig) / REAL(nexca)
    590591       cpl_snow(ig,cpl_index) = cpl_snow(ig,cpl_index) + &
    591             precip_snow(ig) / FLOAT(nexca)
     592            precip_snow(ig) / REAL(nexca)
    592593       cpl_evap(ig,cpl_index) = cpl_evap(ig,cpl_index) + &
    593             evap(ig)        / FLOAT(nexca)
     594            evap(ig)        / REAL(nexca)
    594595       cpl_tsol(ig,cpl_index) = cpl_tsol(ig,cpl_index) + &
    595             tsurf(ig)       / FLOAT(nexca)
     596            tsurf(ig)       / REAL(nexca)
    596597       cpl_fder(ig,cpl_index) = cpl_fder(ig,cpl_index) + &
    597             fder(ig)        / FLOAT(nexca)
     598            fder(ig)        / REAL(nexca)
    598599       cpl_albe(ig,cpl_index) = cpl_albe(ig,cpl_index) + &
    599             albsol(ig)      / FLOAT(nexca)
     600            albsol(ig)      / REAL(nexca)
    600601       cpl_taux(ig,cpl_index) = cpl_taux(ig,cpl_index) + &
    601             taux(ig)        / FLOAT(nexca)
     602            taux(ig)        / REAL(nexca)
    602603       cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) + &
    603             tauy(ig)        / FLOAT(nexca)     
     604            tauy(ig)        / REAL(nexca)     
    604605       cpl_windsp(ig,cpl_index) = cpl_windsp(ig,cpl_index) + &
    605             windsp(ig)      / FLOAT(nexca)
     606            windsp(ig)      / REAL(nexca)
    606607       cpl_taumod(ig,cpl_index) =   cpl_taumod(ig,cpl_index) + &
    607           SQRT ( taux(ig)*taux(ig)+tauy(ig)*tauy(ig) ) / FLOAT (nexca)
     608          SQRT ( taux(ig)*taux(ig)+tauy(ig)*tauy(ig) ) / REAL (nexca)
    608609
    609610       IF (carbon_cycle_cpl) THEN
    610611          cpl_atm_co2(ig,cpl_index) = cpl_atm_co2(ig,cpl_index) + &
    611                co2_send(knindex(ig))/ FLOAT(nexca)
     612               co2_send(knindex(ig))/ REAL(nexca)
    612613       END IF
    613614     ENDDO
    … …  
    777778    DO ig = 1, knon
    778779       cpl_sols(ig,cpl_index) = cpl_sols(ig,cpl_index) + &
    779             swdown(ig)      / FLOAT(nexca)
     780            swdown(ig)      / REAL(nexca)
    780781       cpl_nsol(ig,cpl_index) = cpl_nsol(ig,cpl_index) + &
    781             (lwdown(ig) + fluxlat(ig) +fluxsens(ig)) / FLOAT(nexca)
     782            (lwdown(ig) + fluxlat(ig) +fluxsens(ig)) / REAL(nexca)
    782783       cpl_rain(ig,cpl_index) = cpl_rain(ig,cpl_index) + &
    783             precip_rain(ig) / FLOAT(nexca)
     784            precip_rain(ig) / REAL(nexca)
    784785       cpl_snow(ig,cpl_index) = cpl_snow(ig,cpl_index) + &
    785             precip_snow(ig) / FLOAT(nexca)
     786            precip_snow(ig) / REAL(nexca)
    786787       cpl_evap(ig,cpl_index) = cpl_evap(ig,cpl_index) + &
    787             evap(ig)        / FLOAT(nexca)
     788            evap(ig)        / REAL(nexca)
    788789       cpl_tsol(ig,cpl_index) = cpl_tsol(ig,cpl_index) + &
    789             tsurf(ig)       / FLOAT(nexca)
     790            tsurf(ig)       / REAL(nexca)
    790791       cpl_fder(ig,cpl_index) = cpl_fder(ig,cpl_index) + &
    791             fder(ig)        / FLOAT(nexca)
     792            fder(ig)        / REAL(nexca)
    792793       cpl_albe(ig,cpl_index) = cpl_albe(ig,cpl_index) + &
    793             albsol(ig)      / FLOAT(nexca)
     794            albsol(ig)      / REAL(nexca)
    794795       cpl_taux(ig,cpl_index) = cpl_taux(ig,cpl_index) + &
    795             taux(ig)        / FLOAT(nexca)
     796            taux(ig)        / REAL(nexca)
    796797       cpl_tauy(ig,cpl_index) = cpl_tauy(ig,cpl_index) + &
    797             tauy(ig)        / FLOAT(nexca)     
     798            tauy(ig)        / REAL(nexca)     
    798799       cpl_taumod(ig,cpl_index) = cpl_taumod(ig,cpl_index) + &
    799             SQRT ( taux(ig)*taux(ig)+tauy(ig)*tauy(ig) ) / FLOAT(nexca)
     800            SQRT ( taux(ig)*taux(ig)+tauy(ig)*tauy(ig) ) / REAL(nexca)
    800801    ENDDO
    801802
    … …  
    944945!*************************************************************************************   
    945946!$OMP MASTER
    946     cpl_rriv2D(:,:) = cpl_rriv2D(:,:) + rriv2D(:,:) / FLOAT(nexca)
    947     cpl_rcoa2D(:,:) = cpl_rcoa2D(:,:) + rcoa2D(:,:) / FLOAT(nexca)
     947    cpl_rriv2D(:,:) = cpl_rriv2D(:,:) + rriv2D(:,:) / REAL(nexca)
     948    cpl_rcoa2D(:,:) = cpl_rcoa2D(:,:) + rcoa2D(:,:) / REAL(nexca)
    948949!$OMP END MASTER
    949950
    … …  
    998999!*************************************************************************************   
    9991000!$OMP MASTER
    1000     cpl_rlic2D(:,:) = cpl_rlic2D(:,:) + rlic2D(:,:) / FLOAT(nexca)
     1001    cpl_rlic2D(:,:) = cpl_rlic2D(:,:) + rlic2D(:,:) / REAL(nexca)
    10011002!$OMP END MASTER
    10021003

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv30_routines.F

    r879 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    749749
    750750#include "cv30param.h"
     751      include 'iniprint.h'
    751752
    752753c inputs:
    … …  
    778779c local variables:
    779780      integer i,k,nn,j
     781
     782      CHARACTER (LEN=20) :: modname='cv30_compress'
     783      CHARACTER (LEN=80) :: abort_message
    780784
    781785
    … …  
    820824
    821825      if (nn.ne.ncum) then
    822          print*,'strange! nn not equal to ncum: ',nn,ncum
    823          stop
     826         write(lunout,*)'strange! nn not equal to ncum: ',nn,ncum
     827         abort_message = ''
     828         CALL abort_gcm (modname,abort_message,1)
    824829      endif
    825830

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3_inicp.F

    r966 r1299  
    1313c
    1414      INTEGER iflag_clos
     15      CHARACTER (LEN=20) :: modname='cv3_inicp'
     16      CHARACTER (LEN=80) :: abort_message
    1517c
    1618c --   Mixing probability distribution functions
    … …  
    105107        if (abs(aire-1.0) .gt. 0.02) then
    106108            print *,'WARNING:: AREA OF MIXING PDF IS::', aire
    107             stop
     109            abort_message = ''
     110            CALL abort_gcm (modname,abort_message,1)
    108111        else
    109112            print *,'Area, mean & std deviation are ::', aire,mu,sigma

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3_inip.F

    r1146 r1299  
    1212c
    1313c      INTEGER iflag_mix
     14      include 'iniprint.h'
     15
     16      CHARACTER (LEN=20) :: modname='cv3_inip'
     17      CHARACTER (LEN=80) :: abort_message
     18
    1419c
    1520c --   Mixing probability distribution functions
    … …  
    104109c
    105110        if (abs(aire-1.0) .gt. 0.02) then
    106             print *,'WARNING:: AREA OF MIXING PDF IS::', aire
    107             stop
     111            write(lunout,*)'WARNING:: AREA OF MIXING PDF IS::', aire
     112            abort_message = ''
     113            CALL abort_gcm (modname,abort_message,1)
    108114        else
    109115            print *,'Area, mean & std deviation are ::', aire,mu,sigma

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3_routines.F

    r1277 r1299  
    11!
    2 ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/cv3_routines.F,v 1.16 2008-11-06 16:29:35 lmdzadmin Exp $
     2! $Id$
    33!
    44c
    … …  
    3535      integer nd
    3636      real delt ! timestep (seconds)
     37
     38      CHARACTER (LEN=20) :: modname='cv3_param'
     39      CHARACTER (LEN=80) :: abort_message
    3740
    3841c noff: integer limit for convection (nd-noff)
    … …  
    767770
    768771#include "cv3param.h"
     772      include 'iniprint.h'
    769773
    770774c inputs:
    … …  
    797801      integer i,k,nn,j
    798802
     803      CHARACTER (LEN=20) :: modname='cv3_compress'
     804      CHARACTER (LEN=80) :: abort_message
    799805
    800806      do 110 k=1,nl+1
    … …  
    839845
    840846      if (nn.ne.ncum) then
    841          print*,'strange! nn not equal to ncum: ',nn,ncum
    842          stop
     847         write(lunout,*)'strange! nn not equal to ncum: ',nn,ncum
     848         abort_message = ''
     849         CALL abort_gcm (modname,abort_message,1)
    843850      endif
    844851

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3a_compress.F

    r1146 r1299  
    7676      integer i,k,nn,j
    7777
     78      CHARACTER (LEN=20) :: modname='cv3a_compress'
     79      CHARACTER (LEN=80) :: abort_message
     80
    7881
    7982      do 110 k=1,nl+1
    … …  
    127130
    128131      if (nn.ne.ncum) then
    129          print*,'WARNING nn not equal to ncum: ',nn,ncum
    130          stop
     132        print*,'WARNING nn not equal to ncum: ',nn,ncum
     133        abort_message = ''
     134        CALL abort_gcm (modname,abort_message,1)
    131135      endif
    132136
    … …  
    157161      if (nn.ne.ncum) then
    158162         print*,'WARNING nn not equal to ncum: ',nn,ncum
    159          stop
     163         abort_message = ''
     164         CALL abort_gcm (modname,abort_message,1)
    160165      endif
    161166

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3p1_closure.F

    r973 r1299  
     1!
     2! $Id$
     3!
    14      SUBROUTINE cv3p1_closure(nloc,ncum,nd,icb,inb
    25     :                      ,pbase,plcl,p,ph,tv,tvp,buoy
    … …  
    7477      real wb,sigmax
    7578      data wb /2./, sigmax /0.1/
     79
     80      CHARACTER (LEN=20) :: modname='cv3p1_closure'
     81      CHARACTER (LEN=80) :: abort_message
    7682c
    7783c      print *,' -> cv3p1_closure, Ale ',ale(1)
    … …  
    509515       cbmf1(il) = alp2(il)/(0.5*wb*wb-Cin(il))
    510516       if(cbmf1(il).EQ.0.AND.alp2(il).NE.0.) THEN
    511         print*,'cv3p1_closure cbmf1=0 and alp NE 0 il alp2 alp cin ',il,
     517        write(lunout,*)
     518     &  'cv3p1_closure cbmf1=0 and alp NE 0 il alp2 alp cin ',il,
    512519     . alp2(il),alp(il),cin(il)
    513         STOP
     520        abort_message = ''
     521        CALL abort_gcm (modname,abort_message,1)
    514522       endif
    515523       cbmfmax(il) = sigmax*wb2(il)*100.*p(il,icb(il))

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv_routines.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE cv_param(nd)
    … …  
    3838#include "cvparam.h"
    3939      integer nd
     40      CHARACTER (LEN=20) :: modname='cv_routines'
     41      CHARACTER (LEN=80) :: abort_message
    4042
    4143c noff: integer limit for convection (nd-noff)
    … …  
    429431c local variables:
    430432      integer i,k,nn
     433      CHARACTER (LEN=20) :: modname='cv_compress'
     434      CHARACTER (LEN=80) :: abort_message
     435
     436      include 'iniprint.h'
    431437
    432438
    … …  
    456462
    457463      if (nn.ne.ncum) then
    458          print*,'strange! nn not equal to ncum: ',nn,ncum
    459          stop
     464         write(lunout,*)'strange! nn not equal to ncum: ',nn,ncum
     465         abort_message = ''
     466         CALL abort_gcm (modname,abort_message,1)
    460467      endif
    461468

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cva_driver.F

    r1279 r1299  
    106106#include "dimensions.h"
    107107ccccc#include "dimphy.h"
     108      include 'iniprint.h'
     109
    108110c
    109111c Input
    … …  
    419421      logical, save :: first=.true.
    420422c$OMP THREADPRIVATE(first)
     423      CHARACTER (LEN=20) :: modname='cva_driver'
     424      CHARACTER (LEN=80) :: abort_message
    421425
    422426c
    … …  
    563567c test niveaux couche alimentation KE
    564568       if(sig1feed1.eq.sig2feed1) then
    565                print*,'impossible de choisir sig1feed=sig2feed'
    566                print*,'changer la valeur de sig2feed dans physiq.def'
    567        stop
     569         write(lunout,*)'impossible de choisir sig1feed=sig2feed'
     570         write(lunout,*)'changer la valeur de sig2feed dans physiq.def'
     571         abort_message = ''
     572         CALL abort_gcm (modname,abort_message,1)
    568573       endif
    569574c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/fisrtilp.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    135135         PRINT*, 'fisrtilp, evap_prec:', evap_prec
    136136         PRINT*, 'fisrtilp, cpartiel:', cpartiel
    137          IF (ABS(dtime/FLOAT(ninter)-360.0).GT.0.001) THEN
     137         IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN
    138138          PRINT*, 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
    139139          PRINT*, 'Je prefere un sous-intervalle de 6 minutes'
    … …  
    436436         zfice(i) = zfice(i)**nexpo
    437437         zneb(i) = MAX(rneb(i,k), seuil_neb)
    438          radliq(i,k) = zoliq(i)/FLOAT(ninter+1)
     438         radliq(i,k) = zoliq(i)/REAL(ninter+1)
    439439      ENDIF
    440440      ENDDO
    … …  
    453453                zcl   =cld_lc_con
    454454                zct   =1./cld_tau_con
    455                 zfroi    = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
     455                zfroi    = dtime/REAL(ninter)/zdz(i)*zoliq(i)
    456456     .              *fallvc(zrhol(i)) * zfice(i)
    457457             else
    458458                zcl   =cld_lc_lsc
    459459                zct   =1./cld_tau_lsc
    460                 zfroi    = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
     460                zfroi    = dtime/REAL(ninter)/zdz(i)*zoliq(i)
    461461     .              *fallvs(zrhol(i)) * zfice(i)
    462462             endif
    463              zchau    = zct   *dtime/FLOAT(ninter) * zoliq(i)
     463             zchau    = zct   *dtime/REAL(ninter) * zoliq(i)
    464464     .         *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl   )**2)) *(1.-zfice(i))
    465465             ztot    = zchau    + zfroi
    … …  
    468468         ztot    = MIN(ztot,zoliq(i))
    469469         zoliq(i) = MAX(zoliq(i)-ztot   , 0.0)
    470          radliq(i,k) = radliq(i,k) + zoliq(i)/FLOAT(ninter+1)
     470         radliq(i,k) = radliq(i,k) + zoliq(i)/REAL(ninter+1)
    471471      ENDIF
    472472      ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/fisrtilp_tr.F

    r766 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    140140         PRINT*, 'fisrtilp, evap_prec:', evap_prec
    141141         PRINT*, 'fisrtilp, cpartiel:', cpartiel
    142          IF (ABS(dtime/FLOAT(ninter)-360.0).GT.0.001) THEN
     142         IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN
    143143          PRINT*, 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime
    144144          PRINT*, 'Je prefere un sous-intervalle de 6 minutes'
    … …  
    335335         zfice(i) = zfice(i)**nexpo
    336336         zneb(i) = MAX(rneb(i,k), seuil_neb)
    337          radliq(i,k) = zoliq(i)/FLOAT(ninter+1)
     337         radliq(i,k) = zoliq(i)/REAL(ninter+1)
    338338      ENDIF
    339339      ENDDO
    … …  
    342342      DO i = 1, klon
    343343      IF (rneb(i,k).GT.0.0) THEN
    344          zchau(i) = ct*dtime/FLOAT(ninter) * zoliq(i)
     344         zchau(i) = ct*dtime/REAL(ninter) * zoliq(i)
    345345     .          * (1.0-EXP(-(zoliq(i)/zneb(i)/cl)**2)) *(1.-zfice(i))
    346346         zrhol(i) = zrho(i) * zoliq(i) / zneb(i)
    347          zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
     347         zfroi(i) = dtime/REAL(ninter)/zdz(i)*zoliq(i)
    348348     .              *fallv(zrhol(i)) * zfice(i)
    349349         ztot(i) = zchau(i) + zfroi(i)
    … …  
    351351         ztot(i) = MIN(MAX(ztot(i),0.0),zoliq(i))
    352352         zoliq(i) = MAX(zoliq(i)-ztot(i), 0.0)
    353          radliq(i,k) = radliq(i,k) + zoliq(i)/FLOAT(ninter+1)
     353         radliq(i,k) = radliq(i,k) + zoliq(i)/REAL(ninter+1)
    354354      ENDIF
    355355      ENDDO

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/hines_gwd.F

    r1279 r1299  
    847847C  Use horizontal isotropy to calculate azimuthal variances at bottom level.
    848848C
    849       AZFAC = 1. / FLOAT(NAZ)
     849      AZFAC = 1. / REAL(NAZ)
    850850      DO 20 N = 1,NAZ
    851851        DO 10 I = IL1,IL2

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_bilKP_ave.h

    r766 r1299  
    11c
    2 c $Header$
     2c $Id$
    33c
    44      IF (ok_journe) THEN
    … …  
    1717cym         ENDDO
    1818         DO ll=1,klev
    19             znivsig(ll)=float(ll)
     19            znivsig(ll)=REAL(ll)
    2020         ENDDO
    2121cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_bilKP_ins.h

    r766 r1299  
    11c
    2 c $Header$
     2c $Id$
    33c
    44      IF (ok_journe) THEN
    … …  
    1717cym         ENDDO
    1818         DO ll=1,klev
    19             znivsig(ll)=float(ll)
     19            znivsig(ll)=REAL(ll)
    2020         ENDDO
    2121cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_histISCCP.h

    r1045 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      IF (ok_isccp) THEN
    … …  
    4949c
    5050        DO l=1, ncol(n)
    51           vertlev(l,n)=float(l)
     51          vertlev(l,n)=REAL(l)
    5252        ENDDO !ncol
    5353c

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_histday_seri.h

    r776 r1299  
    11c
    2 c $Header$
     2c $Id$
    33c
    44cym Ne fonctionnera pas en mode parallele
    … …  
    1919         ENDDO
    2020         DO ll=1,klev
    21             znivsig(ll)=float(ll)
     21            znivsig(ll)=REAL(ll)
    2222         ENDDO
    2323         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/ini_histmthNMC.h

    r776 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      IF (ok_mensuel) THEN
    … …  
    2020cym         ENDDO
    2121         DO ll=1,klev
    22             znivsig(ll)=float(ll)
     22            znivsig(ll)=REAL(ll)
    2323         ENDDO
    2424cym         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/inifis.F

    r987 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE inifis(ngrid,nlayer,
    … …  
    4545cym#include "dimphy.h"
    4646
     47      INCLUDE 'iniprint.h'
    4748      REAL prad,pg,pr,pcpp,punjours
    4849 
    … …  
    5253 
    5354      REAL ptimestep
     55      CHARACTER (LEN=20) :: modname='inifis'
     56      CHARACTER (LEN=80) :: abort_message
     57
    5458 
    5559      IF (nlayer.NE.klev) THEN
    … …  
    5862         PRINT*,'nlayer     = ',nlayer
    5963         PRINT*,'klev   = ',klev
    60          STOP
     64         abort_message = ''
     65         CALL abort_gcm (modname,abort_message,1)
    6166      ENDIF
    6267
    … …  
    6671         PRINT*,'ngrid     = ',ngrid
    6772         PRINT*,'klon   = ',klon
    68          STOP
     73         abort_message = ''
     74         CALL abort_gcm (modname,abort_message,1)
    6975      ENDIF
    7076
    7177      RETURN
    72 9999  STOP'Cette version demande les fichier rnatur.dat et surf.def'
     789999  continue
     79      abort_message = 'Cette version demande les fichier rnatur.dat
     80     & et surf.def'
     81      CALL abort_gcm (modname,abort_message,1)
     82
    7383      END

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/iniphysiq.F

    r879 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    6161 
    6262      REAL ptimestep
     63      CHARACTER (LEN=20) :: modname='iniphysiq'
     64      CHARACTER (LEN=80) :: abort_message
    6365 
    6466      IF (nlayer.NE.klev) THEN
    … …  
    6769         PRINT*,'nlayer     = ',nlayer
    6870         PRINT*,'klev   = ',klev
    69          STOP
     71         abort_message = ''
     72         CALL abort_gcm (modname,abort_message,1)
    7073      ENDIF
    7174
    … …  
    7578         PRINT*,'ngrid     = ',ngrid
    7679         PRINT*,'klon   = ',klon_glo
    77          STOP
     80         abort_message = ''
     81         CALL abort_gcm (modname,abort_message,1)
    7882      ENDIF
    7983c$OMP PARALLEL PRIVATE(ibegin,iend)
    … …  
    96100
    97101      RETURN
    98 9999  STOP'Cette version demande les fichier rnatur.dat et surf.def'
     1029999  CONTINUE
     103      abort_message ='Cette version demande les fichier rnatur.dat
     104     & et surf.def'
     105      CALL abort_gcm (modname,abort_message,1)
     106
    99107      END

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/initphysto.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44C
    … …  
    1212       USE IOIPSL
    1313       USE iophy
     14       USE control_mod
     15
    1416      implicit none
    1517
    … …  
    5254#include "serre.h"
    5355#include "indicesol.h"
    54 #include "control.h"
    5556cym#include "dimphy.h"
    5657
    … …  
    108109C
    109110        DO l=1,llm
    110             nivsigs(l)=float(l)
     111            nivsigs(l)=REAL(l)
    111112         ENDDO
    112113

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/initrrnpb.F90

    r1279 r1299  
    11!
    2 ! $Id $
     2! $Id$
    33!
    44SUBROUTINE  initrrnpb(ftsol,pctsrf,masktr,fshtr,hsoltr,tautr,vdeptr,scavtr)
    … …  
    3939  REAL                                  :: s
    4040
     41  CHARACTER (LEN=20) :: modname='initrrnpb'
     42  CHARACTER (LEN=80) :: abort_message
     43
     44
    4145  WRITE(*,*)'PASSAGE initrrnpb ...'
    4246!
    4347! Radon it = 1
    4448!----------------
    45   IF ( nbtr .LE. 0 ) STOP '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def'
     49  IF ( nbtr .LE. 0 ) then
     50    abort_message = '**PHYTRAC:initrrnpb:** nbtr < 0; verifier RN dans traceur.def'
     51    CALL abort_gcm (modname,abort_message,1)
     52  ENDIF
    4653  it = 1
    4754  s = 1.E4             ! Source: atome par m2
    … …  
    6875! 210Pb it = 2
    6976!----------------
    70   IF ( nbtr .LE. 1 ) STOP '**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def'
     77  IF ( nbtr .LE. 1 ) THEN
     78    abort_message='**PHYTRAC**:initrrnpb:** nbtr <= 1; verifier PB dans traceur.def'
     79    CALL abort_gcm (modname,abort_message,1)
     80  ENDIF
    7181  it = 2
    7282  s = 0.                ! Pas de source

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/mod_phys_lmdz_omp_data.F90

    r1001 r1299  
    11!
    2 !$Header$
     2!$Id$
    33!
    44MODULE mod_phys_lmdz_omp_data
    … …  
    2727    INTEGER :: i
    2828
     29    CHARACTER (LEN=20) :: modname='Init_phys_lmdz_omp_data'
     30    CHARACTER (LEN=80) :: abort_message
     31
     32
    2933#ifdef CPP_OMP   
    3034    INTEGER :: OMP_GET_NUM_THREADS
    … …  
    5155     is_omp_root=.TRUE.
    5256   ELSE
    53      PRINT *,'ANORMAL : OMP_MASTER /= 0'
    54      STOP
     57     abort_message = 'ANORMAL : OMP_MASTER /= 0'
     58     CALL abort_gcm (modname,abort_message,1)
    5559   ENDIF
    5660!$OMP END MASTER

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/o3cm.F

    r524 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE o3cm (amb, bmb, sortie, ntab)
    … …  
    1919c======================================================================
    2020      external mbtozm
     21      CHARACTER (LEN=20) :: modname=''
     22      CHARACTER (LEN=80) :: abort_message
    2123c======================================================================
    2224c la fonction en ligne w(x) donne le profil de l'ozone en fonction
    … …  
    2729      w(x) = wp/h * EXP((x-xp)/h)/ (con+EXP((x-xp)/h))**2
    2830c======================================================================
    29       IF (ntab .GT. 499) STOP 'BIG ntab'
    30       xincr = (bmb-amb) / FLOAT(ntab)
     31      IF (ntab .GT. 499) THEN
     32        abort_message = 'BIG ntab'
     33        CALL abort_gcm (modname,abort_message,1)
     34      ENDIF
     35      xincr = (bmb-amb) / REAL(ntab)
    3136      xtab(1) = amb
    3237      DO n = 2, ntab

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/orografi.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE drag_noro (nlon,nlev,dtime,paprs,pplay,
    … …  
    14971497     *       ZDVDT(KLON)
    14981498      REAL ZHCRIT(KLON,KLEV)
     1499      CHARACTER (LEN=20) :: modname='orografi'
     1500      CHARACTER (LEN=80) :: abort_message
    14991501C-----------------------------------------------------------------------
    15001502C
    … …  
    15041506      LIFTHIGH=.FALSE.
    15051507
    1506       IF(NLON.NE.KLON.OR.NLEV.NE.KLEV)STOP
     1508      IF(NLON.NE.KLON.OR.NLEV.NE.KLEV)THEN
     1509        abort_message = 'pb dimension'
     1510        CALL abort_gcm (modname,abort_message,1)
     1511      ENDIF
    15071512      ZCONS1=1./RD
    15081513cym      KLEVM1=KLEV-1

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/orografi_strato.F

    r1001 r1299  
    8989      REAL pt(klon,klev), pu(klon,klev), pv(klon,klev)
    9090      REAL papmf(klon,klev),papmh(klon,klev+1)
     91      CHARACTER (LEN=20) :: modname='orografi_strato'
     92      CHARACTER (LEN=80) :: abort_message
    9193c
    9294c INITIALIZE OUTPUT VARIABLES
    … …  
    16801682      logical lifthigh
    16811683      real zcons1,ztmst
     1684      CHARACTER (LEN=20) :: modname='orolift_strato'
     1685      CHARACTER (LEN=80) :: abort_message
     1686
    16821687
    16831688C-----------------------------------------------------------------------
    … …  
    16881693      lifthigh=.false.
    16891694
    1690       if(nlon.ne.klon.or.nlev.ne.klev)stop
     1695      if(nlon.ne.klon.or.nlev.ne.klev) then
     1696        abort_message = 'pb dimension'
     1697        CALL abort_gcm (modname,abort_message,1)
     1698      ENDIF
    16911699      zcons1=1./rd
    16921700      ztmst=ptsphy

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/pbl_surface_mod.F90

    r1282 r1299  
    2222  USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
    2323  USE coef_diff_turb_mod,  ONLY : coef_diff_turb
     24  USE control_mod
     25
    2426
    2527  IMPLICIT NONE
    … …  
    257259    INCLUDE "YOETHF.h"
    258260    INCLUDE "temps.h"
    259     INCLUDE "control.h"
    260261! Input variables
    261262!****************************************************************************************
    … …  
    657658          tabindx(:)=0.
    658659          DO i=1,knon
    659              tabindx(i)=FLOAT(i)
     660             tabindx(i)=REAL(i)
    660661          END DO
    661662          debugtab(:,:) = 0.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phyredem.F

    r1298 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    1414      USE traclmdz_mod, ONLY : traclmdz_to_restart
    1515      USE infotrac
     16      USE control_mod
     17
    1618
    1719      IMPLICIT none
    … …  
    2426#include "dimsoil.h"
    2527#include "clesphys.h"
    26 #include "control.h"
    2728#include "temps.h"
    2829#include "thermcell.h"

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phys_state_var_mod.F90

    r1279 r1299  
    289289SUBROUTINE phys_state_var_init(read_climoz)
    290290use dimphy
     291USE control_mod
    291292use aero_mod
    292293IMPLICIT NONE
    … …  
    301302
    302303#include "indicesol.h"
    303 #include "control.h"
    304304      ALLOCATE(rlat(klon), rlon(klon))
    305305      ALLOCATE(pctsrf(klon,nbsrf))
    … …  
    411411SUBROUTINE phys_state_var_end
    412412use dimphy
     413use control_mod
    413414IMPLICIT NONE
    414415#include "indicesol.h"
    415 #include "control.h"
    416416
    417417      deallocate(rlat, rlon, pctsrf, ftsol, falb1, falb2)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/physiq.F

    r1298 r1299  
    3939      use conf_phys_m, only: conf_phys
    4040      use radlwsw_m, only: radlwsw
     41      USE control_mod
     42
    4143
    4244      IMPLICIT none
    … …  
    9799#include "dimsoil.h"
    98100#include "clesphys.h"
    99 #include "control.h"
    100101#include "temps.h"
    101102#include "iniprint.h"
    … …  
    13481349         ENDIF
    13491350c
    1350          IF (dtime*FLOAT(radpas).GT.21600..AND.cycle_diurne) THEN
     1351         IF (dtime*REAL(radpas).GT.21600..AND.cycle_diurne) THEN
    13511352           WRITE(lunout,*)'Nbre d appels au rayonnement insuffisant'
    13521353           WRITE(lunout,*)"Au minimum 4 appels par jour si cycle diurne"
    … …  
    15181519         CALL VTb(VTinca)
    15191520!         iii = MOD(NINT(xjour),360)
    1520 !         calday = FLOAT(iii) + jH_cur
    1521          calday = FLOAT(days_elapsed) + jH_cur
     1521!         calday = REAL(iii) + jH_cur
     1522         calday = REAL(days_elapsed) + jH_cur
    15221523         WRITE(lunout,*) 'initial time chemini', days_elapsed, calday
    15231524
    … …  
    18011802!   solarlong0
    18021803        if (solarlong0<-999.) then
    1803            CALL orbite(FLOAT(days_elapsed+1),zlongi,dist)
     1804           CALL orbite(REAL(days_elapsed+1),zlongi,dist)
    18041805        else
    18051806           zlongi=solarlong0  ! longitude solaire vraie
    … …  
    18121813!  Avec ou sans cycle diurne
    18131814      IF (cycle_diurne) THEN
    1814         zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
     1815        zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
    18151816        CALL zenang(zlongi,jH_cur,zdtime,rlat,rlon,rmu0,fract)
    18161817      ELSE
    … …  
    19611962
    19621963      IF (iflag_con.EQ.1) THEN
    1963           stop'reactiver le call conlmd dans physiq.F'
     1964        abort_message ='reactiver le call conlmd dans physiq.F'
     1965        CALL abort_gcm (modname,abort_message,1)
    19641966c     CALL conlmd (dtime, paprs, pplay, t_seri, q_seri, conv_q,
    19651967c    .             d_t_con, d_q_con,
    … …  
    22052207          za = 0.0
    22062208          DO i = 1, klon
    2207             za = za + airephy(i)/FLOAT(klon)
     2209            za = za + airephy(i)/REAL(klon)
    22082210            zx_t = zx_t + (rain_con(i)+
    2209      .                   snow_con(i))*airephy(i)/FLOAT(klon)
     2211     .                   snow_con(i))*airephy(i)/REAL(klon)
    22102212          ENDDO
    22112213          zx_t = zx_t/za*dtime
    … …  
    25992601         za = 0.0
    26002602         DO i = 1, klon
    2601             za = za + airephy(i)/FLOAT(klon)
     2603            za = za + airephy(i)/REAL(klon)
    26022604            zx_t = zx_t + (rain_lsc(i)
    2603      .                  + snow_lsc(i))*airephy(i)/FLOAT(klon)
     2605     .                  + snow_lsc(i))*airephy(i)/REAL(klon)
    26042606        ENDDO
    26052607         zx_t = zx_t/za*dtime
    … …  
    28232825         CALL VTe(VTphysiq)
    28242826         CALL VTb(VTinca)
    2825          calday = FLOAT(days_elapsed + 1) + jH_cur
     2827         calday = REAL(days_elapsed + 1) + jH_cur
    28262828
    28272829         call chemtime(itap+itau_phy-1, date0, dtime)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phystokenc.F

    r1146 r1299  
    1313      USE infotrac, ONLY : nqtot
    1414      USE iophy
     15      USE control_mod
     16
    1517      IMPLICIT none
    1618
    … …  
    2426#include "tracstoke.h"
    2527#include "indicesol.h"
    26 #include "control.h"
    2728c======================================================================
    2829

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phytrac.F90

    r1279 r1299  
    3333  USE traclmdz_mod
    3434  USE tracinca_mod
     35  USE control_mod
     36
    3537
    3638
    … …  
    4345  INCLUDE "temps.h"
    4446  INCLUDE "paramet.h"
    45   INCLUDE "control.h"
    4647  INCLUDE "thermcell.h"
    4748!==========================================================================

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/read_pstoke.F

    r1146 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    1818C******************************************************************************
    1919
    20         use netcdf
    21        USE dimphy
     20      use netcdf
     21      USE dimphy
     22      USE control_mod
     23
    2224       IMPLICIT NONE
    2325
    … …  
    3335#include "serre.h"
    3436#include "indicesol.h"
    35 #include "control.h"
    3637cccc#include "dimphy.h"
    3738       

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/read_pstoke0.F

    r1146 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44c
    … …  
    1919        use netcdf
    2020       USE dimphy
     21       USE control_mod
     22
    2123       IMPLICIT NONE
    2224
    … …  
    3234#include "serre.h"
    3335#include "indicesol.h"
    34 #include "control.h"
    3536cccc#include "dimphy.h"
    3637         

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/readaerosol.F90

    r1279 r1299  
    128128                 DO i = 1, klon
    129129                    pt_out(i,k,it) = &
    130                          pt_out(i,k,it) - FLOAT(iyr_in-iyr1)/FLOAT(iyr2-iyr1) * &
     130                         pt_out(i,k,it) - REAL(iyr_in-iyr1)/REAL(iyr2-iyr1) * &
    131131                         (pt_out(i,k,it) - pt_2(i,k,it))
    132132                 END DO
    … …  
    135135              DO i = 1, klon
    136136                 psurf(i,it) = &
    137                       psurf(i,it) - FLOAT(iyr_in-iyr1)/FLOAT(iyr2-iyr1) * &
     137                      psurf(i,it) - REAL(iyr_in-iyr1)/REAL(iyr2-iyr1) * &
    138138                      (psurf(i,it) - psurf2(i,it))
    139139
    140140                 load(i,it) = &
    141                       load(i,it) - FLOAT(iyr_in-iyr1)/FLOAT(iyr2-iyr1) * &
     141                      load(i,it) - REAL(iyr_in-iyr1)/REAL(iyr2-iyr1) * &
    142142                      (load(i,it) - load2(i,it))
    143143              END DO
    … …  
    493493                spole = spole + varyear(i,jjm+1,k,imth)
    494494             END DO
    495              npole = npole/FLOAT(iim)
    496              spole = spole/FLOAT(iim)
     495             npole = npole/REAL(iim)
     496             spole = spole/REAL(iim)
    497497             varyear(:,1,    k,imth) = npole
    498498             varyear(:,jjm+1,k,imth) = spole

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/readaerosol_interp.F90

    r1279 r1299  
    126126  IF(mpi_rank == 0 .AND. debug)then
    127127     ! 0.02 is about 0.5/24, namly less than half an hour
    128      OLDNEWDAY = (r_day-FLOAT(iday) < 0.02)
     128     OLDNEWDAY = (r_day-REAL(iday) < 0.02)
    129129     ! Once per day, update aerosol fields
    130130     lmt_pas = NINT(86400./pdtphys)
    131      PRINT*,'r_day-FLOAT(iday) =',r_day-FLOAT(iday)
     131     PRINT*,'r_day-REAL(iday) =',r_day-REAL(iday)
    132132     PRINT*,'itap =',itap
    133133     PRINT*,'pdtphys =',pdtphys
    … …  
    233233!
    234234     DO i = 2, 13
    235        month_len(i) = float(ioget_mon_len(year_cur, i-1))
     235       month_len(i) = REAL(ioget_mon_len(year_cur, i-1))
    236236       CALL ymds2ju(year_cur, i-1, 1, 0.0, month_start(i))
    237237     ENDDO
    238      month_len(1) = float(ioget_mon_len(year_cur-1, 12))
     238     month_len(1) = REAL(ioget_mon_len(year_cur-1, 12))
    239239     CALL ymds2ju(year_cur-1, 12, 1, 0.0, month_start(1))
    240      month_len(14) = float(ioget_mon_len(year_cur+1, 1))
     240     month_len(14) = REAL(ioget_mon_len(year_cur+1, 1))
    241241     CALL ymds2ju(year_cur+1, 1, 1, 0.0, month_start(14))
    242242     month_mid(:) = month_start (:) + month_len(:)/2.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/surf_ocean_mod.F90

    r1146 r1299  
    141141!****************************************************************************************
    142142    IF ( MINVAL(rmu0) == MAXVAL(rmu0) .AND. MINVAL(rmu0) == -999.999 ) THEN
    143        CALL alboc(FLOAT(jour),rlat,alb_eau)
     143       CALL alboc(REAL(jour),rlat,alb_eau)
    144144    ELSE  ! diurnal cycle
    145145       CALL alboc_cd(rmu0,alb_eau)

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/sw_aeroAR4.F90

    r1279 r1299  
    186186!$OMP THREADPRIVATE(AEROSOLFEEDBACK_ACTIVE) 
    187187
     188      CHARACTER (LEN=20) :: modname='sw_aeroAR4'
     189      CHARACTER (LEN=80) :: abort_message
     190
    188191  IF ((.not. ok_ade) .and. (AEROSOLFEEDBACK_ACTIVE .ge. 2)) THEN
    189      print*,'Error: direct effect is not activated but assumed to be active - see sw_aeroAR4.F90'
    190      stop
     192     abort_message ='Error: direct effect is not activated but assumed to be active - see sw_aeroAR4.F90'
     193     CALL abort_gcm (modname,abort_message,1)
    191194  ENDIF
    192195  AEROSOLFEEDBACK_ACTIVE=MIN(MAX(AEROSOLFEEDBACK_ACTIVE,0),3)
    193196  IF  (AEROSOLFEEDBACK_ACTIVE .gt. 3) THEN
    194      print*,'Error: AEROSOLFEEDBACK_ACTIVE options go only until 3'
    195      stop
     197     abort_message ='Error: AEROSOLFEEDBACK_ACTIVE options go only until 3'
     198     CALL abort_gcm (modname,abort_message,1)
    196199  ENDIF
    197200

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell.F

    r987 r1299  
     1!
     2! $Id$
     3!
    14      SUBROUTINE calcul_sec(ngrid,nlay,ptimestep
    25     s                  ,pplay,pplev,pphi,zlev
    … …  
    132135      character*10 str10
    133136
     137      character (len=20) :: modname='calcul_sec'
     138      character (len=80) :: abort_message
     139
     140
    134141!      LOGICAL vtest(klon),down
    135142
    … …  
    530537c      write(8,*)zmax(ig),lmax(ig),lentr(ig),lmin(ig)
    531538      enddo
    532 con stoppe après les calculs de zmax et wmax
     539con stope après les calculs de zmax et wmax
    533540      RETURN
    534541
    … …  
    776783         do ig=1,ngrid
    777784            if(fracd(ig,l).lt.0.1) then
    778                stop'fracd trop petit'
     785              abort_message = 'fracd trop petit'
     786              CALL abort_gcm (modname,abort_message,1)
     787
    779788            else
    780789c    vitesse descendante "diagnostique"
    … …  
    860869cRC
    861870      if (w2di.eq.1) then
    862          fm0=fm0+ptimestep*(fm-fm0)/float(tho)
    863          entr0=entr0+ptimestep*(entr-entr0)/float(tho)
     871         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
     872         entr0=entr0+ptimestep*(entr-entr0)/REAL(tho)
    864873      else
    865874         fm0=fm

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcellV0_main.F90

    r1294 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE thermcellV0_main(itap,ngrid,nlay,ptimestep  &
    … …  
    165165      character*2 str2
    166166      character*10 str10
     167
     168      character (len=20) :: modname='thermcellV0_main'
     169      character (len=80) :: abort_message
    167170
    168171      EXTERNAL SCOPY
    … …  
    484487! Test valable seulement en 1D mais pas genant
    485488      if (.not. (f0(1).ge.0.) ) then
    486            stop 'Dans thermcell_main'
     489        abort_message = 'Dans thermcell_main f0(1).lt.0 '
     490        CALL abort_gcm (modname,abort_message,1)
    487491      endif
    488492
    … …  
    827831                  write(6,'(i3,7f10.3)') k,pplay(i,k),ztv(i,k),1000*po(i,k),ztva(i,k),1000*zqla(i,k),f_star(i,k),zw2(i,k)
    828832               enddo
    829 !              stop
    830833           endif
    831834         enddo
    … …  
    863866
    864867      REAL f(ngrid)
     868
     869      character (len=20) :: modname='thermcellV0_main'
     870      character (len=80) :: abort_message
    865871
    866872      do ig=1,ngrid
    … …  
    885891                print*,'zmax_sec',zmax_sec(ig)
    886892                print*,'wmax_sec',wmax_sec(ig)
    887                 stop
     893                abort_message = 'zdenom<1.e-14'
     894                CALL abort_gcm (modname,abort_message,1)
    888895             endif
    889896             if ((zmax_sec(ig).gt.1.e-10).and.(iflag_thermals_ed.eq.0)) then
    … …  
    13381345!        if(detr_star(ig,l).GT.1.) THEN
    13391346!         print*,'th_plumeAFT ig l detr_star f0 1/f0',ig,l,detr_star(ig,l),f0(ig), &
    1340 !   &     float(1)/f0(ig)
     1347!   &     REAL(1)/f0(ig)
    13411348!        endif
    13421349!IM 060508 end
    … …  
    16561663
    16571664            if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then
    1658 !               stop'On tombe sur le cas particulier de thermcell_dry'
    16591665                print*,'On tombe sur le cas particulier de thermcell_plume'
    16601666                zw2(ig,l+1)=0.
    … …  
    18321838
    18331839            if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then
    1834 !               stop'On tombe sur le cas particulier de thermcell_dry'
    18351840!               print*,'On tombe sur le cas particulier de thermcell_dry'
    18361841                zw2(ig,l+1)=0.

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_dq.F90

    r983 r1299  
    3131      real ztimestep
    3232      integer niter,iter
     33      CHARACTER (LEN=20) :: modname='thermcell_dq'
     34      CHARACTER (LEN=80) :: abort_message
    3335
    3436
    … …  
    4244            if (entr(ig,k).gt.zzm) then
    4345               print*,'entr dt > m ',entr(ig,k)*ptimestep,masse(ig,k)
    44                stop
     46               abort_message = ''
     47               CALL abort_gcm (modname,abort_message,1)
    4548            endif
    4649         enddo

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_dry.F90

    r1294 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44       SUBROUTINE thermcell_dry(ngrid,nlay,zlev,pphi,ztv,alim_star,  &
    … …  
    3939       REAL linter(ngrid),zlevinter(ngrid)
    4040       INTEGER lmix(ngrid),lmax(ngrid),lmin(ngrid)
     41      CHARACTER (LEN=20) :: modname='thermcell_dry'
     42      CHARACTER (LEN=80) :: abort_message
    4143
    4244!initialisations

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_flux.F90

    r1146 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44
    … …  
    5151!$OMP THREADPRIVATE(fomass_max,alphamax)
    5252
     53      character (len=20) :: modname='thermcell_flux'
     54      character (len=80) :: abort_message
     55
    5356      fomass_max=0.5
    5457      alphamax=0.7
    … …  
    9295                    print*,'alim_star(ig,l)',alim_star(ig,l)
    9396                    print*,'detr_star(ig,l)',detr_star(ig,l)
    94 !                   stop
    9597               endif
    9698            else
    … …  
    100102                    print*,'alim_star(ig,l)',alim_star(ig,l)
    101103                    print*,'detr_star(ig,l)',detr_star(ig,l)
    102                     stop
     104                    abort_message = ''
     105                    CALL abort_gcm (modname,abort_message,1)
    103106               endif
    104107            endif
    … …  
    264267            if (entr(ig,l)<0.) then
    265268               print*,'N1 ig,l,entr',ig,l,entr(ig,l)
    266                stop 'entr negatif'
     269               abort_message = 'entr negatif'
     270               CALL abort_gcm (modname,abort_message,1)
    267271            endif
    268272            if (detr(ig,l).gt.fm(ig,l)) then
    … …  
    292296               print*,'entr(ig,l)',entr(ig,l)
    293297               print*,'fm(ig,l)',fm(ig,l)
    294                stop 'probleme dans thermcell flux'
     298               abort_message = 'probleme dans thermcell flux'
     299               CALL abort_gcm (modname,abort_message,1)
    295300            endif
    296301         enddo
    … …  
    319324               print*,'detr(ig,l)',detr(ig,l)
    320325               print*,'fm(ig,l)',fm(ig,l)
    321                stop 'probleme dans thermcell flux'
     326               abort_message = 'probleme dans thermcell flux'
     327               CALL abort_gcm (modname,abort_message,1)
    322328            endif
    323329        enddo
    … …  
    420426                         print*,'fm(ig,l+1)',fm(ig,l+1)
    421427                         print*,'fm(ig,l)',fm(ig,l)
    422                          stop 'probleme dans thermcell_flux'
     428                         abort_message = 'probleme dans thermcell_flux'
     429                         CALL abort_gcm (modname,abort_message,1)
    423430                      endif
    424431                      entr(ig,l+1)=entr(ig,l+1)-ddd
    … …  
    478485      character*3 descr
    479486
     487      character (len=20) :: modname='thermcell_flux'
     488      character (len=80) :: abort_message
     489
    480490      lm=lmax(igout)+5
    481491      if(lm.gt.klev) lm=klev
    … …  
    500510          print*,'detr(igout,l)',detr(igout,l)
    501511          print*,'fm(igout,l)',fm(igout,l)
    502           stop
     512          abort_message = ''
     513          CALL abort_gcm (modname,abort_message,1)
    503514          endif
    504515      enddo

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_flux2.F90

    r1146 r1299  
     1!
     2! $Id$
     3!
    14      SUBROUTINE thermcell_flux2(ngrid,klev,ptimestep,masse, &
    25     &       lalim,lmax,alim_star,  &
    … …  
    4750      save fomass_max,alphamax
    4851
     52      character (len=20) :: modname='thermcell_flux2'
     53      character (len=80) :: abort_message
     54
    4955      fomass_max=0.5
    5056      alphamax=0.7
    … …  
    8894                    print*,'alim_star(ig,l)',alim_star(ig,l)
    8995                    print*,'detr_star(ig,l)',detr_star(ig,l)
    90 !                   stop
    9196               endif
    9297            else
    … …  
    96101                    print*,'alim_star(ig,l)',alim_star(ig,l)
    97102                    print*,'detr_star(ig,l)',detr_star(ig,l)
    98                     stop
     103                    abort_message = ''
     104                    CALL abort_gcm (modname,abort_message,1)
    99105               endif
    100106            endif
    … …  
    256262            if (entr(ig,l)<0.) then
    257263               print*,'N1 ig,l,entr',ig,l,entr(ig,l)
    258                stop 'entr negatif'
     264               abort_message = 'entr negatif'
     265               CALL abort_gcm (modname,abort_message,1)
    259266            endif
    260267            if (detr(ig,l).gt.fm(ig,l)) then
    … …  
    285292               print*,'entr(ig,l)',entr(ig,l)
    286293               print*,'fm(ig,l)',fm(ig,l)
    287                stop 'probleme dans thermcell flux'
     294               abort_message = 'probleme dans thermcell flux'
     295               CALL abort_gcm (modname,abort_message,1)
    288296            endif
    289297         enddo
    … …  
    312320               print*,'detr(ig,l)',detr(ig,l)
    313321               print*,'fm(ig,l)',fm(ig,l)
    314                stop 'probleme dans thermcell flux'
     322               abort_message = 'probleme dans thermcell flux'
     323               CALL abort_gcm (modname,abort_message,1)
    315324            endif
    316325        enddo
    … …  
    413422                         print*,'fm(ig,l+1)',fm(ig,l+1)
    414423                         print*,'fm(ig,l)',fm(ig,l)
    415                          stop 'probleme dans thermcell_flux'
     424                         abort_message = 'probleme dans thermcell_flux'
     425                         CALL abort_gcm (modname,abort_message,1)
    416426                      endif
    417427                      entr(ig,l+1)=entr(ig,l+1)-ddd

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_main.F90

    r1294 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE thermcell_main(itap,ngrid,nlay,ptimestep  &
    … …  
    170170      character*2 str2
    171171      character*10 str10
     172
     173      character (len=20) :: modname='thermcell_main'
     174      character (len=80) :: abort_message
    172175
    173176      EXTERNAL SCOPY
    … …  
    466469! Test valable seulement en 1D mais pas genant
    467470      if (.not. (f0(1).ge.0.) ) then
    468            stop'Dans thermcell_main'
     471              abort_message = '.not. (f0(1).ge.0.)'
     472              CALL abort_gcm (modname,abort_message,1)
    469473      endif
    470474
    … …  
    781785      if (prt_level.ge.1) print*,'thermcell_main FIN  OK'
    782786
    783 !     if(icount.eq.501) stop'au pas 301 dans thermcell_main'
    784787      return
    785788      end
    … …  
    817820                  write(6,'(i3,7f10.3)') k,pplay(i,k),ztv(i,k),1000*po(i,k),ztva(i,k),1000*zqla(i,k),f_star(i,k),zw2(i,k)
    818821               enddo
    819 !              stop
    820822           endif
    821823         enddo

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_old.F

    r987 r1299  
    112112      character (len=10) :: str10
    113113
     114      character (len=20) :: modname='thermcell2002'
     115      character (len=80) :: abort_message
     116
    114117      LOGICAL vtest(klon),down
    115118
    … …  
    336339            if(w2di.eq.2) then
    337340               entr(ig,k)=entr(ig,k)+
    338      s         ptimestep*(zzz-entr(ig,k))/float(tho)
     341     s         ptimestep*(zzz-entr(ig,k))/REAL(tho)
    339342            else
    340343               entr(ig,k)=zzz
    … …  
    379382c     print*,'ig,l+1,ztv(ig,l+1)'
    380383c     print*, ig,l+1,ztv(ig,l+1)
    381 c        stop'dans thermiques'
    382384c     endif
    383385               zw2(ig,l+1)=2.*RG*(ztv(ig,l)-ztv(ig,l+1))/ztv(ig,l+1)
    … …  
    395397c     print*,'Tv ',(ztv(ig,ll),ll=1,klev)
    396398c     print*,'Entr ',(entr(ig,ll),ll=1,klev)
    397 c        stop'dans thermiques'
    398399c     endif
    399400               ztva(ig,l)=(fmc(ig,l)*ztva(ig,l-1)+entr(ig,l)*ztv(ig,l))
    … …  
    517518         do ig=1,ngrid
    518519            if(fracd(ig,l).lt.0.1) then
    519                stop'fracd trop petit'
    520             else
     520              abort_message = 'fracd trop petit'
     521              CALL abort_gcm (modname,abort_message,1)
     522           else
    521523c    vitesse descendante "diagnostique"
    522524               wd(ig,l)=fm(ig,l)/(fracd(ig,l)*rhobarz(ig,l))
    … …  
    588590
    589591      if (w2di.eq.1) then
    590          fm0=fm0+ptimestep*(fm-fm0)/float(tho)
    591          entr0=entr0+ptimestep*(entr-entr0)/float(tho)
     592         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
     593         entr0=entr0+ptimestep*(entr-entr0)/REAL(tho)
    592594      else
    593595         fm0=fm
    … …  
    10001002      character*2 str2
    10011003      character*10 str10
     1004
     1005      character (len=20) :: modname='thermcell_cld'
     1006      character (len=80) :: abort_message
    10021007
    10031008      LOGICAL vtest(klon),down
    … …  
    18551860       if (l.eq.klev) then
    18561861          print*,'THERMCELL PB ig=',ig,'   l=',l
    1857           stop
     1862          abort_message = 'THERMCELL PB'
     1863          CALL abort_gcm (modname,abort_message,1)
    18581864       endif
    18591865!       if ((zw2(ig,l+1).gt.1.e-10).and.(zw2(ig,l).gt.1.e-10).and.
    … …  
    21642170         do ig=1,ngrid
    21652171            if(fracd(ig,l).lt.0.1.and.(test(ig).eq.1)) then
    2166                stop'fracd trop petit'
     2172              abort_message = 'fracd trop petit'
     2173              CALL abort_gcm (modname,abort_message,1)
    21672174            else
    21682175c    vitesse descendante "diagnostique"
    … …  
    22622269
    22632270      if (w2di.eq.1) then
    2264          fm0=fm0+ptimestep*(fm-fm0)/float(tho)
    2265          entr0=entr0+ptimestep*(alim+entr-entr0)/float(tho)
     2271         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
     2272         entr0=entr0+ptimestep*(alim+entr-entr0)/REAL(tho)
    22662273      else
    22672274         fm0=fm
    … …  
    27472754      character*10 str10
    27482755
     2756      character (len=20) :: modname='thermcell_eau'
     2757      character (len=80) :: abort_message
     2758
    27492759      LOGICAL vtest(klon),down
    27502760      LOGICAL Zsat(klon)
    … …  
    34103420         do ig=1,ngrid
    34113421            if(fracd(ig,l).lt.0.1) then
    3412                stop'fracd trop petit'
     3422              abort_message = 'fracd trop petit'
     3423              CALL abort_gcm (modname,abort_message,1)
    34133424            else
    34143425c    vitesse descendante "diagnostique"
    … …  
    34813492
    34823493      if (w2di.eq.1) then
    3483          fm0=fm0+ptimestep*(fm-fm0)/float(tho)
    3484          entr0=entr0+ptimestep*(entr-entr0)/float(tho)
     3494         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
     3495         entr0=entr0+ptimestep*(entr-entr0)/REAL(tho)
    34853496      else
    34863497         fm0=fm
    … …  
    38483859      character*10 str10
    38493860
     3861      character (len=20) :: modname='thermcell'
     3862      character (len=80) :: abort_message
     3863
    38503864      LOGICAL vtest(klon),down
    38513865
    … …  
    43944408         do ig=1,ngrid
    43954409            if(fracd(ig,l).lt.0.1) then
    4396                stop'fracd trop petit'
     4410              abort_message = 'fracd trop petit'
     4411              CALL abort_gcm (modname,abort_message,1)
    43974412            else
    43984413c    vitesse descendante "diagnostique"
    … …  
    44774492cRC
    44784493      if (w2di.eq.1) then
    4479          fm0=fm0+ptimestep*(fm-fm0)/float(tho)
    4480          entr0=entr0+ptimestep*(entr-entr0)/float(tho)
     4494         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
     4495         entr0=entr0+ptimestep*(entr-entr0)/REAL(tho)
    44814496      else
    44824497         fm0=fm
    … …  
    52575272      character*10 str10
    52585273
     5274      character (len=20) :: modname='thermcell_sec'
     5275      character (len=80) :: abort_message
     5276
    52595277      LOGICAL vtest(klon),down
    52605278
    … …  
    58225840         do ig=1,ngrid
    58235841            if(fracd(ig,l).lt.0.1) then
    5824                stop'fracd trop petit'
     5842              abort_message = 'fracd trop petit'
     5843              CALL abort_gcm (modname,abort_message,1)
    58255844            else
    58265845c    vitesse descendante "diagnostique"
    … …  
    59055924cRC
    59065925      if (w2di.eq.1) then
    5907          fm0=fm0+ptimestep*(fm-fm0)/float(tho)
    5908          entr0=entr0+ptimestep*(entr-entr0)/float(tho)
     5926         fm0=fm0+ptimestep*(fm-fm0)/REAL(tho)
     5927         entr0=entr0+ptimestep*(entr-entr0)/REAL(tho)
    59095928      else
    59105929         fm0=fm

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_plume.F90

    r1294 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE thermcell_plume(itap,ngrid,klev,ptimestep,ztv,zthl,po,zl,rhobarz,  &
    … …  
    436436            if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then
    437437!               stop'On tombe sur le cas particulier de thermcell_dry'
    438                 print*,'On tombe sur le cas particulier de thermcell_plume'
     438                write(lunout,*)                                         &
     439     &          'On tombe sur le cas particulier de thermcell_plume'
    439440                zw2(ig,l+1)=0.
    440441                linter(ig)=l+1

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/tracinca_mod.F90

    r1279 r1299  
    4545    USE vampir
    4646    USE comgeomphy
     47    USE control_mod
     48
    4749   
    4850    IMPLICIT NONE
    4951   
    5052    INCLUDE "indicesol.h"
    51     INCLUDE "control.h"
    5253    INCLUDE "dimensions.h"
    5354    INCLUDE "paramet.h"
    … …  
    125126    CALL VTb(VTinca)
    126127   
    127     calday = FLOAT(julien) + gmtime
     128    calday = REAL(julien) + gmtime
    128129    ncsec  = NINT (86400.*gmtime)
    129130     

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/undefSTD.F

    r1279 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE undefSTD(nlevSTD,itap,tlevSTD,
    … …  
    1313c I. Musat : 09.2004
    1414c
    15 c Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef))
     15c Calcul * du nombre de pas de temps (REAL(ecrit_XXX)-tnondef))
    1616c          ou la variable tlevSTD est bien definie (.NE.1.E+20),
    1717c et

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/write_histISCCP.h

    r1045 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      IF (ok_isccp) THEN
    … …  
    7777     .                 meantaucld(:,n))
    7878c
    79         zx_tmp_fi2d(1:klon)=float(seed(1:klon,n))
     79        zx_tmp_fi2d(1:klon)=REAL(seed(1:klon,n))
    8080c
    8181c       print*,'n=',n,' write_ISCCP avant seed'

  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/yamada4.F

    r938 r1299  
    11!
    2 ! $Header$
     2! $Id$
    33!
    44      SUBROUTINE yamada4(ngrid,dt,g,rconst,plev,temp
    … …  
    106106c$OMP THREADPRIVATE(rino,smyam,styam,lyam,knyam,w2yam,t2yam)
    107107      logical,save :: firstcall=.true.
     108
     109      character (len=20) :: modname='yamada4'
     110      character (len=80) :: abort_message
     111
    108112c$OMP THREADPRIVATE(firstcall)       
    109113      frif(ri)=0.6588*(ri+0.1776-sqrt(ri*ri-0.3221*ri+0.03156))
    … …  
    128132
    129133      if (.not.(iflag_pbl.ge.6.and.iflag_pbl.le.9)) then
    130            stop'probleme de coherence dans appel a MY'
     134              abort_message = 'probleme de coherence dans appel a MY'
     135              CALL abort_gcm (modname,abort_message,1)
    131136      endif
    132137
Note: See TracChangeset for help on using the changeset viewer.