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/dyn3d
Files:
1 added
43 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)
Note: See TracChangeset for help on using the changeset viewer.