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/dyn3dpar
Files:
1 added
41 edited

Legend:

Unmodified
Added
Removed

  • 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
Note: See TracChangeset for help on using the changeset viewer.