Changeset 468


Ignore:
Timestamp:
01/26/24 17:14:43 (4 months ago)
Author:
aquiquet
Message:

Cleaning branch: removing weird useless boolean

Location:
branches/GRISLIv3/SOURCES
Files:
7 edited

Legend:

Unmodified
Added
Removed

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

    r467 r468  
    225225 
    226226  ! ===================== Booleens ======================================== 
    227   logical :: SHELFY                     !<  
    228   logical :: MARINE                     !<  
    229  
    230227  logical,dimension(nx,ny) :: FLOT      !< vrai si flottant (test d'archimede) 'o' 
    231228  logical,dimension(nx,ny) :: FLOTMX    !< if flotting (one neighboord  floats) '>' 

  • branches/GRISLIv3/SOURCES/conserv-mass-adv-diff_sept2009_mod.f90

    r467 r468  
    2020  use module3D_phy, only: V_limit,num_param,num_rep_42,dx1,mk0,i_Hp,Hp,H,uxbar,uybar,testdiag,& 
    2121       dtmax,dt,dtmin,time,dtt,isynchro,diffmx,diffmy,sdx,sdy,hmx,hmy,flgzmx,flgzmy,flot,timemax,& 
    22        marine,dtdx2,dtdx,bm,bmelt,igrdline,ibmelt_inv,ice,ablbord,hdot 
     22       dtdx2,dtdx,bm,bmelt,igrdline,ibmelt_inv,ice,ablbord,hdot 
    2323  use geography, only: nx,ny,dx,geoplace 
    24   use runparam, only: itracebug,num_tracebug,tgrounded 
     24  use runparam, only: itracebug,num_tracebug 
    2525  use reso_adv_diff_2D_vect, only: init_reso_adv_diff_2D,resol_adv_diff_2D_vect 
    2626!  use prescribe_H, only:prescribe_fixed_points,prescribe_paleo_gl_shelf,prescribe_present_H_gl,& 
    … …  
    371371 
    372372    timemax=time 
    373  
    374     ! etait avant dans step 
    375     if (time.lt.TGROUNDED) then 
    376        MARINE=.false. 
    377     else 
    378        MARINE=.true. 
    379     endif 
    380     ! fin vient de step 
    381373 
    382374    ! les variables dtdx et dtdx2 sont globales 

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

    r467 r468  
    104104     
    105105    use runparam, only:itracebug,nt 
    106     use module3D_phy, only:shelfy,igrdline,mk_init,flot,H,sealevel_2d,Bsoc,S,H,B,& 
     106    use module3D_phy, only:igrdline,mk_init,flot,H,sealevel_2d,Bsoc,S,H,B,& 
    107107         ice,front,& 
    108108         iceberg,uxbar,uybar,mk,gzmx,gzmy,flotmx,flotmy,hmx,hmy,isynchro,ilemx,ilemy,& 
    109          flgzmx,flgzmy,marine,fbm,bm,bmelt,debug_3D,dt 
     109         flgzmx,flgzmy,fbm,bm,bmelt,debug_3D,dt 
    110110    use param_phy_mod, only:row,ro 
    111111    use module_choix, only:mstream_dragging ! subroutine qui depend du dragging 
    … …  
    120120     
    121121    if (itracebug.eq.1)  call tracebug(' Entree dans routine flottab') 
    122  
    123     SHELFY = .FALSE. 
    124122 
    125123 
    … …  
    235233 
    236234          if (flot(i,j)) then 
    237              shelfy = .true. 
    238  
    239235             surnet=H(i,j)*(1.-ro/row) 
    240236             S(i,j)=surnet+sealevel_2d(i,j) !afq -- WARNING: est-ce qu'on veut vraiment mettre S a la valeur locale du niveau marin? 
    … …  
    335331    end do domain_y 
    336332    !$OMP END DO 
    337  
    338  
    339     !------------------------------------------------------------------------------------- 
    340     ! attention : pour expériences Heino 
    341     !             gzmy(i,j)=gzmy_heino(i,j) 
    342     !             shelfy=.true. 
    343     !             shelfy=.false. 
    344     !____________________________________________________________________________________  
    345  
    346333 
    347334!!$ 
    … …  
    536523 
    537524    !$OMP WORKSHARE 
    538     flgzmx(:,:)=(marine.and.(flotmx(:,:).or.gzmx(:,:).or.ilemx(:,:)))   & 
    539          .or.(.not.marine.and.flotmx(:,:)) 
     525    flgzmx(:,:)=(flotmx(:,:).or.gzmx(:,:).or.ilemx(:,:)) 
    540526    where (hmx(:,:).eq.0.) 
    541527       flgzmx(:,:) = .false. 
    542528    endwhere 
    543     flgzmy(:,:)=(marine.and.(flotmy(:,:).or.gzmy(:,:).or.ilemy(:,:)))   & 
    544          .or.(.not.marine.and.flotmy(:,:)) 
     529    flgzmy(:,:)=(flotmy(:,:).or.gzmy(:,:).or.ilemy(:,:)) 
    545530    where (hmy(:,:).eq.0.) 
    546531       flgzmy(:,:) = .false. 

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

    r467 r468  
    3333 
    3434    use iso_declar,only: nlith,dt_iso,tausoc,dl,rl,lbloc,we,charge 
    35     use module3D_phy, only: icouple,marine,err,h0,bsoc0,sealevel_2d 
     35    use module3D_phy, only: icouple,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 
    … …  
    4444    integer :: i,j 
    4545 
    46     if (GEOPLACE.eq.'eismint') then 
    47        NBED=0 
    48        NLITH=0 
    49     else if (GEOPLACE(1:6).eq.'marine') then 
    50        NBED=1 
    51        NLITH=1 
    52     else if ((geoplace.eq.'anteis1').or.(geoplace.eq.'ant20km')) then 
    53  
    54        if (icouple.eq.2) then 
    55           NBED=0 
    56           NLITH=0 
    57           if (marine) then 
    58              NBED=1 
    59              NLITH=1 
    60           endif 
    61        else 
    62           NBED=1 
    63           NLITH=1 
    64        endif 
    65     else 
    66        NBED=1 
    67        !       switch sur la lithosphere 
    68        NLITH=1 
    69     endif 
     46    NBED=1 
     47    NLITH=1 
    7048 
    7149    !      temps de reaction en annees 

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

    r467 r468  
    131131subroutine grisli_init 
    132132 
    133   USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen,marine,num_sealevel, & 
     133  USE module3D_phy, only: itemp,icouple,isynchro,icompteur,iglen,num_sealevel, & 
    134134                          num_ts_ritz,num_ic_vo,num_ic_by,num_ic_dm,num_ic_dc,num_ic_df, & 
    135135                          s,h,b,bsoc,flot,mk,mk0,uxbar,uybar,hwater,time,timemax,ndebug,ndebug_max 
    136   use runparam, only: nt,tbegin,tgrounded,dtprofile,dtcpt,dirnameout,runname,itracebug 
     136  use runparam, only: nt,tbegin,dtprofile,dtcpt,dirnameout,runname,itracebug 
    137137  use geography, only: nx,ny,geoplace 
    138138  use deformation_mod_2lois, only:n1poly,n2poly 
    … …  
    169169  !     sortie profile tous les dtprofile 
    170170  DTPROFILE=50000. 
    171   marine=.true. 
    172171  !     ----------------------------------fin des modifs run les plus usuelles 
    173172  DIRNAMEOUT='../RESULTATS/' 
    … …  
    207206  if (iglob_ncdf .EQ. 1) call sortie_ncdf_cat 
    208207 
    209  
    210 !cdc supprime pour initialisation propre 
    211 !~   if (iter_beta.eq.0) then 
    212  
    213 !~      if (itracebug.eq.1)  call tracebug(' Avant appel routine icethick3') 
    214 !~      call icethick3 
    215 !~      debug_3D(:,:,88) = S(:,:) 
    216 !~      if (itracebug.eq.1)  call tracebug(' Apres appel routine icethick3') 
    217 !~   end if 
    218  
    219  
    220   !     Tgrounded, temps pendant lequel la calotte est terrestre 
    221   tgrounded=tbegin-10. 
    222   !if (tgrounded.le.tbegin) then  
    223   marine=.true. ! Cas la calotte est terrestre  
    224   !end if 
    225  
    226208  ! test vincent car certains H(i,j)=0 dans fichier de reprise 
    227209  do j=1,ny 
    … …  
    232214 
    233215 
    234   ! call firstoutput()           ! ouverture fichier temporel et premieres ecritures 
    235  
    236216  call forclim                   !  initialisation BM et TS          
    237217  call ablation 
    … …  
    242222 
    243223  tcpt:if (ICOMPTEUR.eq.0) then 
    244  
    245  
    246      if ((GEOPLACE.ne.'eismint').and.(GEOPLACE(1:6).ne.'marine')) then 
    247         !       ITEMP=1 => calcul de T lliboutry; ITEMP=2 => reprise d'un fichier cptr 
    248         !       ITEMP=0 => on ne prend pas en compte T Lliboutry 
    249         !       ITEMP=3 => on prend les temperatures d'un fichier cptr 
    250224 
    251225 
    … …  
    285259           call SIA_velocities() 
    286260        endif 
    287  
    288  
    289      endif 
    290      !     fin du test geoplace 
    291261 
    292262  else  ! tcpt     on reprend un fichier compteur (ICOMPTEUR.eq.1) 

  • branches/GRISLIv3/SOURCES/runparam_mod.f90

    r432 r468  
    2525  integer :: itracebug                  !< si 1 fait une impression au début des routines 
    2626  integer :: num_tracebug               !< numero de l'unite itracebug 
    27   real tgrounded 
    2827  real dtsortie 
    2928  real dtcpt 

  • branches/GRISLIv3/SOURCES/steps_time_loop.f90

    r465 r468  
    166166subroutine step_thermomeca 
    167167   
    168   use module3d_phy, only: ispinup,isynchro,timemax,time,marine,iglen,shelfy,icompteur,& 
     168  use module3d_phy, only: ispinup,isynchro,timemax,time,iglen,icompteur,& 
    169169       inv_beta,dtmin,flot,mk,mk0 
    170170  use runparam, only: itracebug,nt,num_tracebug 
    … …  
    196196       '   nt=',nt 
    197197 
    198   !-------------------------------------------------------------------- 
    199   !     lecture dans un petit fichier "runname.limit" au cas ou les 
    200   !     parametres du run seraient changes   
    201   !                                                    
    202   !-------------------------------------------------------------------- 
    203  
    204   !  call limit_file(nt,real(time),dt,tend,dtsortie,dtcpt,testdiag,dtt,runname) 
    205  
    206   !if (time.lt.tgrounded) then 
    207   !   marine=.false. 
    208   !else 
    209   marine=.true. 
    210   !endif 
    211  
    212198  ! update the regions (floating, ice ...) 
    213199  !------------------------------------------- 
    … …  
    348334     !-------------   ------------------------------------------------------- 
    349335 
    350      isync_2:    if ((shelfy).and.(marine).and.(isynchro.eq.1)) then                 
     336     isync_2:    if (isynchro.eq.1) then                 
    351337 
    352338        !les 3 lignes ci-dessous pour avoir les champs avant diagno (en cas de pb colineaire) 
Note: See TracChangeset for help on using the changeset viewer.