Changeset 65 for trunk/SOURCES


Ignore:
Timestamp:
05/18/16 17:11:29 (8 years ago)
Author:
dumas
Message:

Deleting unused variables and move old sources

Location:
trunk/SOURCES
Files:
29 edited
2 moved

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/Fichiers-parametres/Makefile.tof-lsce3130.inc

    r64 r65  
    3939        IFORT= ifort 
    4040 
    41         ARITHMi    = -O2 -fp-model precise  -heap-arrays -traceback -mcmodel=medium -diag-disable warn #  -warn -traceback   -CB  -g  # options pour une meilleure arithmetique 
     41        ARITHMi    = -O2 -fp-model precise  -heap-arrays -traceback -mcmodel=medium #  -warn -traceback   -CB  -g  # options pour une meilleure arithmetique 
    4242                                                                      # (normalement reproductible) 
    4343    ifeq ($(debug), 1) 
    4444                ARITHM    = $(ARITHMi) -CB -g -traceback -warn all 
    4545        else 
    46                 ARITHM    = $(ARITHMi) 
     46                ARITHM    = $(ARITHMi) -diag-disable warn 
    4747        endif 
    4848                 

  • trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

    r59 r65  
    379379  implicit none 
    380380 
    381   integer i,j,ift 
     381  integer i,ift 
    382382   
    383383  select case (massb_time) 

  • trunk/SOURCES/Hemin40_files/output_hemin40_mod.f90

    r13 r65  
    157157     integer inp(13) ! Surface posee (nb de noeuds) 
    158158     integer inf(13) ! surface flottante 
    159       integer npab, npcal 
    160  
    161       REAL long, lat, rd, cx, cy, cxmin, cymin, r2 
     159 
    162160!Variables pour sommer       
    163161      real isvol(13),isvolf(13) ! volume posé et flottants 
    … …  
    165163      REAL ISABLBORD(13),ABLATOT(13),TACC(13),TBM(13) 
    166164      REAL ITJJA(13) 
    167       REAL ABLAMEAN 
     165!      REAL ABLAMEAN 
    168166!moyennes utilisées en output 
    169167      REAL HMAX_(13) , HMEAN_(13)  
    … …  
    178176 
    179177      REAL HDOTMEAN_G 
    180       REAL ABLA(NX,NY) 
     178!      REAL ABLA(NX,NY) 
    181179      REAL DELTAVOL 
    182180 
    … …  
    384382905 format(f9.1,1x, e11.4,1x,e11.5, 1x,i5, 1x,e10.4,1x,i5 , 6(1x,e12.5), & 
    385383         12(2(1x,e10.4,1x,i5), 9(1x,e11.4) ) )  
    386 900   format(i5,1x,f9.1,1x,f5.1,1x,e10.5,1x,i4,1x,i4,1x,f6.2,1X, & 
    387              f6.2,1X,f6.2,1X,f6.2,1X,i5,8(1x,f6.2)) 
    388 903   format(f9.1,1x,f4.1,1x,f5.1,1x,e10.5,1x,i5,1x,i4,1x,i5,1X, & 
    389              f7.3,1x,f6.2,1X,i6,4(1x,e7.1),1X,f6.2,1x,f7.3) 
    390384  
    391 904  format(F9.1,1x,e10.5,1x,e11.5,1x,i5,6(1x,e12.6),3(1x,e10.5,        & 
    392          1x,i5,6(1X,e12.6)),1X,i4,1X,i5,1x,f7.3,1x,                      & 
    393          f6.2,1x,i6,2(1x,e7.1),2(1x,f8.4),1X,f6.2,2(1x,i5),1x,e7.1,      & 
    394          8(1x,e12.6),4(1x,f15.7),9(1x,e10.5,1x,i5,2(1x,e12.6),1x,f15.7), & 
    395          26(1x,f8.3)) 
    396385     
    397386940   format('%%%% ',a,'   time=',f8.0,' %%%%') 

  • trunk/SOURCES/Makefile

    r64 r65  
     1 
    12# Makefile pour construire les executables grisli en fonction de la machin 
    23# Pour compiler : make -f Makefile-template Grice2sea 
    … …  
    1415# Choice for librairies (mkl_c = 1 = MKL,  = 0 = BLAS) 
    1516# NOTA: mkl_c = 1 with ifort = 0 ignores the MKL and uses BLAS anyhow 
    16 mkl_c ?= 0 
     17mkl_c ?= 1 
    1718 
    1819#ifeq ($(ifort),1) 

  • trunk/SOURCES/Makefile.grisli.inc

    r64 r65  
    7575        calving_frange.o no_calving.o no_lakes.o \ 
    7676        out_profile_mod.o printtable_mod.o mix-SIA-L1_mod.o \ 
    77         relaxation_mod-0.3.o relaxation_water_mod-0.4.o relaxation_water_diffusion.o \ 
     77        relaxation_water_diffusion.o \ 
    7878        prescribe-H-i2s_mod.o  \ 
    7979        resol_adv_diff_2D-sept2009.o  \ 
    … …  
    9191        calving_frange.o no_calving.o no_lakes.o \ 
    9292        out_profile_mod.o printtable_mod.o mix-SIA-L1_mod.o \ 
    93         relaxation_mod-0.3.o relaxation_water_mod-0.4.o relaxation_water_diffusion.o \ 
     93        relaxation_water_diffusion.o \ 
    9494        prescribe-H-i2s_mod.o   \ 
    9595        resol_adv_diff_2D-sept2009.o  \ 

  • trunk/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90

    r60 r65  
    157157    integer :: num_file=22 
    158158    integer :: i1,i2,i3 
    159     integer :: i,j,k,l      
     159    integer :: i,j,k      
    160160    character(len=20) :: name1 
    161161    character(len=20) :: name2 
    … …  
    383383    integer   :: ifin 
    384384    integer   :: npr 
    385     integer  :: i,j,k !< indices de travail 
     385    integer  :: i !< indices de travail 
    386386 
    387387    if (itracebug.eq.1)  call tracebug(' Entree dans routine testsort_time_ncdf') 

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

    r4 r65  
    8686subroutine diagnoshelf !      Resolution numerique des equations diagnostiques 
    8787 
    88   integer :: itest 
    8988 
    9089  if (itracebug.eq.1)  call tracebug(' Entree dans diagnoshelf') 

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

    r59 r65  
    8787!-------------------- 
    8888 
    89 real    :: epsilon 
    9089real    :: dx2=dx*dx     ! variable de travail 
    9190real    :: beta          ! pour le frottement 
    92 real    :: moteur        ! pour le frottement 
    9391real    :: scal          ! pour le conditionnement (diagonale=1) 
     92 
     93character(len=80) :: fileout ! pour sortie fichier verification 
    9494 
    9595 
    … …  
    233233 
    234234if (iter_beta.eq.1) then 
    235    call calc_beta(Vcol_x, Vcol_y)  
    236  
    237    call write_datfile3(nx,ny,betamx,betamy,beta_centre,'beta-estime.dat') 
    238    if (itracebug.eq.1)  call tracebug(' remplimat apres calc_beta') 
    239    iter_beta = 0             ! pour faire un calcul de vitesse le coup d'apres 
    240    corr_def = .true. 
    241    goto 888 
     235        call calc_beta(Vcol_x, Vcol_y)  
     236        fileout='beta-estime.dat' 
     237        call write_datfile3(nx,ny,betamx,betamy,beta_centre,fileout) 
     238        if (itracebug.eq.1)  call tracebug(' remplimat apres calc_beta') 
     239        iter_beta = 0             ! pour faire un calcul de vitesse le coup d'apres 
     240        corr_def = .true. 
     241        goto 888 
    242242 
    243243else if ((iter_beta.eq.0).and.(corr_def)) then 

  • trunk/SOURCES/Recul_force_grounding_line/toy_retreat_mod.f90

    r22 r65  
    115115  use declar_toy_retreat 
    116116  implicit none 
    117   integer                           :: som_voisins  !< pour des tests 
     117!  integer                           :: som_voisins  !< pour des tests 
    118118 
    119119 
    … …  
    637637  implicit none 
    638638 
    639   integer                           :: som_voisins  !< pour des tests 
     639!  integer                           :: som_voisins  !< pour des tests 
    640640  real                              :: Hf           !< Hf = 0 si socle > 0 sinon H_float 
    641641  real                              :: dyy         ! variable de travail  longueur de la maille en diagonale        
    … …  
    752752  implicit none 
    753753 
    754   integer                           :: som_voisins  !< pour des tests 
     754!  integer                           :: som_voisins  !< pour des tests 
    755755 
    756756  ! les nouveaux points a traiter 
    … …  
    873873integer    :: num_job           ! numero du job 
    874874integer    :: lenrun            ! longueur du run name 
    875 integer    :: lenb              ! longueur du nom de fichier  
    876 integer    :: pos               ! position d'un sous_string 
     875!integer    :: lenb              ! longueur du nom de fichier  
     876!integer    :: pos               ! position d'un sous_string 
    877877 
    878878real*8, dimension(:,:),   pointer  :: tab         !< tableau 2d real ecrit dans le fichier 
    879879character(len=4)                   :: char_num    ! pour la lecture 
    880880character(len=80)                  :: file_B_test !< file name of the bedrock map to test topo instability 
    881 character(len=80)                  :: file_voisin !< file name of the bedrock map min or max to test topo instability 
     881!character(len=80)                  :: file_voisin !< file name of the bedrock map min or max to test topo instability 
    882882 
    883883 
    … …  
    14451445integer,intent(in)                 :: ii,jj         ! indice du point considere 
    14461446 
    1447 real                               :: dH_x          ! variable de travail 
    1448 real                               :: dH_y          ! variable de travail 
     1447!real                               :: dH_x          ! variable de travail 
     1448!real                               :: dH_y          ! variable de travail 
    14491449 
    14501450DelHdt_sanity(ii,jj) = 5000. 

  • trunk/SOURCES/climat-forcage-insolation_mod.f90

    r21 r65  
    149149 
    150150  implicit none 
    151   integer :: mo,ti,tj   
    152   character(len=8) :: control   ! label to check clim. forc. file (filin) is usable 
    153   integer :: l                  ! In snapshot files:the first column is the mask, read but not used 
    154   real :: T_surf_ref            ! variable de travail calcul temp a l'instant t a la surface S 
     151   
    155152  integer :: intr 
    156153  integer :: igtr 
    157154  integer :: ictr 
    158   integer :: k,d1,d2 
     155  integer :: d1,d2 
    159156  integer :: i,j 
    160   character(len=100) :: file_ncdf      !< fichier netcdf  
     157!  character(len=100) :: file_ncdf      !< fichier netcdf  
    161158  real*8, dimension(:,:,:), pointer :: data_3D => null() ! donnees lues dans le netcdf 
    162159  real*8, dimension(:,:),pointer :: data_2D => null()    ! donnees lues dans le netcdf 
    … …  
    450447  write(num_rep_42,'(A,2(f7.1,","))')   'palier_ice = ', palier_ice(:,:) 
    451448  write(num_rep_42,'(A,A)')   'surf_ice = ', surf_ice 
    452   write(num_rep_42,'(A,2(f3.1,","))')   'palier_CO2 = ', palier_CO2(:) 
     449  write(num_rep_42,'(A,2(f6.1,","))')   'palier_CO2 = ', palier_CO2(:) 
    453450  write(num_rep_42,'(A,A)')   'orb_file = ', orb_file 
    454451  write(num_rep_42,'(A,A)')   'co2_file = ', co2_file 
    … …  
    494491implicit none 
    495492 
    496 real COEFT,COEFP       ! 
    497493!integer l              ! dumm index for loops on snapshots files  l=ITR,NTR-1 
    498494!cdc integer itr            ! index of the current snapshot file (change with time) 
    … …  
    517513!real,dimension(nx,ny) :: coeff1,coeff2 
    518514real,dimension(nx,ny) :: coeff2 
    519 integer :: k1,k2,m1,m2 
     515integer :: m1,m2 
    520516integer,dimension(nx,ny) :: Mask,Mask2 
    521517integer :: update_H 
    … …  
    523519 
    524520!Definition des variables pour le calcul de la fonction de CO2 
    525 double precision :: coefCO2 
     521!double precision :: coefCO2 
    526522 
    527523!***************************************************************************** 

  • trunk/SOURCES/climat-forcage-insolation_mod_oneway.f90

    r18 r65  
    8585 
    8686  implicit none 
    87   integer :: mo,ti,tj   
    88   character(len=8) :: control   ! label to check clim. forc. file (filin) is usable 
    89   integer :: l                  ! In snapshot files:the first column is the mask, read but not used 
    90   real :: T_surf_ref            ! variable de travail calcul temp a l'instant t a la surface S 
     87 
    9188  integer :: intr 
    9289  integer :: igtr 
    9390  integer :: ictr 
    94   integer :: i,j 
    95   character(len=100) :: file_ncdf      !< fichier netcdf issue des fichiers .dat 
     91!  character(len=100) :: file_ncdf      !< fichier netcdf issue des fichiers .dat 
    9692  real*8, dimension(:,:,:), pointer :: data_3D => null() ! donnees lues dans le netcdf 
    9793  real*8, dimension(:,:),pointer :: data_2D => null()    ! donnees lues dans le netcdf 
    … …  
    231227implicit none 
    232228 
    233 real COEFT,COEFP       ! 
    234229!integer l              ! dumm index for loops on snapshots files  l=ITR,NTR-1 
    235230!cdc integer itr            ! index of the current snapshot file (change with time) 
    236231 
    237232integer mo 
    238 integer :: ictr 
    239 integer :: igtr 
    240233integer :: i,j 
    241234 

  • trunk/SOURCES/climat-perturb_mod-0.4.f90

    r27 r65  
    197197!  use module3d_phy 
    198198  implicit none 
    199   integer :: i,j,ift 
     199  integer :: i,ift 
    200200 
    201201  !       time en dehors des limites du fichier forcage 

  • trunk/SOURCES/climat_forcage_mois_mod.f90

    r27 r65  
    6464 
    6565  implicit none 
    66   integer :: i,j,mo,ti,tj   
    67   character(len=8) :: control   ! label to check clim. forc. file (filin) is usable 
    68   integer :: l                  ! In snapshot files:the first column is the mask, read but not used 
    69   real :: T_surf_ref            ! variable de travail calcul temp a l'instant t a la surface S 
    70   character(len=100) :: file_ncdf      !< fichier netcdf issue des fichiers .dat 
     66!  integer :: i,j,mo,ti,tj   
     67!  integer :: l                         ! In snapshot files:the first column is the mask, read but not used 
     68!  real :: T_surf_ref            ! variable de travail calcul temp a l'instant t a la surface S 
     69!  character(len=100) :: file_ncdf      !< fichier netcdf issue des fichiers .dat 
    7170  real*8, dimension(:,:,:), pointer :: data_3D => null() ! donnees lues dans le netcdf 
    7271 

  • trunk/SOURCES/conserv-mass-adv-diff_sept2009_mod.f90

    r4 r65  
    8383real,dimension(nx,ny) :: Dminx,Dminy 
    8484real,dimension(nx,ny) :: Uxdiff,Uydiff           ! vitesse due a la diffusion 
    85 integer :: it1,it2,it3 
    86  
    8785real aux                                         ! pour le calcul du critere 
    8886real maxdia                                      ! sur le pas de temps 

  • trunk/SOURCES/diffusiv-polyn-0.6.f90

    r59 r65  
    4444REAL :: INV_4DX ! inverse de dx pour eviter les divisions =1/(4*dx) 
    4545REAL :: INV_4DY ! inverse de dy pour eviter les divisions =1/(4*dy) 
    46 INTEGER :: ii 
    47 integer :: ll 
    4846 
    4947if (itracebug.eq.1)  call tracebug(' Entree dans routine diffusiv') 

  • trunk/SOURCES/initial-0.3.f90

    r59 r65  
    3030  use sorties_ncdf_grisli 
    3131  use util_recovery 
    32   character(len=80) :: filinit 
    3332 
    3433  !------------------------------------------------------------------------------------- 

  • trunk/SOURCES/litho-0.4.f90

    r4 r65  
    4545      implicit none 
    4646 
    47       INTEGER :: IP,JP,LPX,LPY,II,SOM1,SOM2 
     47      INTEGER :: II,SOM1,SOM2 
    4848      REAL, dimension(:,:), allocatable :: WLOC 
    4949 

  • trunk/SOURCES/main3D-0.4-40km.f90

    r24 r65  
    116116 
    117117 
    118  
    119   integer :: m 
    120  
    121118  ! good luck 
    122119 

  • trunk/SOURCES/neffect-0.4.f90

    r4 r65  
    2424      implicit none 
    2525 
    26 !!!!  real pwater,N1,N2,XFRAC,SLOPEM,uxmy,uymx 
    27  
    28       real XFRAC,SLOPEM,uxmy,uymx 
    29 !     real,dimension(nx,ny) :: pwater !!! global var 
    30       real pentem,toslid,topente,neffoto,fonction_N 
    31       real,dimension(nx,ny) :: table_interm 
     26!      real,dimension(nx,ny) :: table_interm 
    3227      real :: Nefmin=1.e5                  ! Pression effective minimum (~ 10 m de glace) 
    3328 

  • trunk/SOURCES/next-time.f90

    r4 r65  
    2121  real  (kind=kind(0.d0)) :: time 
    2222  real (kind=kind(0.d0)) :: time_loc 
    23   real (kind=kind(0.d0)) :: test_time 
    2423  real (kind=kind(0.d0)) :: time2_loc 
    2524  real (kind=kind(0.d0)) :: timesup 

  • trunk/SOURCES/no_spinup_mod.f90

    r4 r65  
    1414module no_spinup 
    1515 
     16use module3d_phy 
    1617 
    1718contains 
    … …  
    2223subroutine init_spinup 
    2324 
    24 use module3d_phy 
    2525implicit none 
    2626 

  • trunk/SOURCES/out_cptr_mod.f90

    r8 r65  
    387387  character(len=20),dimension(4)    :: dimnames3d      !< pour 3d troisieme dim est nz 
    388388  character(len=20),dimension(4)    :: dimnames3dT     !< pour 3d troisieme dim est nz+nzm 
    389   integer                           :: nrecs=1         !< compteur pour les enregistrements temps des variables   
    390   integer                           :: idef=0          !< pour savoir si la variable a ete definie ou non 
     389!  integer                           :: nrecs=1         !< compteur pour les enregistrements temps des variables   
     390!  integer                           :: idef=0          !< pour savoir si la variable a ete definie ou non 
    391391  integer                           :: isortie         !< pour le choix de type de fichier en sortie  
    392392 

  • trunk/SOURCES/out_horiz_mod.f90

    r10 r65  
    7676    integer :: i2 
    7777    integer :: i3 
    78     integer  :: i,j,k !< indices de travail 
     78    integer  :: k !< indices de travail 
    7979 
    8080    integer :: num_dat = 21 
    … …  
    194194integer   :: ifin 
    195195integer   :: npr 
    196 integer  :: i,j,k ! indices de travail 
     196integer  :: i ! indices de travail 
    197197 
    198198 
    … …  
    281281real :: maxcol 
    282282real :: coef 
    283  integer  :: i,j,k ! indices de travail 
    284283 
    285284npos=npos+1 

  • trunk/SOURCES/out_profile_mod.f90

    r4 r65  
    146146integer :: numtime 
    147147!integer :: nmax 
    148 integer :: n,l 
     148integer :: n 
    149149integer, dimension(:), allocatable :: nbr_ligne_prof 
    150150integer :: nbr_ligne_total 

  • trunk/SOURCES/relaxation_water_diffusion.f90

    r41 r65  
    1010!< 
    1111module relaxation_waterdif_mod 
     12 
     13 
    1214 
    1315 
    … …  
    4850    REAL :: RESTE,DELH,VH 
    4951    INTEGER :: ntour 
    50     INTEGER :: mbord 
    5152    REAL  :: DTSRGDX,dtwdx2 
    5253    LOGICAL :: STOPP 

  • trunk/SOURCES/resol_adv_diff_2D-sept2009.f90

    r4 r65  
    146146 
    147147 
    148 integer :: it1,it2,jt1,jt2           ! pour des tests d'asymétrie 
    149  
    150  
    151148if (itracebug.eq.1)  call tracebug(' Entree dans routine resolution_diffusion') 
    152149 

  • trunk/SOURCES/snaptime.f90

    r4 r65  
    2424double precision :: tunit 
    2525integer :: ipuissance 
    26 integer :: itime,m 
     26integer :: itime 
    2727character(len=1) :: signe,unit,mm  
    2828character(len=1),dimension (4) :: unittab 

  • trunk/SOURCES/spinup_mod.f90

    r29 r65  
    516516 
    517517    implicit none 
    518     real,dimension(nx,ny) :: corr_ghf        !< correction du gflux geothermique 
    519518    real,dimension(nx,ny) :: coefdef_maj     !< coefficient deformation 
    520519    real :: increment_ghf 

  • trunk/SOURCES/write_datfile.f90

    r4 r65  
    2222 
    2323integer  :: i,j                ! working integers 
    24 integer  :: lx,ly              ! nxx, nyy read in the file 
    25  
    2624 
    2725open(22,file=trim(filename)) 
    … …  
    5553 
    5654integer  :: i,j                ! working integers 
    57 integer  :: lx,ly              ! nxx, nyy read in the file 
    58  
    5955 
    6056open(22,file=trim(filename)) 
Note: See TracChangeset for help on using the changeset viewer.