Changeset 481 for branches


Ignore:
Timestamp:
02/09/24 10:58:55 (4 months ago)
Author:
aquiquet
Message:

Cleaning branch: unused variables removed following strict compilation

Location:
branches/GRISLIv3/SOURCES
Files:
34 edited

Legend:

Unmodified
Added
Removed

  • branches/GRISLIv3/SOURCES/3D-physique-gen_mod.f90

    r474 r481  
    3232  integer ::  i_resolmeca               !< defini le type d'association SIA-L1 
    3333  integer ::  inv_beta                  !< 0 : run std, 1 : inversion du beta 
    34   integer ::  ICOUPLE                   !<  
    3534  integer ::  ICOMPTEUR                 !< reprise fichier cptr  
    3635  integer ::  IOUT                      !< type de fichier de sortie pour reprise 
    … …  
    4039  integer ::  NTMAX                     !<  
    4140 
    42   integer ::  ITEMP                     !<  
    4341  integer ::  niter_nolin               !< number of velocity iterations for non-linear dragging laws 
    4442  integer ::  ndebug                    !< pour faire des sorties de debug a plusieurs pas de temps 

  • branches/GRISLIv3/SOURCES/Ant16_files/bmelt-seuil-profondeur_initmip_inv_mod.f90

    r446 r481  
    6969    real*8, dimension(:,:),   pointer    :: tab               !< tableau 2d real pointer 
    7070    character(len=100) :: file_number_shelves ! fichier avec les zones ice-shelves 
    71     character(len=100) :: file_bmelt_anom 
    7271    real,dimension(nx,ny) :: dist_talus   ! distance au talus 
    73     character(len=100) :: file_dist_talus ! pas utilise ici mais pour compatibilite lecture namelist 
     72    character(len=100) :: file_dist_talus  
    7473    real :: bmelt_dist0                   ! pas utilise ici mais pour compatibilite lecture namelist 
    7574     

  • branches/GRISLIv3/SOURCES/Climate_modules/climat-perturb_mod-0.4.f90

    r467 r481  
    193193 
    194194 
    195   character(len=80):: filin 
    196   integer :: i,j 
    197   integer :: ivo,jvo 
    198  
    199195namelist/clim_pert/coefT,rappact,retroac,rapbmshelf,mincoefbmelt,maxcoefbmelt,filforc 
    200196 
    … …  
    217213write(num_rep_42,*)'/'                             
    218214write(num_rep_42,*) 
    219  
    220  
    221  
    222 !cdc Commente pour être compatible avec lecture fichiers Cat Schoofing 
    223  
    224 !~ !!!!!!!! ATTENTION AJOUTE POUR TEST MAIS A REMETTRE AU PROPRE PLUS TARD C. DUMAS !!!!!!!!!! 
    225 !~ !!!!!!!! ancien input_climat_ref de lect_clim_act_anteis 
    226 !~ !     accumulation de Philippe 
    227 !~       filin='accumHUY40km.dat' 
    228 !~       call lect_eis(nx,ny,precip,filin,DIRNAMEINP) 
    229 !~ !====================================== La reponse est 42 =========== 
    230 !~       write(num_rep_42,*) 'fichier accum : ', filin  
    231  
    232 !~ !     cas particulier de Vostok 
    233 !~       ivo=101 
    234 !~       jvo=62 
    235 !~       do j=jvo-1,jvo+1 
    236 !~         do i=ivo-1,ivo+1 
    237 !~             precip(i,j)=0.02         ! valeur plus faible a Vostok. 
    238 !~         end do 
    239 !~       end do 
    240 !~       acc(:,:)=precip(:,:) 
    241  
    242 !~ !    temperature en surface : 
    243 !~ !    parametrisation de Fortuin pour la temperature annuelle. 
    244 !~       do j=1,ny 
    245 !~         do i=1,nx 
    246  
    247 !~             if (s0(i,j).le.200.) then   ! shelfs 
    248 !~                tann(i,j)=49.642-0.943*abs(ylat(i,j)) 
    249 !~             else if ((s0(i,j).gt.200.).and.(s0(i,j).lt.1500.)) then ! pente 
    250 !~                tann(i,j)=36.689-0.005102*s0(i,j)-0.725*abs(ylat(i,j)) 
    251 !~             else if (s0(i,j).ge.1500.) then        ! plateau 
    252 !~                tann(i,j)=7.405-0.014285*s0(i,j)-0.180*abs(ylat(i,j)) 
    253 !~             endif 
    254  
    255 !~             ta0(i,j)=tann(i,j) 
    256 !~ !           pour la temperature d'ete, idem parametrisation huybrechts 
    257 !~             tjuly(i,j)=tann(i,j)-17.65+0.00222*s0(i,j)& 
    258 !~                          +0.40802*abs(ylat(i,j)) 
    259 !~         end do 
    260 !~       end do 
    261 !~ !!!!!!!! FIN MODIF  TEMPORAIRE !!!!!!!!!! 
    262 !cdc fin Commente pour être compatible avec lecture fichiers Cat Schoofing 
    263215 
    264216return 

  • branches/GRISLIv3/SOURCES/Climate_modules/climat_InitMIP_years_perturb_mod.f90

    r451 r481  
    8282 
    8383  integer            :: err                ! recuperation d'erreur 
    84   integer            :: i,k 
    85  
    86   !aurel for Tafor: 
    87   character(len=8)   :: control      !label to check clim. forc. file (filin) is usable 
    88   character(len=80)  :: filin 
    89    
     84  integer            :: k 
     85  
    9086  real               :: time_depart_snaps  !> temps du debut premier snapshot 
    9187 
    9288  namelist/clim_smb_T_gen/smb_file,coef_smb_unit,temp_annual_file 
    93 ! version Aurelien Groenland avec file_lapse : 
    94 !  namelist/smb_anom_initMIP/file_smb_anom,file_lapse,coef_smb_anom_unit,nb_snap,time_depart_snaps,massb_time 
    95 ! version utilisee pour Antarctique :   
    9689  namelist/smb_anom_initMIP/file_smb_anom,file_lapse,coef_smb_anom_unit,nb_snap,time_depart_snaps,massb_time 
    9790   
    98 !!!!!  namelist/clim_snap/nb_snap,file_smb_snap,massb_time_snap 
    9991 
    10092428 format(A) 
    … …  
    319311  implicit none 
    320312 
    321   integer i,ift 
    322313  real             :: coefanomtime 
    323314   

  • branches/GRISLIv3/SOURCES/Climate_modules/climat_forcage_mod.f90

    r451 r481  
    8080  implicit none 
    8181  character(len=8) :: control     !< label to check clim. forc. file (filin) is usable 
    82   integer :: l                    !< in snapshot files:the first column is the mask, read but not used 
    8382  integer :: err                  !< recuperation erreur 
    84   integer :: i,j,k,m 
     83  integer :: i,k,m 
    8584  real*8, dimension(:,:), pointer   :: tab2d => null()   !< tableau 2d pointer 
    8685  real*8, dimension(:,:,:), pointer :: tab3d => null()   !< tableau 3d pointer 
    … …  
    292291   ttr(k)=ttr_temp(k) 
    293292  enddo 
    294     
    295  
    296 !~   write(num_rep_42,'(A)')'!___________________________________________________________'  
    297 !~   write(num_rep_42,'(A)') '&clim_forcage                                ! nom du bloc ' 
    298 !~   write(num_rep_42,*) 
    299 !~   write(num_rep_42,'(A,A,A)') 'clim_ref_file = "',trim(clim_ref_file),'"' 
    300 !~   write(num_rep_42,*) 'ttr = ', ttr 
    301 !~ !  write(num_rep_42,'(A,A,A)') 'forcage_file= "',trim(forcage_file),'"' 
    302 !~   write(num_rep_42,'(A,A,A)') 'forcage_file1 = "',trim(forcage_file1),'"' 
    303 !~   write(num_rep_42,'(A,A,A)') 'forcage_file2 = "',trim(forcage_file2),'"' 
    304 !~   write(num_rep_42,*) 'typerun = ', typerun 
    305 !~   write(num_rep_42,*) 'lapserate = ', lapserate 
    306 !~   write(num_rep_42,*) 'rappact = ', rappact 
    307 !~   write(num_rep_42,'(A,A,A)') ' filforc      = "',trim(filforc),'"' 
    308 !~   write(num_rep_42,*) 'pertbmb = ', pertbmb 
    309 !~   write(num_rep_42,*) 'coefT   = ', coefT 
    310 !~   write(num_rep_42,*) 'coefbmb = ', coefbmb 
    311 !~   write(num_rep_42,*) 'r_atmvar = ', r_atmvar 
    312 !~   write(num_rep_42,*)'/'                             
    313 !~   write(num_rep_42,*) 
    314293 
    315294  PYY=2.*PI/50. 
    … …  
    330309  integer :: l                       !< dumm index for loops on snapsots files  l=itr,ntr-1 
    331310  integer :: itr                     !< nb of the current snapshot file (change with time) 
    332   integer :: i,j,k,m 
     311  integer :: i,m 
    333312  integer :: ift                     !< indice correspondant au pas de temps du fichier de forcage 
    334   real :: TEMP                       !< calcul du nbr de jour < psolid 
    335313  real,dimension(nx,ny) :: deltaZs   ! diff temp par rapport a topo orog0  
    336314  real,dimension(nx,ny,12) :: prmois ! precip sur topo Zs 

  • branches/GRISLIv3/SOURCES/Climate_modules/climat_smb_perturb_mod.f90

    r452 r481  
    6565  character(len=100) :: smb_file           ! fichier smb 
    6666  character(len=100) :: temp_annual_file   ! temperature annuelles 
    67   character(len=100) :: file_smb_snap      !> nom du fichier dans lequel sont les snapshots 
    6867 
    6968  integer            :: err                ! recuperation d'erreur 
    70   integer            :: i,j,k 
     69  integer            :: i 
    7170 
    7271  !aurel for Tafor: 

  • branches/GRISLIv3/SOURCES/Climate_modules/climat_transient_GCM_mod.f90

    r451 r481  
    4747  real*8, dimension(:,:), pointer        :: tab2D             !< tableau 2d real pointer 
    4848  real*8, dimension(:,:,:), pointer      :: tab3D             !< tableau 3d real pointer 
    49   real*8, dimension(:,:,:,:), pointer    :: tab4D             !< tableau 4d real pointer 
    5049   
    5150  real               :: time_depart_snaps  !> temps du debut premier snapshot 

  • branches/GRISLIv3/SOURCES/Draggings_modules/dragging_beta_iter_vitbil_mod.f90

    r449 r481  
    2323 
    2424    use io_netcdf_grisli, only: read_ncdf_var 
    25     use runparam, only: itracebug,num_tracebug 
     25    use runparam, only: itracebug 
    2626    use geography, only: dirnameinp,nx,ny 
    2727    use module3D_phy, only: num_param,num_rep_42,inv_beta,betamax,betamax_2d,mk_init,flot,beta_centre, & 

  • branches/GRISLIv3/SOURCES/Draggings_modules/dragging_coulomb_friction_mod.f90

    r473 r481  
    185185 
    186186use geography,    only: nx,ny 
    187 use module3d_phy, only: betamax,beta_centre,betamx,betamy,neffmx,neffmy,hwater,flot,ux,uy, & 
     187use module3d_phy, only: betamax,beta_centre,betamx,betamy,neffmx,neffmy,hwater,flot, & 
    188188                        uxbar,uybar 
    189189 

  • branches/GRISLIv3/SOURCES/Draggings_modules/dragging_coulomb_friction_simplhydro_mod.f90

    r408 r481  
    199199 
    200200use geography,    only: nx,ny 
    201 use module3d_phy, only: betamax,beta_centre,betamx,betamy,hwater,flot,ux,uy,H,B,bsoc,sealevel_2d, & 
     201use module3d_phy, only: betamax,beta_centre,betamx,betamy,hwater,flot,H,bsoc,sealevel_2d, & 
    202202                        uxbar,uybar 
    203203use param_phy_mod,only: G,RO,ROFRESH 

  • branches/GRISLIv3/SOURCES/Draggings_modules/dragging_prescr_beta_mod.f90

    r408 r481  
    5252                                ! nouvelle version : lit les fichiers nc 
    5353 
    54     use runparam, only: itracebug,num_tracebug 
     54    use runparam, only: itracebug 
    5555    use geography, only: dirnameinp 
    5656    use module3d_phy, only: num_param,num_rep_42,inv_beta,betamax,betamax_2d,mk_init,h,bsoc, & 

  • branches/GRISLIv3/SOURCES/Fichiers-parametres/Makefile.tof-lsce3130.inc

    r361 r481  
    2323 
    2424# mettre debug à 1 pour compiler avec options strictes 
    25 debug ?= 0 
     25debug ?= 1 
    2626 
    2727# compilation avec ifort : 

  • branches/GRISLIv3/SOURCES/Netcdf-routines/io_netcdf_GRISLI.f90

    r449 r481  
    839839  subroutine Write_Ncdf_var2d_Real(varname,dimname,file,tabvar,typevar) 
    840840    use netcdf,only: nf90_open,nf90_close,nf90_write,nf90_noerr,nf90_inq_varid, & 
    841                      nf90_def_var,nf90_put_var,nf90_inq_dimid,nf90_enddef,nf90_redef 
     841                     nf90_def_var,nf90_put_var,nf90_inq_dimid,nf90_enddef,nf90_redef, & 
     842                     nf90_double,nf90_float 
    842843    !       
    843844    !      implicit none 
    … …  
    14271428  !!@param[in]     typevar : type of the variable to write 
    14281429  !!@return idef 
    1429   subroutine Write_Ncdf_var1d_Int_t(varname,dimname,file,tabvar,time,idef,typevar) 
     1430  subroutine Write_Ncdf_var1d_Int_t(varname,dimname,file,tabvar,time,idef) 
    14301431    use netcdf,only: nf90_open,nf90_close,nf90_write,nf90_noerr,nf90_inq_varid,nf90_double, & 
    14311432                     nf90_float,nf90_def_var,nf90_put_var,nf90_inq_dimid,nf90_enddef,nf90_redef 
    … …  
    14331434    implicit none 
    14341435    !        
    1435     Character(*),intent(in) :: varname,file,typevar 
     1436    Character(*),intent(in) :: varname,file 
    14361437    Character(*),dimension(2),intent(in) :: dimname 
    14371438    Integer :: time 

  • branches/GRISLIv3/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90

    r471 r481  
    665665    use geography,    only: nx,ny,nz,nzm,dx,dy 
    666666    use runparam,     only: itracebug,xmin,ymin,dirsource 
    667     use module3d_phy, only: s,s0,h,h0,b,sealevel,mk,mk_init,hdot,dhdt,tann,tjuly,t,tpmp, & 
    668                             acc,abl,ghf,phid,bmelt,uxbar,uybar,ux,uy,frotmx,frotmy,tobmx,tobmy,    & 
     667    use module3d_phy, only: s,s0,h,b,sealevel,mk,mk_init,hdot,dhdt,tann,tjuly,t,tpmp, & 
     668                            acc,ghf,phid,bmelt,uxbar,uybar,ux,uy,frotmx,frotmy,tobmx,tobmy,    & 
    669669                            taushelf,epsxx,epsyy,epsxy,eps,abar,pvi,pvm,betamx,betamy,beta_centre, & 
    670670                            ablbord_dtt,ice,front,gr_line_schoof,hwater,hdotwater,bsoc,bm,    & 

  • branches/GRISLIv3/SOURCES/New-remplimat/diagno-L2_mod.f90

    r471 r481  
    683683subroutine mismip_boundary_cond 
    684684 
    685 use module3d_phy, only: imx_diag,imy_diag,uxbar,uybar 
     685use module3d_phy, only: imx_diag,imy_diag,uybar 
    686686use runparam, only: itracebug 
    687687 

  • branches/GRISLIv3/SOURCES/New-remplimat/remplimat-shelves-tabTu.f90

    r478 r481  
    107107real, dimension(nx,ny) :: Hmy_oppos      !  
    108108 
    109 integer i,j,err 
     109integer i,j 
    110110 
    111111if (itracebug.eq.1)  call tracebug(' Subroutine rempli_L2') 

  • branches/GRISLIv3/SOURCES/Temperature-routines/icetemp_mod.f90

    r467 r481  
    1717  !> 
    1818 
    19   Subroutine init_icetemp(Num_rep_42_m) 
     19  Subroutine init_icetemp 
     20 
    2021    Use Icetemp_declar, only: Nfracq,Iq,Ecart_phid 
    2122    Use module3D_phy,   only: tbdot 
    … …  
    2324    Implicit None 
    2425 
    25     !< Arguments 
    26     Integer,                       Intent(In)::Num_rep_42_m !< Id Of Reponse 
    27     ! Character(Len=80) :: Filin 
    2826    !________________________________________________________________ 
    2927    ! La Temperature Du Point De Congelation De L'Eau De Mer 
    … …  
    3634    Nfracq=1 
    3735 
    38 !    Write(Num_rep_42_m,*) 'Calcul Des Temperatures' 
    39 !    Write(Num_rep_42_m,*) 'Nfracq=',Nfracq 
    40  
    4136    Iq=1 
    42 !    If (Iq.Eq.1) Write(Num_rep_42_m,*) 'Chaleur Demi Maille, Iq=',Iq 
    43 !    If (Iq.Eq.2) Write(Num_rep_42_m,*) 'Chaleur Centree, Iq=',Iq 
    44 !    If (Iq.Eq.3) Write(Num_rep_42_m,*) 'Chaleur Pente, Iq=',Iq 
    45 !    If (Iq.Eq.4) Write(Num_rep_42_m,*) 'Chaleur U_taub, Iq=',Iq 
    46 !    If (Iq.Eq.5) Write(Num_rep_42_m,*) 'Chaleur U_taub_stag, Iq=',Iq 
    47 !    If (Iq.Eq.6) Write(Num_rep_42_m,*) 'Chaleur Q_sia_stag,Iq=',Iq 
    48 !    If (Iq.Eq.7) Write(Num_rep_42_m,*) 'Chaleur Q_all_stag,Iq=',Iq 
    49 !    Write(Num_rep_42_m,*) 
    5037 
    5138    Ecart_phid=0.5 
    52 !    Write(Num_rep_42_m,*) ' Si Base Froide, La Chaleur De Glissement Est En Exp(Delta T * Ecart_phid)' 
    53 !    Write(Num_rep_42_m,*) 'Ecart_phid=',Ecart_phid 
    54 !    Write(Num_rep_42_m,*) '---------------------------------------------------------------------------' 
    5539 
    5640    tbdot(:,:) = 0.  
    … …  
    7054  subroutine icetemp 
    7155    !$ use OMP_LIB  
    72     use module3D_phy, only: uxbar,uybar,dtt,nzm,time,ux,uy,uzr,ts,h,b, & 
     56    use module3D_phy, only: uxbar,uybar,dtt,nzm,time,ux,uy,uzr,ts,h, & 
    7357                             ghf,t,tpmp,phid,flot,ibase,debug_3D,tbdot 
    7458    use runparam, only: itracebug,num_tracebug 

  • branches/GRISLIv3/SOURCES/bmelt-beckmann-gcm_mod.f90

    r472 r481  
    167167    call Read_Ncdf_var('thetao',file_inputs,tab3d) 
    168168    temp_ocean(:,:,:)  = tab3d(:,:,:)  
    169          
     169 
    170170    file_inputs=TRIM(DIRNAMEINP)//'Snapshots-GCM/'//TRIM(SO_file) 
    171171    call Read_Ncdf_var('so',file_inputs,tab3d) 
    … …  
    339339!   print*, 'salinity_ocean', salinity_ocean(1,1,:) 
    340340 
    341         !temp_ocean(:,:,:) = temp_ocean(:,:,:) - 273.15 
    342          
     341    !temp_ocean(:,:,:) = temp_ocean(:,:,:) - 273.15 
     342 
    343343    !coef_OM = K_t * 0.01420418516 
    344344    coef_OM = K_t * (row * cpw / (cl * ro))**2 
    … …  
    512512    ! This will be the netCDF ID for the file and data variable. 
    513513    integer :: ncid, varid, status 
    514     integer, dimension(3) :: count 
    515     integer, dimension(3) :: start 
    516514 
    517515    ! Open the file.  

  • branches/GRISLIv3/SOURCES/celltest_tracer.f90

    r446 r481  
    55 
    66      !!  Find the cell origin index for the previous time_tra step 
    7       use module3d_phy, only: e,uxbar,uybar,flotmx,flotmy,ux,uy,gzmx,gzmy 
     7      use module3d_phy, only: e,uybar,flotmy,uy,gzmy 
    88      use geography, only: nx,ny,nz,dx,dy 
    99      use tracer_vars, only: xgrid,ygrid 
    … …  
    8787            print *,"CellTest(y+;",time_tra,"): Out-of-range increment (",ic,",",jc,",",kc,")", v_y 
    8888!           print*, ytracer, ygrid(jm) 
    89             print*, uybar(ic,jc), uybar(ic,jc+1), flotmy(ic,jc), flotmy(ic, jc+1) 
    90             print *, uy(ic,jc,:) 
    91             print *, uy(ic,jc+1,:) 
    92             print*, gzmy(ic, jc), gzmy(ic,jc+1) 
     89            print*, uybar(ic,jc), uybar(ic,jc+1), flotmy(ic,jc), flotmy(ic, jc+1) 
     90            print *, uy(ic,jc,:) 
     91            print *, uy(ic,jc+1,:) 
     92            print*, gzmy(ic, jc), gzmy(ic,jc+1) 
    9393!           STOP 
    9494             jm = jc-1 
    … …  
    114114  if(E_k < E(1)) then 
    115115        E_k = E(1) 
    116         km=1      
     116        km=1 
    117117    
    118118  else if(E_k > E(nz)) then 
    119119        E_k = E(nz) 
    120         km=nz-1 
     120        km=nz-1 
    121121 
    122122  else   ! this used to be a loop on its own, now shortcuts if E_k is special  
    123123 
    124124      if(v_E>=0.0) then   !! Downward flow of ice in E system 
    125         if (kc==1) then 
    126           !! No action, pending tactical discussion with GKCC 
     125        if (kc==1) then 
     126          !! No action, pending tactical discussion with GKCC 
    127127          km = 1 
    128128        else 
    129           km=kc - 1 
     129          km=kc - 1 
    130130          do 
    131131            if( (E_k >= E(km)) .and. (E_k <= E(km+1)) ) exit 

  • branches/GRISLIv3/SOURCES/flottab2-0.7.f90

    r477 r481  
    6868    use module3D_phy, only:igrdline,mk_init,flot,H,sealevel_2d,Bsoc,S,H,B,& 
    6969         ice,front,& 
    70          iceberg,uxbar,uybar,mk,gzmx,gzmy,flotmx,flotmy,hmx,hmy,isynchro,& 
     70         iceberg,uxbar,uybar,mk,gzmx,gzmy,flotmx,flotmy,hmx,hmy,& 
    7171         flgzmx,flgzmy,fbm,bm,bmelt,debug_3D,dt 
    7272    use param_phy_mod, only:row,ro 
    … …  
    8383    real,dimension(nx,ny) ::  uys1   !< uybar a l'entree de flottab 
    8484     
    85     integer ::  nb_pt 
    8685    real    ::  petit_H=0.001 ! pour test ice sur zone flottante 
    8786     
    88     integer :: i,j,ii,jj 
     87    integer :: i,j 
    8988     
    9089    if (itracebug.eq.1)  call tracebug(' Entree dans routine flottab') 

  • branches/GRISLIv3/SOURCES/furst_schoof_mod.f90

    r470 r481  
    130130                      frot_coef = betamx(i+1,j) 
    131131                   endif 
    132                    call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),m_weert,bfx,phi_prescr) 
     132                   call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),bfx,phi_prescr) 
    133133                else 
    134134                   print*,'ATTENTION FLUX AUTRE QUE TSAI OU SCHOOF NON IMPLEMENTE' 
    … …  
    220220                      frot_coef = betamx(i+1,j) 
    221221                   endif 
    222                    call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),m_weert,bfx,phi_prescr) 
     222                   call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),bfx,phi_prescr) 
    223223                else 
    224224                   print*,'ATTENTION FLUX AUTRE QUE TSAI OU SCHOOF NON IMPLEMENTE' 
    … …  
    307307                      frot_coef = betamy(i,j+1) 
    308308                   endif 
    309                    call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),m_weert,bfy,phi_prescr) 
     309                   call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),bfy,phi_prescr) 
    310310                else 
    311311                   print*,'ATTENTION FLUX AUTRE QUE TSAI OU SCHOOF NON IMPLEMENTE' 
    … …  
    391391                      frot_coef = betamy(i,j+1) 
    392392                   endif 
    393                    call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),m_weert,bfy,phi_prescr) 
     393                   call flux_Schoof4Schoof (Hgl,Abar(i,j),frot_coef,alpha_flot,glen(1),bfy,phi_prescr) 
    394394                else 
    395395                   print*,'ATTENTION FLUX AUTRE QUE TSAI OU SCHOOF NON IMPLEMENTE' 
    … …  
    573573  ! Dans un premier temps, pas de terme horizontal (regarder dans toy_recul une tentative de confinement) 
    574574 
    575   subroutine flux_Schoof4Schoof (Hgl,A_mean,C_frot,alpha,n_glen,m_weert,back_force_coef,phi_schoof) 
     575  subroutine flux_Schoof4Schoof (Hgl,A_mean,C_frot,alpha,n_glen,back_force_coef,phi_schoof) 
    576576 
    577577    use param_phy_mod, only:rog 
    … …  
    584584    real,intent(in) :: alpha 
    585585    real,intent(in) :: n_glen 
    586     real,intent(in) :: m_weert 
    587586    real,intent(in) :: back_force_coef 
    588587    real,intent(out) :: phi_schoof 

  • branches/GRISLIv3/SOURCES/initial-0.3.f90

    r474 r481  
    2222subroutine initial 
    2323 
    24   use module3D_phy, only:time,H,hmx,hmy,T,tpmp,bmelt,icompteur,num_rep_42 
     24  use module3D_phy, only:time,H,hmx,hmy,T,tpmp,bmelt,icompteur 
    2525  use geography, only: geoplace 
    2626  use runparam, only: tbegin 
    … …  
    8585 
    8686  if (geoplace(1:5).ne.'mism3') then 
    87      call init_icetemp(num_rep_42)        ! initialise le calcul de temperature 
     87     call init_icetemp                    ! initialise le calcul de temperature 
    8888     call thermal_prop_icetemp            ! appel aux proprietes thermiques 
    8989     call init_bmelt                      ! Appel aux tableaux d'initialisation des ice shelves 

  • branches/GRISLIv3/SOURCES/interpolate_tracer.f90

    r409 r481  
    99  integer, intent(IN) :: ic, jc, kc, im, jm, km 
    1010  real, intent(INOUT) :: fx, fy, fz, E_k 
    11   real :: fx_1, fy_1, fz_1, a_xy, a_xy_p, f_zt, a_tmp, ac , al 
     11  real :: fx_1, fy_1, fz_1, a_xy, a_xy_p, a_tmp, ac , al 
    1212  real :: fra1, fra2, fra0, accucumul1, accucumul2, accucumul0, & 
    1313       thinning, thinning2, f_k, a_xy_m, fzz, fprime, fsec 

  • branches/GRISLIv3/SOURCES/isostasie_mod-0.3.f90

    r468 r481  
    3333 
    3434    use iso_declar,only: nlith,dt_iso,tausoc,dl,rl,lbloc,we,charge 
    35     use module3D_phy, only: icouple,err,h0,bsoc0,sealevel_2d 
     35    use module3D_phy, only: err,h0,bsoc0,sealevel_2d 
    3636    use geography, only: geoplace,nx,ny,dx,dy 
    3737    use param_phy_mod, only: ro,row,rog,rowg,romg 

  • branches/GRISLIv3/SOURCES/lake_rsl_mod.f90

    r449 r481  
    3535  implicit none 
    3636 
    37   !~ !* General & discrete parameters 
    38   !~ real :: undef=-9999. !undef=1.e20 
    39  
    40  
    41   !~ !* Parameters and arrays used in the Initialization step 
    42   real, dimension(nx,ny) :: orog, orog_lake, lake_level 
    43   integer, dimension(nx,ny) :: omsk, cmsk 
    44   integer, dimension(nx,ny) :: cont!, conttrack, contsave, runoff_msk, mask_search, exut 
    45   integer, dimension(nx,ny) :: ocean 
    46   integer, dimension(nx,ny) :: pot_lake 
    47   integer :: hi_res_runoff ! select 1 : interpolation on hi resolution topo or 0 : no interpolation 
    48  
    49   integer :: countcont ! nbr de continents 
    50   !~ integer :: nbocean, countocean 
    51  
    52   !~ !* Parameters and arrays of the runoff calculation 
    53   integer, dimension(nx,ny) :: rtm 
    54  
    55   !~ !* Debug mode & nc files 
    56   character(len=100) :: file_topo_hi_res, coord_topo_hi_res 
    57   !~ character(len=15) :: invar 
    58   !~ integer :: ncid, xdimid, ydimid, varid 
    59   !~ integer :: latid, lonid, rtmid, orogid, cmskid, contid, oroglakeid, lake_levelid 
     37  integer :: hi_res_runoff ! select 1 : interpolation on hi resolution topo or 0: no interpolation 
    6038 
    6139  !* Routing code parameter 
    6240  integer, dimension(8,2) :: inc 
    63  
    64   !~ !cdc variable ajoutees 
    65  
    66   !~ integer, dimension(8) :: loccont ! cont des 8 points autour du point ii jj 
    67  
    68   !~ integer :: sumcont ! somme des cont des 8 points autour du point ii jj 
    69  
    70  
    71   ! 
    72   !~ integer :: maxloop ! pour debug evite boucle infinies 
    73   !~ integer,dimension(2) :: ijcoord ! coordonnees ij du point le plus bas 
    74   !~ integer,dimension(2) :: ijcoord_cmsk ! coordonnees ij du point cotier le plus bas 
    75   !~ real, dimension(8) :: locorog ! orog des 8 points voisins 
    76   integer,dimension(nx,ny) :: lake ! 1 si lac, 0 sinon 
    77   !~ integer,dimension(8) :: locexut ! 1 si exutoire => permet de gerer les lacs 
    78  
    79   !~ !integer :: contprint=0 ! pour debug et arret du print 
    80   !~ integer :: contexut ! compte les points sous le niveau du lac 
    81  
    82   !~ logical :: stop_lake 
    83   !~ integer, dimension(2) :: ij_lake 
    84   !~ integer,dimension(8) :: locomsk ! omsk des 8 voisins 
    85   !~ integer :: kloc 
    86  
    8741  integer,dimension(80,2) :: inclake ! 80 voisins (+4 -4) en i,j 
    88   integer, dimension(nx,ny) :: mask_extlake ! extension du lac sous la glace 
    89   real, dimension(nx,ny) :: orog_extlake ! niveau du lac sur les zones glace a proximité (rayon de 4 pts de grille) 
    90   !~ integer, dimension(nx,ny) :: clake ! point lac avec voisin point englace 
    9142  real, dimension(nx,ny) :: sealevel_2D_loc 
    9243 
    … …  
    11263  real, dimension(nxhi,nyhi) :: anom_Bsoc_hi, anom_S_hi, anom_H_hi ! anomalie de topo resol std vs resol_hi 
    11364  real, dimension(nxhi,nyhi) :: Bsoc_hi, S_hi, H_hi ! topo haute resolution a l'instant t 
    114   real, dimension(nxhi,nyhi) :: orog_hi ! orographie hi_res 
    115   real, dimension(nxhi,nyhi) :: S_interp_hi, H_interp_hi ! interpolation de S et H low res vers hi_res sans anomalie (pour faire topo englacee non bruitee) 
    116   integer, dimension(nxhi,nyhi) :: omsk_hi, cmsk_hi 
    117   integer, dimension(nxhi,nyhi) :: cont_hi ! masque des points continent hi_res 
    118   integer, dimension(nxhi,nyhi) :: ocean_hi ! masque des points ocean hi_res 
    119   integer, dimension(nxhi,nyhi) :: pot_lake_hi 
    120   integer, dimension(nxhi,nyhi) :: lake_hi ! 1 si lac, 0 sinon 
    121   integer, dimension(nxhi,nyhi) :: rtm_hi ! runoff direction (1 to 8) 
    122   real, dimension(nxhi,nyhi) :: lake_level_hi ! niveau du lac sur grille hi_res 
    123   real, dimension(nxhi,nyhi) :: orog_lake_hi  ! orographie a la surface des lacs 
    124  
    125   integer :: countcont_hi ! nbr de continents sur grille hi_res 
     65 
    12666 
    12767contains 
    … …  
    14080 
    14181    namelist/lake_rsl/hi_res_runoff, file_topo_hi_res, coord_topo_hi_res 
     82 
     83    !~ !* Debug mode & nc files 
     84    character(len=100) :: file_topo_hi_res, coord_topo_hi_res 
     85 
    14286 
    14387    integer :: ii, i, j, k 
    … …  
    158102    write(num_rep_42,*) 
    159103 
    160     !~   print*,'Debut input_rsl' 
    161  
    162     !* Continent number (outfile) 
    163     !ncfilecont = "continent_number_test-0.05.nc" 
    164     !ncfilecmsk = "coastal_mask_test-0.05.nc" 
    165     !* Runoff & corrected topo 
    166     !ncfilernff = "runoff_20Ma_Fred.nc" 
    167     !ncfiletopoout = "Topo.20Ma_Fred.OK_runoff.nc" 
    168     !ncfilernff = "runoff_NorthAm_15ka_test-0.05.nc" 
    169     !ncfiletopoout = "Topo.NorthAm_15ka_test.OK_runoff-0.05.nc" 
    170104 
    171105    !*************************** 
    … …  
    237171       enddo 
    238172       close(2004) 
    239  
    240        !  xmin=xcc(1,1)/1000. 
    241        !  ymin=ycc(1,1)/1000. 
    242        !  xmax=xcc(nx,ny)/1000. 
    243        !  ymax=ycc(nx,ny)/1000. 
    244173 
    245174 
    … …  
    302231    implicit none 
    303232 
    304     integer :: i, j, k, l 
    305     integer :: ii, jj 
    306     ! pour les lectures ncdf 
    307     real*8, dimension(:,:),   pointer    :: tab => null()       !< tableau 2d real pointer 
     233    integer :: i, j, k 
     234 
     235    !~ !* Parameters and arrays used in the Initialization step 
     236    real, dimension(nx,ny) :: orog, orog_lake, lake_level 
     237    integer, dimension(nx,ny) :: omsk, cmsk 
     238    integer, dimension(nx,ny) :: cont!, conttrack, contsave, runoff_msk,mask_search, exut 
     239    integer, dimension(nx,ny) :: ocean 
     240    integer, dimension(nx,ny) :: pot_lake 
     241 
     242    integer :: countcont ! nbr de continents 
     243    !~ integer :: nbocean, countocean 
     244 
     245    !~ !* Parameters and arrays of the runoff calculation 
     246    integer, dimension(nx,ny) :: rtm 
     247 
     248    integer,dimension(nx,ny) :: lake ! 1 si lac, 0 sinon 
     249    integer, dimension(nx,ny) :: mask_extlake ! extension du lac sous la glace 
     250    real, dimension(nx,ny) :: orog_extlake ! niveau du lac sur les zones glace a proximité (rayon de 4 pts de grille) 
    308251 
    309252    real, dimension(nx,ny) :: orog_lake_tmp ! pour calculer orog_lake a partir de orog_lake_hi 
     253 
     254    real, dimension(nxhi,nyhi) :: orog_hi ! orographie hi_res 
     255    real, dimension(nxhi,nyhi) :: S_interp_hi, H_interp_hi ! interpolation de S et H low res vers hi_res sans anomalie (pour faire topo englacee non bruitee) 
     256    integer, dimension(nxhi,nyhi) :: omsk_hi, cmsk_hi 
     257    integer, dimension(nxhi,nyhi) :: cont_hi ! masque des points continent hi_res 
     258    integer, dimension(nxhi,nyhi) :: ocean_hi ! masque des points ocean hi_res 
     259    integer, dimension(nxhi,nyhi) :: pot_lake_hi 
     260    integer, dimension(nxhi,nyhi) :: lake_hi ! 1 si lac, 0 sinon 
     261    integer, dimension(nxhi,nyhi) :: rtm_hi ! runoff direction (1 to 8) 
     262    real, dimension(nxhi,nyhi) :: lake_level_hi ! niveau du lac sur grille hi_res 
     263    real, dimension(nxhi,nyhi) :: orog_lake_hi  ! orographie a la surface des lacs 
     264 
     265    integer :: countcont_hi ! nbr de continents sur grille hi_res 
     266 
    310267    ! pour debug : ecriture resultats hi_res 
    311     character(len=100) :: ncfilecont, ncfilernff, ncfiletopoout,ncfilecmsk    ! nom fichier netcdf de sortie 
     268    character(len=100) :: ncfilecont    ! nom fichier netcdf de sortie 
    312269    character(len=1),dimension(2) :: dimname ! nom des dimensions du fichier netcdf de sortie 
    313270    integer :: ncid, xdimid, ydimid, S_hivarid, H_hivarid, Bsoc_hivarid, rtm_hivarid, cont_hivarid, ocean_hivarid, cmsk_hivarid, orog_lake_hivarid, lake_hivarid 
    … …  
    398355    ! mise à jour du niveau des lacs et du routage : 
    399356    if (mod(abs(TIME),500.).lt.dtmin) then 
    400        call mask_orog_ocean(nx,ny,orog,omsk,cmsk,cont,ocean,countcont,pot_lake) 
     357       call mask_orog_ocean(nx,ny,omsk,cmsk,cont,ocean,countcont,pot_lake) 
    401358 
    402359       !~   print*,'3/ rsl' 
    403360       if (hi_res_runoff.eq.1) then 
    404           call mask_orog_ocean(nxhi,nyhi,orog_hi,omsk_hi,cmsk_hi,cont_hi,ocean_hi,countcont_hi,pot_lake_hi) 
     361          call mask_orog_ocean(nxhi,nyhi,omsk_hi,cmsk_hi,cont_hi,ocean_hi,countcont_hi,pot_lake_hi) 
    405362          call runoff(nxhi,nyhi,countcont_hi,cont_hi,cmsk_hi,omsk_hi,pot_lake_hi,lake_hi,rtm_hi,lake_level_hi,orog_lake_hi) 
    406363          ! interpolation des variables sur grille std : 
    … …  
    444401       endwhere 
    445402 
    446        if (hi_res_runoff.eq.1) then 
     403       if (hi_res_runoff.eq.1) then   !afq-- this block is only for debug as it is called each timestep!! 
    447404 
    448405          !* Save mode 
    … …  
    494451  ! Subroutine de recherche des points terre, cotiers et ocean 
    495452  !---------------------------------------------------------------------------- 
    496   subroutine mask_orog_ocean(nx_loc,ny_loc,orog_loc,omsk_loc,cmsk_loc,cont_loc,ocean_loc,countcont_loc,pot_lake_loc) 
     453  subroutine mask_orog_ocean(nx_loc,ny_loc,omsk_loc,cmsk_loc,cont_loc,ocean_loc,countcont_loc,pot_lake_loc) 
    497454 
    498455    implicit none 
    499456 
    500457    integer,intent(in) :: nx_loc,ny_loc 
    501     real,dimension(nx_loc,ny_loc), intent(in) :: orog_loc 
    502458    integer,dimension(nx_loc,ny_loc), intent(inout) :: omsk_loc 
    503459    integer,dimension(nx_loc,ny_loc), intent(out) :: cmsk_loc, cont_loc 
    … …  
    744700    real,dimension(nx_loc,ny_loc), intent(out) :: lake_level_loc ! lake level 
    745701    real,dimension(nx_loc,ny_loc), intent(inout) :: orog_lake_loc  ! orographie a la surface des lacs 
    746  
    747     real,dimension(nx_loc,ny_loc) :: orog_extlake ! niveau du lac sur les zones glace a proximité (rayon de 4 pts de grille) 
    748702 
    749703    integer,dimension(nx_loc,ny_loc) :: exut ! > 0 si exutoire, 0 sinon 
    … …  
    917871 
    918872 
    919  

  • branches/GRISLIv3/SOURCES/lineartemp-0.2.f90

    r469 r481  
    1313 
    1414  use geography, only:nx,ny,nz,nzm 
    15   use module3D_phy, only:H,B,bmelt,T,TS,tpmp,ghf 
     15  use module3D_phy, only:H,bmelt,T,TS,tpmp,ghf 
    1616  use icetemp_declar, only:Dzm,Cm 
    1717   

  • branches/GRISLIv3/SOURCES/main3D-0.4-40km.f90

    r476 r481  
    131131subroutine grisli_init 
    132132 
    133   USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen, & 
     133  USE module3D_phy, only: isynchro,icompteur,iglen, & 
    134134                          ice,bm,bmelt,ablbord,ablbord_dtt,dt,    & 
    135135                          s,h,b,bsoc,flot,mk,mk0,uxbar,uybar,hwater,time,timemax,ndebug,ndebug_max 
    136   use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,runname,itracebug 
     136  use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,itracebug 
    137137  use geography, only: nx,ny,geoplace 
    138138  use deformation_mod_2lois, only:n1poly,n2poly 
    … …  
    156156  !      1 ou 2 passe (se met a 0 tout seul si on prend un fichier .cptr) 
    157157 
    158   ITEMP=0 
    159  
    160   !      switch couple physique faible =>  CP et CT independant T 
    161   !               0     pas de trait. vert.  A FAIRE           niveau L0 
    162   !               1     pas de couplage , faible physique      niveau L1 
    163   !               2     couplage, faible physique              niveau L2 
    164   !               3     couplage, physique complete sans CBT   niveau L3 
    165   !               4     idem 3 mais loi de def. Duval          niveau L4 
    166   ICOUPLE=4 
    167  
    168158  nt=-1   ! utilisee dans initialisation flottab 
    169159  !     sortie profile tous les dtprofile 
    … …  
    209199  tcpt:if (ICOMPTEUR.eq.0) then 
    210200 
    211  
    212  
    213         if ((ITEMP.eq.0).or.(ITEMP.eq.3)) then 
    214201           call masque(flot,mk,mk0,itracebug) 
    215202 
    … …  
    224211           !       ======================================================== 
    225212 
    226            if (ITEMP.eq.0) call lineartemp() 
     213           call lineartemp() 
    227214 
    228215           call bmelt_grounded  
    … …  
    244231           call diffusiv() 
    245232           call SIA_velocities() 
    246         endif 
    247233 
    248234  else  ! tcpt     on reprend un fichier compteur (ICOMPTEUR.eq.1) 

  • branches/GRISLIv3/SOURCES/mix-SIA-L1_mod.f90

    r416 r481  
    5656    !$ USE OMP_LIB 
    5757    use geography, only:nx,ny,nz 
    58     use module3D_phy, only:debug_3D,ubx,uby,num_rep_42,flgzmx,flgzmy,uxbar,uybar,ux,uy,uzr,ice,i_resolmeca,& 
     58    use module3D_phy, only:debug_3D,ubx,uby,flgzmx,flgzmy,uxbar,uybar,ux,uy,uzr,ice,i_resolmeca,& 
    5959         uxflgz,uyflgz,uxdef,uydef,isynchro,secyear 
    6060    use runparam, only:itracebug 

  • branches/GRISLIv3/SOURCES/readinput.f90

    r449 r481  
    182182       real,dimension(:,:),intent(inout) :: Tab 
    183183       character(*),intent(in) :: filename 
    184        real*8, dimension(:,:), pointer :: tabvar 
    185184     end subroutine lect_ncfile 
    186185  end interface !<<<<<<<<<<<<<<<<<<< 

  • branches/GRISLIv3/SOURCES/spinup_mod.f90

    r467 r481  
    161161  subroutine lect_vitbil_Lebrocq 
    162162 
    163     use module3D_phy, only: num_rep_42,mstream_mx,num_param,pi 
     163    use module3D_phy, only: num_rep_42,num_param,pi 
    164164    use io_netcdf_grisli, only: Read_Ncdf_var 
    165165    use geography, only: dirnameinp 

  • branches/GRISLIv3/SOURCES/steps_time_loop.f90

    r477 r481  
    1414subroutine step_time_loop         
    1515 
    16   use module3d_phy, only: ispinup,isynchro,timemax,time,flot,S,B,Bsoc,H,sealevel_2d,hmx,hmy,& 
     16  use module3d_phy, only: ispinup,isynchro,time,flot,S,B,Bsoc,H,sealevel_2d,hmx,hmy,& 
    1717       ice,bm,bmelt,ablbord,ablbord_dtt,dt, & 
    1818       uxbar,uybar,hdot,dtmax,dtmin,iout,flot_marais,debug_3D 

  • branches/GRISLIv3/SOURCES/tab-litho-0.3.f90

    r469 r481  
    7272 
    7373 
    74  
    75   !open(55,file='iso-cat') 
    76   !write(55,*) 'WE avec LBLOC sur 480 km' 
    77   !do i=-lbloc,lbloc 
    78   !   write(55,'25(i4,1x)') (int(WE(i,j)*1.e08),j=-lbloc,lbloc) 
    79   !end do 
    80   !close(55) 
    81  
    8274  ! sommation des elements de WE dans som pour normaliser 
    8375  ! la somme des enfoncements a 1 

  • branches/GRISLIv3/SOURCES/tracebug.f90

    r465 r481  
    2121  implicit none  
    2222 
    23   integer :: i,j 
    2423  character(len=*),intent(in):: chaine ! chaine a imprimer 
    2524 
    26 i=5 
    27 j=5 
     25  write(num_tracebug,*) chaine 
    2826 
    29 write(num_tracebug,*) chaine 
    30  
    31 !  write(num_tracebug,901) chaine,i,j,mk_init(i,j),S(i,j),H(i,j),Bsoc(i,j),B(i,j),pvi(i,j),pvm(i,j) 
    32 901   format(A50,' S0,S,H0,H,B0,B   ',3(i0,1x),2x, 4(f0.1,1x),2(es12.4,1x) ) 
    33  
    34 ! write(num_tracebug,898) chaine,i,j,ibase(i,j),mstream_mx(i,j), betamx(i,j), uxflgz(i,j), gzmx(i,j),ok_umat(i,j) 
    35 898   format(A50,' ibase, mstream,betamx,drag_mx    ',4(i0,1x),3x, 2(es12.4,1x), 2x,2(l5,1x)) 
    36  
    37 !  write(num_tracebug,896) chaine,i,j,Uxbar(i-1:i+1,j),Uxflgz(i-1:i+1,j),Uxbar(i-1:i+1,j)-Uxflgz(i-1:i+1,j),Hdot(i,j),dt  
    38 896   format(A50,'  i,j,Uxbar,Uxflgz,Uxdef ', 2(i0,1x),8(f0.3,1x)) 
    39  
    40 !  write(num_tracebug+1,897) chaine, eps(i,j), pvi(i,j),Hdot(i-1:i+1,j),betamx(i,j) 
    41 897   format(A50,' eps,pvi,Hdot,beta ', 6(es12.4,1x)) 
    42  
    43  
    44  
    45 ! (195,154) un point ile 
    46 ! (30,253) 
    47 !  write(num_tracebug,889) chaine,i,j,S(i,j),H(i,j),Bsoc(i,j),Hp(i,j),Mk_init(i,j),i_hp(i,j) 
    48 ! write(num_tracebug,895) chaine, i,j,i_HP(i,j),i_Hp0(i,j),Hp(i,j),Hp0(i,j),H(i,j) 
    49  
    50 !write(num_tracebug,894) chaine,i,j,H(i,j),uxbar(i,j),uxbar(i+1,j),pvi(i,j),eps(i,j) 
    51  
    52 !write(num_tracebug,*) flotmx(i-1:i+5,j),'     ',flot(i-1:i+5,j) 
    53  
    54  
    55 !  write(num_tracebug,887) chaine, ispinup,iter_beta  
    56 !  write(num_tracebug,891) chaine,Vcol_x(i,j),Vsl_x(i,j),uxbar(i,j),ubx(i,j),uxdef(i,j),betamx(i,j) 
    57 !  write(num_tracebug,880) chaine 
    58  
    59 ! write(num_tracebug,*) chaine,'xmin,ymin',xmin,ymin 
    60  
    61 ! S(30,253),S(178,121) 1 point ile + 1 point shelf 
    62    
    63 887   format(A50,'  ispinup =',i0,'    iter_beta =',i0) 
    64 888   format(A50,'    S(30,253)) = ',f0.3, '    S(178,121) = ',f0.3) 
    65 889   format(A50,'  SHBMk de ', 2(i0,1x),': '4(f0.3,1x),2(i0,1x)) 
    66 890   format(A50,'  S Ux Uy ', 2(i0,1x),': '6(f0.3,1x)) 
    67 891   format(A50,'  uxbal Vsl Ux uxb uxdef beta : '6(f0.2,1x)) 
    68 892   format(A50,'  Uxbar ', 2(i0,1x),': '(f0.3,1x)) 
    69 893   format(A50,'  betamx, drag_mx, beta_c  ', 2(i0,1x),': '4(f0.3,1x)) 
    70 894   format(A50,'  H, Uxbar, pvi,epsxx,eps    ', 2(i0,1x),': '3(f0.2,1x),2(es12.4,1x) ) 
    71 895   format(A50,'  i,j,ihp,ihp0,hp,hp0,h   ', 4(i0,1x),': '3(f0.2,1x) ) 
    7227end subroutine tracebug 

  • branches/GRISLIv3/SOURCES/util_recovery.f90

    r446 r481  
    1313 
    1414module util_recovery 
     15 
    1516  implicit none  
    1617  double precision    ,dimension(:),allocatable :: tab_time 
    … …  
    1819  real*8, dimension(:), pointer                 :: time_out 
    1920  logical                                       :: logic_out  
     21 
    2022contains 
    2123 
    … …  
    2628  subroutine init_recovery 
    2729 
    28     use module3d_phy, only: time 
    2930    use runparam, only: itracebug,dirsource 
    3031 
Note: See TracChangeset for help on using the changeset viewer.