VARphy.F90

Go to the documentation of this file.

MODULE varphy

IMPLICIT NONE

! MAR Constants
                                                                              
      CHARACTER(LEN=1), PARAMETER :: labnum(0:9) =                            &
              (/'0','1','2','3','4','5','6','7','8','9'/)             
!C +...                labnum: Alphanumeric Character                          
                                                                               
      INTEGER, PARAMETER ::  izr=0,iun=1                                             
      REAL, PARAMETER    ::  zero = 0.0e+0, demi = 0.5e+0, unun = 1.0e+0,     &
     &                       epsi = 1.0e-6, eps9 = 1.0e-9, eps12= 1.0e-12,    &
     &                       third=1/3,     thous=1.e3             
      REAL               ::  argmin,argmax      !hj valeurs??
!C +...                    zero = 0.0e+0                                       
!C +                            demi = 0.5e+0                                 
!C +                                 unun = 1.0e+0                             
!C +                                      epsi = 1.0e-6                        
!C +                                           eps9 = 1.0e-9                   
!C +                                                eps12= 1.0e-12             
!C +                                                      third=1/3            
!C +                                                            thous=1.e3     
!C +                       argmin:        Function exp(x) Minimum Argument     
!C +                              argmax: Function exp(x) Maximum Argument     
!C +                                      (Min/Max are Machine Dependant)      
                                                                              
      REAL, PARAMETER    ::  pi = 3.141592653589793238462643e0,               &
     &                       degrad = pi / 180.d0, hourad = pi /  12.d0 
!C +...                pi    = 3.141592653589793238462643e0                   
!C +                   degrad= pi / 180.d0                                     
!C +                   hourad= pi /  12.d0                                    
                                                                               
      REAL, PARAMETER    ::  earthr = 6371.229e+3, earthv = 729.217e-7,       &
     &                       gravit = 9.81e0,      gravi2 = gravit**2 ,       &
     &                       grvinv = 1./gravit                       
!C +...                earthr: Earth Radius                  =6371.d+3    m    
!C +...                earthv: Earth Angular Velocity        = 7.29d-5  s-1   
!C +                   gravit: Earth Gravity Acceleration    = 9.81    m/s2    
!C +                   gravi2: idem (squared)                                  
!C +                   grvinv: idem (inverse)                                  
                                                                              
      REAL, PARAMETER    ::  akmol = 1.35e-5, vonkar = 0.40e0, a_turb = 5.8,  &
     &                       ahturb = 5.4,    asturb = 4.0,    r_turb = 3.0            
!C +...                akmol : Air Viscosity                 = 1.35d-5 m2/s   
!C +                   vonkar: von Karman constant           = 0.4            
!C +                   A_Turb: Stability  Coefficient Moment = 5.8             
!C +                   AhTurb: Stability  Coefficient Heat   = 5.4             
!C +                   AsTurb: Stability  Coefficient BLOW * = 4.0            
!C +                   r_Turb: Turbulent Diffusivities Ratio K*/Km             
                                                                              
!      REAL    ::  RDryAi,cp,CvDrya,racv,cap,pcap,stefan,qv_MIN           
      REAL, PARAMETER    ::  rdryai=287.05967, ra=287.05967, cp=1004.708845,  &
     &                       cvdrya=cp-rdryai, racv=rdryai/cvdrya,            &
     &                       cap=0.28586e0,    pcap=3.73010e0,                &
     &                       stefan=5.67e-8,   qv_min=3.00e-6          
      REAL, PARAMETER    ::  rvapor=461.e0,    epsq=3.00e-6, lv_h2o=2.5008e+6,&
     &                       ls_h2o=2.8345e+6, r_lvcp=2490.04,r_lccp=332.27,  &
     &                       r_lscp=2822.31            
!C +...                RDryAi: perfect gas law constant for dry air (287 J/kg/K) 
!C +                   ra    : perfect gas law constant for dry air (287 J/kg/K) 
!C +                   cp    : dry air specific heat at constant p (1004 J/kg/K) 
!C +                   CvDryA: dry air specific heat at constant V ( 717 J/kg/K) 
!C +                   racv  = ra / cv      = ra /(cp-ra)                        
!C +                   cap   = ra / cp                                           
!C +                   pcap  = 100.[kPa]**cap                                    
!C +...                RVapor: perfect gas law constant Water Vapor (461 J/kg/K) 
!C +                   epsq  : Minimum Water Vapor Content    =    3.00d-6 kg/kg 
!C +                   Lv_H2O: Latent Heat Vaporisation/Water = 2500.00d+3 J/kg  
!C +                   Ls_H2O: Latent Heat Sublimation /Ice   = 2833.60d+3 J/kg  
!C +                   r_LvCp [Lv=2500000J/kg]/[Cp=1004J/kg/K]= 2490.04    K/[kg/
!C +                   r_LvCp [Lc= 333600J/kg]/[Cp=1004J/kg/K]=  332.27    K/[kg/
!C +                   r_LvCp [Ls=2833600J/kg]/[Cp=1004J/kg/K]= 2822.31    K/[kg/
!C +                   stefan: Stefan Constant                =    5.67d-8 W/m2/K
!C +                   qv_MIN: Minimum Specific Humidity      =    3.00d-6 kg/kg 
                                                                                
      REAL, PARAMETER    ::  tfrwat=271.20e0,  siclf=302.00e6,cdsice=2.04e0,  &
     &                       hic0=1.00e0, ro_wat=1000.00e0, c__wat=4186.00e0, &
     &                       fracoh=0.75e0             
!C +...                tfrwat: Sea-Water freezing Temperature =  271.20d+0 K     
!C +                   siclf : Sea-Ice Heat of Fusion         =  302.00d+6 J/m3  
!C +                   cdsice: Sea-Ice Thermal Conductivity   =    2.04d+0 W/mK  
!C +                   hic0  : Sea-Ice Initial Thickness      =    2.00d+0 m     
!C +                   ro_Wat: Density       of Water         = 1000.00d+0 kg/m3 
!C +                   C__Wat: Heat Capacity of Water         = 4186.00d+0 J/kg/K
!C +                   fracoh: 25% of cooling = Oceanic Heat Flux                
!C +***                        (Hibler 1984)    (ANTARCTIC Ocean)                
                                                                                
      REAL, PARAMETER    ::  ro_ice=920.e0, cdice=2.10e0                                 
!C +...                ro_Ice:  Density      of Pure Ice       =  920.00d+0 kg/m3
!C +                   cdice :  Conductivity of Pure Ice       =    2.51d+0 W/m/K
                                                                                
      REAL, PARAMETER    ::  tfsnow=273.15e+0,csnow=2105.00e+0,r0sno=3.00e+1, &
     &                       blsno=3.30e+2,   lf_h2o=3.337e+5                          
!C +...                TfSnow:        Snow melting Temperature=  273.15d+0 K     
!C +                    csnow:Heat Capacity of Snow             2105      J/kg/K 
!C +                         (Loth et al. 1993, JGR 98 D6, 2.2.2 2e para p.10453)
!C +                   r0sno : Fresh  Snow Density                50.00d+0 kg/m3 
!C +                   blsno : Blowed Snow Density               330.00d+0 kg/m3 
!C +                   Lf_H2O: Latent Heat of Fusion of Snow  =    3.34d+5 J/kg  

! MAR Constants
                      
                                                                                

END MODULE varphy