Changeset 180


Ignore:
Timestamp:
02/13/18 16:23:07 (6 years ago)
Author:
aquiquet
Message:

Cleaning-up: unused variables removed

Location:
trunk/SOURCES
Files:
18 edited

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90

    r144 r180  
    6969  character(len=100) :: temp_annual_file   ! temperature annuelles 
    7070  character(len=100) :: file_smb_anom      !> nom du fichier dans lequel il y a l'anomalie smb 
    71   character(len=100) :: file_bmelt_anom    !> nom du fichier dans lequel il y a l'anomalie bmelt 
    7271 
    7372  integer            :: err                ! recuperation d'erreur 
    74   integer            :: i,j,k 
     73  integer            :: i 
    7574 
    7675  !aurel for Tafor: 

  • trunk/SOURCES/Ant40_files/bmelt-ant-regions-oce_mod.f90

    r34 r180  
    5353  integer :: nbr_pts ! nbr de points utilises pour calcul temperature oceanique sur points sans valeurs 
    5454  logical,dimension(nx,ny) :: mask_oce1,mask_oce2 ! masque des zones avec temp ocean (true si existe valeur) 
    55   real :: flpath ! distance de recherche (rayon pour les temp oceanique non presentes dans le fichier) 
    5655  integer :: n,imin,imax,jmin,jmax 
    5756  logical :: test ! pour stoper boucle 

  • trunk/SOURCES/Draggings_modules/dragging_prescr_beta_mod.f90

    r168 r180  
    292292    if (itracebug.eq.1)  call tracebug(' apres criteres flot') 
    293293 
    294     if (itracebug.eq.1) write(num_tracebug,*) 'gzmx', gzmx(127,330)  
     294    if (itracebug.eq.1) write(num_tracebug,*) 'gzmx', gzmx(127,127)  
    295295 
    296296    !--------- autres criteres 
    … …  
    326326 
    327327    if (itracebug.eq.1)  call tracebug(' apres autres criteres ') 
    328     if (itracebug.eq.1) write(num_tracebug,*) 'gzmx', gzmx(127,330)  
     328    if (itracebug.eq.1) write(num_tracebug,*) 'gzmx', gzmx(127,127)  
    329329 
    330330    flgzmx(:,:) = flgzmx(:,:) .or. gzmx(:,:) 

  • trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_mod.f90

    r55 r180  
    323323  implicit none 
    324324  integer             :: k_snap               ! pour calculer les indices de temps 
    325   real                :: time_bis             ! pour repliquer les dernieres annees 
    326325  integer             :: k 
    327326 

  • trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

    r174 r180  
    6464real,dimension(:),allocatable :: bpert          !> basal melt index perturbation 
    6565real :: coefT                    !> pour modifier l'amplitude de la perturb. T 
    66 character(len=120) :: filforc    !> nom du fichier forcage 
     66character(len=150) :: filforc    !> nom du fichier forcage 
    6767integer :: pertsmb               !> boolean: do we modify the smb? 
    6868real    :: rapsmb                !> if we modify the smb this is the equivalent of rappact 
    … …  
    102102  !aurel for Tafor: 
    103103  character(len=8) :: control      !label to check clim. forc. file (filin) is usable 
    104   character(len=80):: filin 
     104  character(len=100):: filin 
    105105 
    106106  namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file 
    … …  
    300300   
    301301  tpert(:)=tpert(:)*coefT 
    302   bpert(:)=1.+bpert(:)*coefbmb 
     302  bpert(:)=max(1.+bpert(:)*coefbmb,0.01) !freezing is not allowed 
    303303 
    304304   
    … …  
    444444        !afq -- as in Pollard and Deconto, melt is more efficient at high temperature... -> 
    445445        if (coefbmshelf.gt.1.) then 
    446            coefbmshelf = 1. + ( coefbmshelf - 1. ) * 5. !afq.... tuning for the degla?... 
     446           coefbmshelf = 1. + ( coefbmshelf - 1. ) * 1.5 !afq.... tuning for the degla?... 
    447447        end if 
    448448         
    … …  
    467467  implicit none 
    468468  integer             :: k_snap               ! pour calculer les indices de temps 
    469   real                :: time_bis             ! pour repliquer les dernieres annees 
    470469  integer             :: k 
    471470 

  • trunk/SOURCES/New-remplimat/diagno-L2_mod.f90

    r176 r180  
    350350!$  logical :: paral 
    351351 
    352   integer :: t1,t2,ir 
    353   real                           :: temps, t_cpu_0, t_cpu_1, t_cpu, norme   
    354  
    355352if (itracebug.eq.1)  call tracebug(' Calc pvi') 
    356353 

  • trunk/SOURCES/New-remplimat/eq_ellipt_sgbsv_mod-0.2.f90

    r72 r180  
    8888! ------------------------------------------------------------------------------- 
    8989 
    90 subroutine resol_ellipt(nx1,nx2,ny1,ny2,uxprec,uyprec,uxnew,uynew,imx,imy,ifail) 
     90subroutine resol_ellipt(nx1,nx2,ny1,ny2,uxprec,uyprec,uxnew,uynew,ifail) 
    9191 
    9292!  prepare la matrice mmat (version "couchee" de L2) 
    … …  
    101101!      uynew(n1,n2) uynew resultat de cette iteration  
    102102! 
    103 !      imx(n1,n2) masque pour imposer les vitesses ou leur dérivee     
    104 !      imy(n1,n2) masque pour imposer les vitesses ou leur dérivée           
    105103!--------------------------------------------------------------------------- 
    106104 
    … …  
    114112real, dimension(nx1:nx2,ny1:ny2),intent(in)  :: uxprec    ! vitesse en entree routine 
    115113real, dimension(nx1:nx2,ny1:ny2),intent(in)  :: uyprec    ! vitesse en entree routine 
    116  
    117 integer, dimension(nx1:nx2,ny1:ny2),intent(in)  :: imx    ! masque en entree routine 
    118 integer, dimension(nx1:nx2,ny1:ny2),intent(in)  :: imy    ! masque en entree routine 
    119  
    120114 
    121115real, dimension(nx1:nx2,ny1:ny2),intent(out) :: uxnew     ! vitesse en sortie de la routine 

  • trunk/SOURCES/New-remplimat/remplimat-shelves-tabTu.f90

    r169 r180  
    9292real    :: scal          ! pour le conditionnement (diagonale=1) 
    9393 
    94 character(len=80) :: fileout ! pour sortie fichier verification 
    95  
    96  
    9794! pour les fronts, on suppose que l'epaisseur est celle du noeud amont 
    9895! dans opposx et opposy 
    … …  
    326323                    uxnew(nx1:nx2,ny1:ny2),  uynew(nx1:nx2,ny1:ny2),     &   ! nouvelles vitesses 
    327324 
    328                     imx(nx1:nx2,ny1:ny2),  imy(nx1:nx2,ny1:ny2),ifail_L2)    ! masques 
     325                    ifail_L2)                                                ! masques 
    329326 
    330327if (itracebug.eq.1) call tracebug ('apres subroutine resol_ellipt') 

  • trunk/SOURCES/calving_frange.f90

    r179 r180  
    114114  subroutine calving 
    115115 
    116     integer :: I_did_something, m ! pour la boucle sur le calving 
     116    integer :: I_did_something  ! pour la boucle sur le calving 
    117117 
    118118    ! initialisation calving 

  • trunk/SOURCES/courbures.f90

    r4 r180  
    3939parameter (npts=60000)      !< nombre de points a traiter (nx*ny) 
    4040 
    41  
    42 real,dimension(mnode) :: Z_sg               !< altitude des points de la sous-grille 
    4341real :: dx_ls                               !< pas de la grille dans le moindre carre 
    44 real,dimension(ncoef) :: acoef              !< matrice de coef pour construction  
    45 integer :: i,j,k,l,pp,iz,jz,iz1,jz1         !< compteurs 
     42integer :: i,j,k,pp,iz,jz,iz1,jz1         !< compteurs 
    4643real :: x,y                                 !< coordonnees dans la sous-grille 
    4744real :: aaa                                 !< pour le calcul de l'azimut 

  • trunk/SOURCES/eaubasale-0.5_mod.f90

    r148 r180  
    112112kond(:,:) = kond(:,:)*SECYEAR  
    113113kondmax   = 1.*SECYEAR  
    114 keffmax=kondmax*hmax_wat 
    115114hdotwater(:,:)=0. 
    116115 
    … …  
    274273!$OMP END PARALLEL 
    275274!!$OMP SINGLE 
    276            call relaxation_waterdif(nxlocal,nylocal,dt,dx,vieuxhwater,limit_hw,klimit,bmelt_w,infiltr,pgx,pgy,keff,keffmax,hwater) 
     275           call relaxation_waterdif(nxlocal,nylocal,dt,dx,vieuxhwater,limit_hw,klimit,bmelt_w,infiltr,pgx,pgy,keff,hwater) 
    277276!!$OMP END SINGLE 
    278277        else 

  • trunk/SOURCES/flottab2-0.7.f90

    r169 r180  
    640640!      integer :: mask_nb = 4 
    641641integer,parameter :: mask_nb = 2   ! version ou on ne compte pas les diagonales 
    642 integer :: vartemp       ! variable temporaire pour reorganiser compt 
    643642!     integer,dimension(mask_nb) :: mask 
    644643integer,dimension(mask_nb) :: mask 

  • trunk/SOURCES/furst_schoof_mod.f90

    r176 r180  
    5454  integer,dimension(nx,ny) :: countx, county   ! how often do we modify ux/uy 
    5555  real :: phi_prescr 
    56   real :: archim 
    5756  real :: toutpetit = 1e-6 
    5857  real :: denom, prodscal 

  • trunk/SOURCES/litho-0.4.f90

    r74 r180  
    4545  implicit none 
    4646 
    47   INTEGER :: II,SOM1,SOM2 
     47  INTEGER :: SOM1,SOM2 
    4848  REAL, dimension(:,:), allocatable :: WLOC 
    4949 

  • trunk/SOURCES/out_cptr_mod.f90

    r145 r180  
    608608!> 
    609609 
    610 subroutine symetry_cptr(iaxe,jaxe) 
     610subroutine symetry_cptr(jaxe) 
    611611  use module3D_phy 
    612612  use tracer_vars 
    613613  implicit none 
    614   integer :: iaxe  
    615614  integer :: jaxe 
    616615  integer :: jsym 
    617   integer :: isym 
    618616 
    619617  !symetrique par rapport à jaxe la référence est le bas. Pour l'instant seulement lui 

  • trunk/SOURCES/relaxation_water_diffusion.f90

    r76 r180  
    1717 
    1818 
    19 subroutine relaxation_waterdif(NXX,NYY,DT,DX,vieuxHWATER,limit_hw,klimit,BMELT,INFILTR,PGMX,PGMY,KOND,KONDMAX,HWATER) 
     19subroutine relaxation_waterdif(NXX,NYY,DT,DX,vieuxHWATER,limit_hw,klimit,BMELT,INFILTR,PGMX,PGMY,KOND,HWATER) 
    2020 
    2121  !$ USE OMP_LIB 
    … …  
    3030    REAL,    intent(in) ::  DX            !< pas en x 
    3131    REAL,    intent(in) :: INFILTR        !< basal infiltration (lose of water)  
    32     REAL,    intent(in) :: KONDMAX        !< maximum hydaulic conductivity (outside ice sheet) 
    3332 
    3433    REAL,dimension(NXX,NYY), intent(in) :: limit_hw    !< conditions aux limites 

  • trunk/SOURCES/steps_time_loop.f90

    r165 r180  
    2828  implicit none 
    2929  integer :: nt_init = 0      ! number of loops for initialisation of ice thickness 
    30   integer :: iter_visco       ! number of iterations for ssa viscosity 
    3130 
    3231  if (itracebug.eq.1)  call tracebug('Entree dans step_time_loop ') 

  • trunk/SOURCES/steps_time_loop_avec_iterbeta.f90

    r161 r180  
    2828  implicit none 
    2929  integer :: nt_init = 0      ! number of loops for initialisation of ice thickness 
    30   integer :: iter_visco       ! number of iterations for ssa viscosity 
    3130 
    3231  if (itracebug.eq.1)  call tracebug('Entree dans step_time_loop ') 
Note: See TracChangeset for help on using the changeset viewer.