source: dynamico_lmdz/simple_physics/phyparam/param/inifrict.F90 @ 4211

Last change on this file since 4211 was 4202, checked in by dubos, 6 years ago

simple_physics : cleanup
File size: 3.1 KB
Line 
1      SUBROUTINE inifrict(timestep)
2      USE comgeomfi, ONLY : ngridmax,nlayermx
3      USE constlim
4      USE phys_const
5      IMPLICIT NONE
6!=======================================================================
7!                                                                       
8!   Calcul des coefficients de ls diffusion verticale                   
9!                                                                       
10!=======================================================================
11!-----------------------------------------------------------------------
12!   Declarations:                                                       
13!   -------------                                                       
14                                                                       
15                                                                       
16!   local:                                                             
17!   ------                                                             
18                                                                       
19      REAL dgrad,cpr,rl2,rrr,timestep
20      INTEGER l,ij
21                                                                       
22!-----------------------------------------------------------------------
23                                                                       
24      call initconstlim
25                                                                       
26      dtradia=timestep
27      PRINT*,'DTPHYS',dtradia
28      lmixmin=100.
29      ais1    = 1.
30      ais2    = ais1 - 1.
31      print*,'ais1',ais1,'ais2',ais2
32      cdzmin  = 1.e-6
33      OPEN(99,file='cdzmin',status='old',err=9999)
34      READ(99,*) cdzmin
35 9999 CLOSE(99)
36      PRINT*,'cdzmin=',cdzmin
37                                                                       
38      cpr     = cpp/ r
39      ccdzh   = 2.5*g
40!!!!  cpgam   = 5.e-3*cpp                                               
41      cpgam=0.
42                                                                       
43!-----------------------------------------------------------------------
44!   coefficient de diffusion dans l'atmosphere:                         
45!   -------------------------------------------                         
46                                                                       
47      rl2=lmixmin**2
48      cdzconst(1)= 0.
49      DO 15 l=1,nlayermx - 1
50         cdzconst(l+1)= dtradia*g*g*cpr*rl2
51         print*,'cdzconst(',l+1,')  =  ',cdzconst(l+1)
52   15 END DO
53                                                                       
54!-----------------------------------------------------------------------
55!   couche limite de surface:                                           
56!   -------------------------                                           
57                                                                       
58      cdrat   = 2.e-3
59      dgrad   = dtradia*g*cpp/r
60      DO 16 ij = 1, ngridmax
61         dgcdrag( ij ) = cdrat * dgrad
62   16 END DO
63                                                                       
64      RETURN
65      END                                           
Note: See TracBrowser for help on using the repository browser.