Changeset 1960


Ignore:
Timestamp:
Feb 7, 2014, 5:20:25 PM (10 years ago)
Author:
fhourdin
Message:

Nettoyage du code 1D pour limiter les warning en mode debug.
Essentiellement des declarations inutiles.

Cleaning of the 1D code in order to limit warnings in debug mode.
Essentially supressing unused declared variables.

Location:
LMDZ5/trunk/libf/phy1d
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phy1d/1DUTILS.h

    r1953 r1960  
    44c
    55c
    6       SUBROUTINE conf_unicol( tapedef )
     6      SUBROUTINE conf_unicol
    77c
    88#ifdef CPP_IOIPSL
    … …  
    1515c-----------------------------------------------------------------------
    1616c     Auteurs :   A. Lahellec  .
    17 c
    18 c     Arguments :
    19 c
    20 c     tapedef   :
    21 
    22        INTEGER tapedef
    2317c
    2418c   Declarations :
    … …  
    367361c   Variables locales pour NetCDF:
    368362c   ------------------------------
    369       INTEGER nid, nvarid
    370       INTEGER idim_s
    371       INTEGER ierr, ierr_file
    372363      INTEGER iq
    373364      INTEGER length
    … …  
    378369      character*80 abort_message
    379370      LOGICAL found
    380 c
    381       INTEGER nb
    382371
    383372      modname = 'dyn1deta0 : '
    … …  
    508497c   ----------
    509498      CHARACTER*(*) fichnom
    510       REAL time
    511499cAl1 plev tronque pour .nc mais plev(klev+1):=0
    512500      real :: plev(klon,klev),play (klon,klev),phi(klon,klev)
    … …  
    520508c   Variables locales pour NetCDF:
    521509c   ------------------------------
    522       INTEGER nid, nvarid
    523       INTEGER idim_s
    524       INTEGER ierr, ierr_file
     510      INTEGER nid
     511      INTEGER ierr
    525512      INTEGER iq,l
    526513      INTEGER length
    … …  
    535522      DATA nb / 0 /
    536523
    537       REAL zan0,zjulian,hours
    538       INTEGER yyears0,jjour0, mmois0
    539       character*30 unites
    540 
    541 cDbg
    542524      CALL open_restartphy(fichnom)
    543525      print*,'redm1 ',fichnom,klon,klev,nqtot
    … …  
    550532      ierr = NF_OPEN(fichnom, NF_WRITE, nid)
    551533      IF (ierr .NE. NF_NOERR) THEN
    552          PRINT*, "Pb. d ouverture "//fichnom
    553          CALL abort
     534         abort_message="Pb. d ouverture "//fichnom
     535         CALL abort_gcm('Modele 1D',abort_message,1)
    554536      ENDIF
    555537
    … …  
    661643!   traitement des point normaux
    662644         DO j=2,jm-1
    663             ig=2+(j-2)*(im-1)
     645            ig=2+(j-2)*(im-1)
    664646            CALL SCOPY(im-1,pfi(ig,ifield),1,pdyn(1,j,ifield),1)
    665             pdyn(im,j,ifield)=pdyn(1,j,ifield)
     647            pdyn(im,j,ifield)=pdyn(1,j,ifield)
    666648         ENDDO
    667649      ENDDO
    … …  
    992974      DO ifield=1,nfield
    993975         DO j=2,jm-1
    994             ig=2+(j-2)*(im-1)
     976            ig=2+(j-2)*(im-1)
    995977            CALL SCOPY(im-1,pdyn(1,j,ifield),1,pfi(ig,ifield),1)
    996978         ENDDO
    … …  
    11511133
    11521134!======================================================================
    1153        SUBROUTINE read_tsurf1d(knon,knindex,sst_out)
     1135       SUBROUTINE read_tsurf1d(knon,sst_out)
    11541136
    11551137! This subroutine specifies the surface temperature to be used in 1D simulations
    … …  
    11581140
    11591141      INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
    1160       INTEGER, DIMENSION(klon), INTENT(IN) :: knindex  ! grid point number for compressed grid
    11611142      REAL, DIMENSION(klon), INTENT(OUT)   :: sst_out  ! tsurf used to force the single-column model
    11621143
    … …  
    12201201
    12211202       SUBROUTINE advect_va(llm,omega,d_t_va,d_q_va,d_u_va,d_v_va,
    1222      !                q,temp,u,v,
    1223      !            play,plev)
     1203     s                q,temp,u,v,play)
    12241204!itlmd
    12251205!----------------------------------------------------------------------
    … …  
    12371217        real  q(llm,3),temp(llm)
    12381218        real  u(llm),v(llm)
    1239         real  play(llm),plev(llm+1)
     1219        real  play(llm)
    12401220! interne
    12411221        integer l
    … …  
    13231303        real dph(llm),dqdp(llm),dtdp(llm)
    13241304! interne
    1325         integer l,k
    1326         real alpha,omdn,omup
     1305        integer k
     1306        real omdn,omup
    13271307
    13281308!        dudp=0.
    … …  
    14031383      character*80 fich_toga
    14041384
    1405       integer no,l,k,ip
    1406       real riy,rim,rid,rih,bid
     1385      integer k,ip
     1386      real bid
    14071387
    14081388      integer iy,im,id,ih
    … …  
    14221402
    14231403       do k = 1, nlev_toga
    1424          read(21,230) plev_toga(k,ip), t_toga(k,ip), q_toga(k,ip) 
     1404         read(21,230) plev_toga(k,ip), t_toga(k,ip), q_toga(k,ip)
    14251405     :       ,u_toga(k,ip), v_toga(k,ip), w_toga(k,ip)
    14261406     :       ,ht_toga(k,ip), vt_toga(k,ip), hq_toga(k,ip), vq_toga(k,ip)
    … …  
    14431423
    14441424  223 format(4i3,6f8.2)
    1445   226 format(f7.1,1x,10f8.2)
    1446   227 format(f7.1,1x,1p,4e11.3)
    14471425  230 format(6f9.3,4e11.3)
    14481426
    … …  
    14621440      character*80 fich_sandu
    14631441
    1464       integer no,l,k,ip
    1465       real riy,rim,rid,rih,bid
    1466 
     1442      integer ip
    14671443      integer iy,im,id,ih
    14681444
    1469        real plev_min
    1470 
    1471        plev_min = 55000.  ! pas de tendance de vap. d eau au-dessus de 55 hPa
     1445      real plev_min
     1446
     1447      print*,'nlev_sandu',nlev_sandu
     1448      plev_min = 55000.  ! pas de tendance de vap. d eau au-dessus de 55 hPa
    14721449
    14731450      open(21,file=trim(fich_sandu),form='formatted')
    … …  
    14821459
    14831460  223 format(4i3,f8.2)
    1484   226 format(f7.1,1x,10f8.2)
    1485   227 format(f7.1,1x,1p,4e11.3)
    1486   230 format(6f9.3,4e11.3)
    14871461
    14881462          return
    … …  
    15041478      character*80 fich_astex
    15051479
    1506       integer no,l,k,ip
    1507       real riy,rim,rid,rih,bid
    1508 
     1480      integer ip
    15091481      integer iy,im,id,ih
    15101482
    15111483       real plev_min
    15121484
     1485      print*,'nlev_astex',nlev_astex
    15131486       plev_min = 55000.  ! pas de tendance de vap. d eau au-dessus de 55 hPa
    15141487
    … …  
    15281501
    15291502  223 format(4i3,e13.2,f7.2,f7.3,f7.2,f7.3,f7.2)
    1530   226 format(f7.1,1x,10f8.2)
    1531   227 format(f7.1,1x,1p,4e11.3)
    1532   230 format(6f9.3,4e11.3)
    15331503
    15341504          return
    … …  
    15511521      character*80 :: fich_twpice
    15521522      real*8 time(ntime)
    1553       real*8 lat, lon, alt, phis       
     1523      real*8 lat, lon, alt, phis
    15541524      real*8 lev(nlevel)
    15551525      real*8 plev(nlevel,ntime)
    … …  
    15611531      real*8 T_adv_h(nlevel,ntime)
    15621532      real*8 T_adv_v(nlevel,ntime), q_adv_h(nlevel,ntime)
    1563       real*8 q_adv_v(nlevel,ntime)     
     1533      real*8 q_adv_v(nlevel,ntime)
    15641534      real*8 s(nlevel,ntime), s_adv_h(nlevel,ntime)
    15651535      real*8 s_adv_v(nlevel,ntime)
    … …  
    19761946         integer ierr
    19771947
    1978          integer i
    19791948         integer timevar,levvar
    19801949         integer timelen,levlen
    … …  
    20502019       real omega_mod(llm),o3mmr_mod(llm)
    20512020
    2052        integer l,k,k1,k2,kp
    2053        real aa,frac,frac1,frac2,fact
     2021       integer l,k,k1,k2
     2022       real frac,frac1,frac2,fact
    20542023
    20552024       do l = 1, llm
    … …  
    21682137       real o3mmr_mod(llm),ql_mod(llm),qt_mod(llm)
    21692138
    2170        integer l,k,k1,k2,kp
    2171        real aa,frac,frac1,frac2,fact
     2139       integer l,k,k1,k2
     2140       real frac,frac1,frac2,fact
    21722141
    21732142       do l = 1, llm
    … …  
    24442413        real ts_prof
    24452414! local:
    2446         integer it_sandu1, it_sandu2,k
     2415        integer it_sandu1, it_sandu2
    24472416        real timeit,time_sandu1,time_sandu2,frac
    24482417! Check that initial day of the simulation consistent with SANDU period:
    … …  
    25112480      character*80 fich_armcu
    25122481
    2513       integer no,l,k,ip
    2514       real riy,rim,rid,rih,bid
     2482      integer ip
    25152483
    25162484      integer iy,im,id,ih,in
     2485
     2486      print*,'nlev_armcu',nlev_armcu
    25172487
    25182488      open(21,file=trim(fich_armcu),form='formatted')
    … …  
    25292499
    25302500  223 format(5i3,5f8.3)
    2531   226 format(f7.1,1x,10f8.2)
    2532   227 format(f7.1,1x,1p,4e11.3)
    2533   230 format(6f9.3,4e11.3)
    25342501
    25352502          return
    … …  
    25712538       real hq_mod(llm),vq_mod(llm)
    25722539 
    2573        integer l,k,k1,k2,kp
    2574        real aa,frac,frac1,frac2,fact
     2540       integer l,k,k1,k2
     2541       real frac,frac1,frac2,fact
    25752542 
    25762543       do l = 1, llm
    … …  
    26842651        real div_prof,ts_prof,ug_prof,vg_prof,ufa_prof,vfa_prof
    26852652! local:
    2686         integer it_astex1, it_astex2,k
     2653        integer it_astex1, it_astex2
    26872654        real timeit,time_astex1,time_astex2,frac
    26882655
    … …  
    29852952
    29862953        integer, parameter :: ilesfile=1
    2987         integer :: ierr,irad,imax,jtot,k,itrac,nt1,nt2
    2988         logical :: lmoist,lcoriol,ltimedep
    2989         real :: xsize,ysize
    2990         real :: ustin,wsvsurf,timerad
    2991         character(80) :: chmess
     2954        integer :: ierr,k,itrac,nt1,nt2
    29922955
    29932956        if(.not.(llesread)) return
    … …  
    30433006      implicit none
    30443007
    3045         integer nlev_max,kmax,kmax2
     3008        integer nlev_max,kmax
    30463009        logical :: llesread = .true.
    30473010
    … …  
    30523015
    30533016        integer, parameter :: ilesfile=1
    3054         integer :: ierr,irad,imax,jtot,k
    3055         logical :: lmoist,lcoriol,ltimedep
    3056         real :: xsize,ysize
    3057         real :: ustin,wsvsurf,timerad
    3058         character(80) :: chmess
     3017        integer :: k,ierr
    30593018
    30603019        if(.not.(llesread)) return
    … …  
    30793038      implicit none
    30803039
    3081         integer nlev_max,kmax,kmax2
     3040        integer nlev_max,kmax
    30823041        logical :: llesread = .true.
    30833042
    … …  
    30883047
    30893048        integer, parameter :: ilesfile=1
    3090         integer :: ierr,irad,imax,jtot,k
    3091         logical :: lmoist,lcoriol,ltimedep
    3092         real :: xsize,ysize
    3093         real :: ustin,wsvsurf,timerad
    3094         character(80) :: chmess
     3049        integer :: ierr,k
    30953050
    30963051        if(.not.(llesread)) return
    … …  
    31173072      implicit none
    31183073
    3119         integer nlev_max,kmax,kmax2
     3074        integer nlev_max,kmax
    31203075        logical :: llesread = .true.
    31213076
    … …  
    31273082        integer, parameter :: ilesfile=1
    31283083        integer, parameter :: ifile=2
    3129         integer :: ierr,irad,imax,jtot,k
    3130         logical :: lmoist,lcoriol,ltimedep
    3131         real :: xsize,ysize
    3132         real :: ustin,wsvsurf,timerad
    3133         character(80) :: chmess
     3084        integer :: ierr,jtot,k
    31343085
    31353086        if(.not.(llesread)) return
    … …  
    31823133
    31833134      integer ntime,nlevel
    3184       integer l,k
    31853135      character*80 :: fich_amma
    3186       real*8 time(ntime)
    3187       real*8 zz(nlevel)
     3136      real*8 zz(nlevel)
    31883137
    31893138      real*8 temp(nlevel),pp(nlevel)
    … …  
    31923141      real*8 dw(nlevel,ntime)
    31933142      real*8 dt(nlevel,ntime)
    3194       real*8 dq(nlevel,ntime)   
     3143      real*8 dq(nlevel,ntime)
    31953144      real*8 flat(ntime),sens(ntime)
    31963145
    … …  
    35223471
    35233472      integer ntime,nlevel
    3524       integer l,k
    35253473      character*80 :: fich_fire
    3526       real*8 time(ntime)
    3527       real*8 zz(nlevel)
     3474      real*8 zz(nlevel)
    35283475
    35293476      real*8 thl(nlevel)
    … …  
    35323479      real*8 ug(nlevel,ntime),vg(nlevel,ntime),wls(nlevel,ntime)
    35333480      real*8 dqtdx(nlevel,ntime),dqtdy(nlevel,ntime)
    3534       real*8 dqtdt(nlevel,ntime),thl_rad(nlevel,ntime) 
     3481      real*8 dqtdt(nlevel,ntime),thl_rad(nlevel,ntime)
    35353482
    35363483      integer nid, ierr

  • LMDZ5/trunk/libf/phy1d/1D_decl_cases.h

    r1907 r1960  
    1616        real    sec_print
    1717!!
    18         integer nn
    19         integer it_toga1, it_toga2
    20         real time_toga1,time_toga2
    21 
    2218        real ts_toga(nt_toga)
    2319        real plev_toga(nlev_toga,nt_toga),w_toga(nlev_toga,nt_toga)
    … …  
    3430        real hq_prof(nlev_toga),vq_prof(nlev_toga)
    3531
    36         real plev_mod(llm),w_mod(llm), t_mod(llm),q_mod(llm)
     32        real w_mod(llm), t_mod(llm),q_mod(llm)
    3733        real u_mod(llm),v_mod(llm), ht_mod(llm),vt_mod(llm)
    3834        real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm)
    … …  
    8783        character*80 :: fich_amma
    8884! Option du cas AMMA ou on impose la discretisation verticale (Ap,Bp)
    89         logical  :: fixe_disvert=.true.
    9085        integer nlev_amma, nt_amma
    9186!       parameter (nlev_amma=29, nt_amma=48)  ! Fleur, juillet 2012
    … …  
    10499!profils initiaux:
    105100        real plev_amma(nlev_amma)
    106         real tv_amma(nlev_amma),rho_amma(nlev_amma)
    107         real thv_amma(nlev_amma)
    108101       
    109102        real z_amma(nlev_amma)
    … …  
    111104        real u_amma(nlev_amma)
    112105        real v_amma(nlev_amma)
    113 
    114         real thvsurf_amma,tvsurf_amma,rhosurf_amma,thsurf
    115106
    116107        real th_ammai(nlev_amma),q_ammai(nlev_amma)
    … …  
    130121
    131122!champs interpoles
    132         real plev_profamma(nlev_amma),vitw_profamma(nlev_amma)
     123        real vitw_profamma(nlev_amma)
    133124        real ht_profamma(nlev_amma)
    134125        real hq_profamma(nlev_amma)
    … …  
    148139        integer year_ini_fire, day_ini_fire, mth_ini_fire
    149140        real heure_ini_fire
    150         real day_ju_ini_fire   ! Julian day of fire first day
    151141        parameter (year_ini_fire=1987)
    152142        parameter (mth_ini_fire=7)
    … …  
    154144        parameter (heure_ini_fire=0.) !0h en secondes
    155145
    156 !profils initiaux:
    157         real z_fire(nlev_fire)
    158         real thl_fire(nlev_fire),qt_fire(nlev_fire)
    159         real u_fire(nlev_fire), v_fire(nlev_fire)
    160         real tke_fire(nlev_fire)
    161        
    162 !forcings
    163         real ugeo_fire(nlev_fire),vgeo_fire(nlev_fire)
    164         real wls_fire(nlev_fire),dqtdx_fire(nlev_fire)
    165         real dqtdy_fire(nlev_fire)
    166         real dqtdt_fire(nlev_fire),thl_rad_fire(nlev_fire)
    167          
    168         real ugeo_mod(llm),vgeo_mod(llm),wls_mod(llm)
    169         real dqtdx_mod(llm),dqtdy_mod(llm),dqtdt_mod(llm)
    170         real thl_rad_mod(llm)
    171146!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    172147! Declarations specifiques au cas GCSSold
    … …  
    180155        real  Ts_gcssold
    181156        real  dtime_frcg
    182 cAl1     logical :: imp_fcg_gcssold
    183 c        logical :: ts_fcg_gcssold
    184 c        logical :: Tp_fcg_gcssold
    185157        logical :: Turb_fcg_gcssold
    186         common /turb_forcing/ dtime_frcg,
    187      $      Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold
     158
     159        common /turb_forcing/
     160     s  dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold
    188161!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    189162! Declarations specifiques au cas Arm_cu
    … …  
    206179        real adv_qt_armcu(nt_armcu)
    207180        real theta_mod(llm),rv_mod(llm),play_mod(llm)
    208         real d_t_dyn_ls(llm),d_q_dyn_ls(llm)
    209181! profc comme "profil armcu"
    210         real h_profc,play_profc,t_profc,th_profc,plev_profc
    211         real u_profc,v_profc,qv_profc,rv_profc
    212182       
    213183! forcages interpoles dans le temps
    … …  
    226196        parameter (mth_ini_sandu=7)
    227197        parameter (day_ini_sandu=15)  ! 196 = 15 juillet 2006
    228         real dt_sandu, tau_sandu
     198        real dt_sandu
    229199        logical  :: trouve_700=.true.
    230200        parameter (dt_sandu=6.*3600.)   ! forcages donnes ttes les 6 heures par ifa_sandu.txt
    231201!       parameter (tau_sandu=3600.)  ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa
    232202!!
    233         integer it_sandu1, it_sandu2
    234         real time_sandu1,time_sandu2
    235 
    236203        real ts_sandu(nt_sandu)
    237204! profs comme "profil sandu"
    … …  
    243210
    244211        real thl_mod(llm),omega_mod(llm),o3mmr_mod(llm),tke_mod(llm)
    245 ! pour relaxer u,v,thl et qt vers les profils initiaux au dessus de 700hPa
    246         real relax_u(llm),relax_v(llm),relax_thl(llm),relax_q(llm,2)
    247212!vertical advection computation
    248213        real d_t_z(llm), d_q_z(llm)
    … …  
    260225        parameter (mth_ini_astex=6)
    261226        parameter (day_ini_astex=13)  ! 165 = 13 juin 1992
    262         real dt_astex, tau_astex
     227        real dt_astex
    263228        parameter (dt_astex=3600.)    ! forcages donnes ttes les heures par ifa_astex.txt
    264         integer it_astex1, it_astex2
    265         real time_astex1,time_astex2
    266229        real ts_astex(nt_astex),div_astex(nt_astex),ug_astex(nt_astex)
    267230        real vg_astex(nt_astex),ufa_astex(nt_astex),vfa_astex(nt_astex)

  • LMDZ5/trunk/libf/phy1d/1Dconv.h

    r1907 r1960  
    2222      REAL hplaym(100) !pression en hPa milieux des couches Meso-NH
    2323
    24       integer i,j,k,ii,ll,in
    25       REAL tsol,qsol
     24      integer i,j,k,ll,in
    2625
    2726      CHARACTER*80 file_forctl,file_fordat
    2827
    29       COMMON/com1_phys_gcss/klev,play,JM,coef1,coef2
    30       COMMON/com2_phys_gcss/nblvlm,playm,hplaym
     28      COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev
     29      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
    3130
    3231c======================================================================
    … …  
    170169             hqTurbbef(i)=hqTurbaft(i)
    171170          enddo
    172                tsbef = tsaft
     171          tsbef = tsaft
    173172          timebef=pasprev*dt
    174173          timeaft=timebef+dt
    … …  
    213212               print*,'hqTurb_mes ',(hqTurb_mes(i),i=1,nblvlm)
    214213                  endif
    215                IF (ts_fcg) print*,'ts_subr', ts_subr
     214          IF (ts_fcg) print*,'ts_subr', ts_subr
    216215c*** on interpole les champs meso_NH sur les niveaux de pression***
    217216c*** gcm . on obtient le nouveau champ after                    ***
    … …  
    263262               hqTurb(ll)=hqTurbaft(ll)
    264263          enddo
    265                ts_subr = tsaft
     264          ts_subr = tsaft
    266265      else ! temps.ge.pasmax
    267266c*** on interpole sur les pas de temps de 10mn du gcm a partir   ***
    … …  
    282281             endif ! Turb_fcg
    283282         enddo
    284               ts_subr = ((timeaft-time)*tsbef + (time-timebef)*tsaft)/dt
     283         ts_subr = ((timeaft-time)*tsbef + (time-timebef)*tsaft)/dt
    285284       endif ! temps.ge.pasmax
    286285c
    … …  
    440439                endif ! Turb_fcg
    441440          enddo
    442                ts_subr = tsaft
     441          ts_subr = tsaft
    443442          close(99)
    444443          close(98)
    … …  
    505504      REAL hplaym(100)!pression en hecto-Pa des milieux de couche Meso-NH
    506505
    507       COMMON/com1_phys_gcss/klev,play,JM,coef1,coef2
    508       COMMON/com2_phys_gcss/nblvlm,playm,hplaym
    509 
    510       integer i,k,klevgcm
     506      COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev
     507      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
     508
     509      integer k,klevgcm
    511510      real playgcm(klevgcm) ! pression en milieu de couche du gcm
    512511      real psolgcm
    … …  
    577576      REAL playm(100)  !pression en Pa milieu de chaque couche Meso-NH
    578577      REAL hplaym(100) !pression en hPa des milieux de couche Meso-NH
    579       COMMON/com2_phys_gcss/nblvlm,playm,hplaym
    580 
    581       INTEGER i,lu,k,mlz,mlzh,j
     578      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
     579
     580      INTEGER i,lu,mlz,mlzh
    582581 
    583582      character*80 file_forctl
    … …  
    644643      real ts
    645644c
    646       INTEGER i, k
     645      INTEGER k
    647646c
    648647      LOGICAL imp_fcg,ts_fcg,Turb_fcg
    … …  
    725724      REAL hplaym(100)!pression en hPa des milieux de couche Meso-NH
    726725
    727       COMMON/com1_phys_gcss/klev,play,JM,coef1,coef2
    728       COMMON/com2_phys_gcss/nblvlm,playm,hplaym
     726      COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev
     727      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
    729728
    730729      REAL psol
    731730      REAL val
    732       INTEGER k, mlz, mlzh
     731      INTEGER k, mlz
    733732
    734733

  • LMDZ5/trunk/libf/phy1d/add_pbl_tend.F

    r1907 r1960  
    1919        real dtime_frcg
    2020        logical  Turb_fcg_gcssold
    21         common /turb_forcing/ dtime_frcg,
    22      $           Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold
     21        common /turb_forcing/
     22     s  dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold
     23
    2324! Arguments :
    2425!------------

  • LMDZ5/trunk/libf/phy1d/compar1d.h

    r1907 r1960  
    2727      logical :: ok_old_disvert
    2828
    29       common/com_par1d/forcing_type,nat_surf,tsurf,rugos,               &
     29      common/com_par1d/
     30     & nat_surf,tsurf,rugos,                                            &
    3031     & qsol,qsurf,psurf,zsurf,albedo,time,time_ini,xlat,xlon,airefi,    &
    3132     & wtsurf,wqsurf,restart_runoff,xagesno,qsolinp,zpicinp,            &
     33     & forcing_type,
    3234     & restart,ok_old_disvert
    3335

  • LMDZ5/trunk/libf/phy1d/lmdz1d.F

    r1948 r1960  
    8181      integer :: an
    8282 
    83 !
    84       real :: paire    = 1.     ! aire de la maille
    85 !**      common /flux_arp/fsens,flat,ok_flux_surf
    86 
    8783!---------------------------------------------------------------------
    8884!  Declarations related to forcing and initial profiles
    … …  
    9288        integer llm700,nq1,nq2
    9389        INTEGER, PARAMETER :: nlev_max=1000, nqmx=1000
    94         real timestep, frac, timeit
     90        real timestep, frac
    9591        real height(nlev_max),tttprof(nlev_max),qtprof(nlev_max),
    9692     .              uprof(nlev_max),vprof(nlev_max),e12prof(nlev_max),
    … …  
    10096     .              qprof(nlev_max,nqmx)
    10197
    102         real    :: fff
    10398c        integer :: forcing_type
    10499        logical :: forcing_les     = .false.
    … …  
    138133      real :: pzero=1.e5
    139134      real :: play (llm),zlay (llm),sig_s(llm),plev(llm+1)
    140       real :: playd(llm),zlayd(llm),ap_amma(llm+1),bp_amma(llm+1),poub
    141135
    142136!---------------------------------------------------------------------
    … …  
    144138!---------------------------------------------------------------------
    145139
    146       integer :: iq
    147140      real :: phi(llm)
    148141      real :: teta(llm),tetal(llm),temp(llm),u(llm),v(llm),w(llm)
    … …  
    152145      real :: sfdt, cfdt
    153146      real :: du_phys(llm),dv_phys(llm),dt_phys(llm)
    154       real :: du_dyn(llm),dv_dyn(llm),dt_dyn(llm)
    155       real :: dt_cooling(llm),d_t_cool(llm),d_th_adv(llm)
    156       real :: dq_cooling(llm),d_q_cool(llm)
    157       real :: tmpvar(llm)
     147      real :: dt_dyn(llm)
     148      real :: dt_cooling(llm),d_th_adv(llm)
    158149      real :: alpha
     150      real :: ttt
    159151
    160152      REAL, ALLOCATABLE, DIMENSION(:,:):: q
    … …  
    205197!  Fichiers et d'autres variables
    206198!---------------------------------------------------------------------
    207       real ttt,bow,q1
    208       integer :: ierr,k,l,i,it=1,mxcalc
     199      integer :: k,l,i,it=1,mxcalc
    209200      integer jjmp1
    210201      parameter (jjmp1=jjm+1-1/jjm)
    … …  
    231222!---------------------------------------------------------------------
    232223cAl1
    233         call conf_unicol(99)
     224        call conf_unicol
    234225cAl1 moves this gcssold var from common fcg_gcssold to
    235226        Turb_fcg_gcssold = xTurb_fcg_gcssold
    … …  
    375366      itau_phy = 0
    376367      call ymds2ju(annee_ref,mois,day_ref,heure,day)
    377       day_ini = day
     368      day_ini = int(day)
    378369      day_end = day_ini + fnday
    379370
    … …  
    859850!
    860851        du_age(1:mxcalc)= -2.*sfdt/timestep*
    861      :          (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) -
    862      :           cfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
     852     s          (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) -
     853     s           cfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
    863854!!     : fcoriolis*(v(1:mxcalc)-vg(1:mxcalc))
    864855!
    865856       dv_age(1:mxcalc)= -2.*sfdt/timestep*
    866      :          (cfdt*(u(1:mxcalc)-ug(1:mxcalc)) +
    867      :           sfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
     857     s          (cfdt*(u(1:mxcalc)-ug(1:mxcalc)) +
     858     s           sfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
    868859!!     : -fcoriolis*(u(1:mxcalc)-ug(1:mxcalc))
    869860!
    … …  
    879870!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    880871        u(1:mxcalc)=u(1:mxcalc) + timestep*(
    881      :              du_phys(1:mxcalc)
    882      :             +du_age(1:mxcalc) )           
     872     s              du_phys(1:mxcalc)
     873     s             +du_age(1:mxcalc) )           
    883874        v(1:mxcalc)=v(1:mxcalc) + timestep*(
    884      :               dv_phys(1:mxcalc)
    885      :             +dv_age(1:mxcalc) )
     875     s              dv_phys(1:mxcalc)
     876     s             +dv_age(1:mxcalc) )
    886877        q(1:mxcalc,:)=q(1:mxcalc,:)+timestep*(
    887      :                dq(1:mxcalc,:)
    888      :               +d_q_adv(1:mxcalc,:) )
     878     s                dq(1:mxcalc,:)
     879     s               +d_q_adv(1:mxcalc,:) )
    889880
    890881        if (prt_level.ge.1) then
Note: See TracChangeset for help on using the changeset viewer.