Changeset 443 for branches


Ignore:
Timestamp:
07/03/23 17:00:03 (11 months ago)
Author:
aquiquet
Message:

ESIA and ESSA can be correlated by a simple function, activated in the param_list

Location:
branches/GRISLIv3
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • branches/GRISLIv3/Param/anteis1_param_list_Cyc-Scho.dat

    r235 r443  
    214214 sf01           =   0.1250000     
    215215 sf03           =   0.1250000     
     216 sffunc         =   0 
    216217 pvimin         =    1500.000     
    217218 / 
    … …  
    219220! sf01 : coefficient viscosite loi lineaire  
    220221! sf03 : coefficient viscosite loi n=3  
     222! sffunc : 1 for visco coef that depends on enh factor 
    221223! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3 
    222224! tres petit par rapport aux valeurs standards ~ 1.e10 

  • branches/GRISLIv3/Param/anteis1_param_list_Cyc-Tsai.dat

    r235 r443  
    214214 sf01           =   0.1250000     
    215215 sf03           =   0.1250000     
     216 sffunc         =   0 
    216217 pvimin         =    1500.000     
    217218 / 
    … …  
    219220! sf01 : coefficient viscosite loi lineaire  
    220221! sf03 : coefficient viscosite loi n=3  
     222! sffunc : 1 for visco coef that depends on enh factor 
    221223! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3 
    222224! tres petit par rapport aux valeurs standards ~ 1.e10 

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

    r401 r443  
    22module diagno_mod                   ! nom pendant les tests 
    33 !$ USE OMP_LIB 
    4 use module3D_phy,only: nx,ny,dx 
     4use module3D_phy 
     5use module_choix 
    56      
    67implicit none 
    … …  
    4950subroutine init_diagno 
    5051 
    51 use module3D_phy,only: num_rep_42,num_param,geoplace,pvimin 
    52  
    53 implicit none 
    54  
    55 namelist/diagno_rheol/sf01,sf03,pvimin 
     52integer sffunc 
     53 
     54namelist/diagno_rheol/sf01,sf03,sffunc,pvimin 
    5655 
    5756! attribution des coefficients de viscosite 
    … …  
    7170write(num_rep_42,*) 'sf01           = ',sf01  
    7271write(num_rep_42,*) 'sf03           = ',sf03 
     72write(num_rep_42,*) 'sffunc         = ',sffunc 
    7373write(num_rep_42,*) 'pvimin         = ',pvimin 
    7474write(num_rep_42,*)'/'      
    … …  
    7676write(num_rep_42,428) '! sf01 : coefficient viscosite loi lineaire ' 
    7777write(num_rep_42,428) '! sf03 : coefficient viscosite loi n=3 ' 
     78write(num_rep_42,428) '! sffunc : 1 for visco coef that depends on enh factor' 
    7879write(num_rep_42,428) '! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3' 
    7980write(num_rep_42,428) '! tres petit par rapport aux valeurs standards ~ 1.e10' 
    … …  
    8182write(num_rep_42,*) 
    8283 
     84if (sffunc.eq.1) then 
     85     sf01 = (1 - (sf(2)-1)/10.)/sf(2)  !the sf are known via the module_choix... 
     86     sf03 = (1 - (sf(1)-1)/10.)/sf(1) 
     87endif 
    8388 
    8489!      Precision utilisee dans de calcul 
    … …  
    100105subroutine diagnoshelf !      Resolution numerique des equations diagnostiques 
    101106 
    102   use module3D_phy, only: itracebug,niter_nolin,taushelf,flot,h,pvi,coef_sflot,rog, & 
    103                           gr_line_schoof,geoplace,uxbar,uybar,pvimin,debug_3D,V_limit, & 
    104                           imx_diag,imy_diag,tobmx,tobmy,neffmx,neffmy,gr_line, schoof, & 
    105                           betamx,betamy,flgzmx,flgzmy,uxflgz,uyflgz,taub 
    106   use furst_schoof_mod, only: back_force_x,back_force_y,gr_select,frot_coef,interpol_glflux 
    107   use module_choix, only: btt,dragging,mstream_dragging  ! afq - 03/2023 - not clean 
    108  
    109   implicit none 
    110107 
    111108  integer :: diagno_grline 
    112109 
    113   integer itour_nolin  
    114  
    115   integer i,j ! loop integers 
    116  
     110   
    117111  if (itracebug.eq.1)  call tracebug(' Entree dans diagnoshelf') 
    118112 
    … …  
    366360subroutine calc_pvi 
    367361 
    368 use module3d_phy, only: nx,ny,nz,itracebug,pvi,pvimin,abar,flot,gzmx,gzmy, & 
    369                         ilemx,ilemy,eps,taushelf,h,ramollo,debug_3d,pvm 
    370 use module_choix, only: btt,visc  ! afq - 03/2023 - not clean 
    371  
    372 implicit none 
    373  
    374362! calcule les viscosites integrees pvi et pvm 
    375363! loi polynomiale + couplage thermomécanique 
    … …  
    397385!$  integer :: rang ,nb_taches 
    398386!$  logical :: paral 
    399  
    400 integer i,j,k !loop integers 
    401387 
    402388if (itracebug.eq.1)  call tracebug(' Calc pvi') 
    … …  
    533519subroutine imx_imy_nx_ny 
    534520 
    535 use module3d_phy, only: imx_diag,imy_diag,flgzmx,flgzmy 
    536  
    537 implicit none 
    538  
    539521! definition des masques 
    540522! pour rempli_L2 : calcule les masques imx et imy qui  
    … …  
    616598 
    617599subroutine imx_imy_nx_ny_reduce(choix) 
    618  
    619   use module3d_phy, only: flot,imx_diag,imy_diag,gr_line 
    620  
    621   implicit none 
    622600 
    623601!afq -- For the backforce computation we do not need to compute the velocities everywhere 
    … …  
    685663 
    686664subroutine mismip_boundary_cond 
    687  
    688  use module3d_phy, only: itracebug,imx_diag,imy_diag,uxbar,uybar 
    689  
    690  implicit none 
    691  
    692  if (itracebug.eq.1)  call tracebug(' Subroutine mismip_boundray_cond') 
     665if (itracebug.eq.1)  call tracebug(' Subroutine mismip_boundray_cond') 
    693666 
    694667! Condition pas de flux sur les bords nord et sud 

  • branches/GRISLIv3/lauching_area/degla_PMIP_euras20_param_list.dat

    r365 r443  
    149149sf01           =        0.125   ! 0.125 
    150150sf03           =        0.125   ! 0.125 
     151sffunc         =        0 
    151152pvimin         =        1.5e3   ! 1.e3 
    152153/   
    … …  
    154155! sf01 : coefficient viscosite loi lineaire  
    155156! sf03 : coefficient viscosite loi n=3  
     157! sffunc : 1 for visco coef that depends on enh factor 
    156158! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3 
    157159! tres petit par rapport aux valeurs standards ~ 1.e10 

  • branches/GRISLIv3/lauching_area/degla_PMIP_hemin40_param_list.dat

    r367 r443  
    149149sf01           =        0.125   ! 0.125 
    150150sf03           =        0.125   ! 0.125 
     151sffunc         =        0 
    151152pvimin         =        1.5e3   ! 1.e3 
    152153/   
    … …  
    154155! sf01 : coefficient viscosite loi lineaire  
    155156! sf03 : coefficient viscosite loi n=3  
     157! sffunc : 1 for visco coef that depends on enh factor 
    156158! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3 
    157159! tres petit par rapport aux valeurs standards ~ 1.e10 

  • branches/GRISLIv3/lauching_area/degla_PMIP_laure16_param_list.dat

    r365 r443  
    149149sf01           =        0.125   ! 0.125 
    150150sf03           =        0.125   ! 0.125 
     151sffunc         =        0 
    151152pvimin         =        1.5e3   ! 1.e3 
    152153/   
    … …  
    154155! sf01 : coefficient viscosite loi lineaire  
    155156! sf03 : coefficient viscosite loi n=3  
     157! sffunc : 1 for visco coef that depends on enh factor 
    156158! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3 
    157159! tres petit par rapport aux valeurs standards ~ 1.e10 
Note: See TracChangeset for help on using the changeset viewer.