Ignore:
Timestamp:
Feb 7, 2014, 5:20:25 PM (11 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.

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