Changeset 4202 for dynamico_lmdz/simple_physics/phyparam
- Timestamp:
- Dec 20, 2019, 2:40:30 PM (5 years ago)
- Location:
- dynamico_lmdz/simple_physics/phyparam
- Files:
-
- 1 deleted
- 4 edited
- 3 moved
-
param/coefdifv.F (modified) (1 diff)
-
param/inifrict.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/inifrict.F) (1 diff)
-
param/iniphyparam.F (modified) (1 diff)
-
param/interpol.F (deleted)
-
param/mellor_yamada.F (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/my_25.F)
-
param/paramdef.F (modified) (1 diff)
-
param/phyparam.F (modified) (1 diff)
-
physics/turbulence.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/physics/vdif_mod.F90) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
dynamico_lmdz/simple_physics/phyparam/param/coefdifv.F
r4193 r4202 3 3 $ pcdzv,pcdzh) 4 4 USE constlim, ONLY : ccdzh, cdzmin, cpgam, cdzconst, dgcdrag 5 USE vdif_mod, ONLY : lmixmin5 USE turbulence, ONLY : lmixmin 6 6 USE phys_const, ONLY : g 7 7 IMPLICIT NONE -
dynamico_lmdz/simple_physics/phyparam/param/inifrict.F90
r4196 r4202 1 SUBROUTINE inifrict(timestep) 2 USE comgeomfi, ONLY : ngridmax,nlayermx 3 USE constlim 4 USE phys_const 5 IMPLICIT NONE 6 c=======================================================================7 c 8 c Calcul des coefficients de ls diffusion verticale 9 c 10 c=======================================================================11 c-----------------------------------------------------------------------12 c Declarations: 13 c ------------- 14 15 16 c local: 17 c ------ 18 19 REAL dgrad,cpr,rl2,rrr,timestep 20 INTEGER l,ij 21 22 c-----------------------------------------------------------------------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 c!!! cpgam = 5.e-3*cpp 41 cpgam=0. 42 43 c-----------------------------------------------------------------------44 c coefficient de diffusion dans l'atmosphere: 45 c ------------------------------------------- 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 CONTINUE53 54 c-----------------------------------------------------------------------55 c couche limite de surface: 56 c ------------------------- 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 CONTINUE63 64 RETURN 65 END 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 -
dynamico_lmdz/simple_physics/phyparam/param/iniphyparam.F
r4192 r4202 14 14 USE planet, ONLY : coefir, coefvis 15 15 USE astronomy 16 USE vdif_mod, ONLY : lmixmin, emin_turb16 USE turbulence, ONLY : lmixmin, emin_turb 17 17 IMPLICIT NONE 18 18 -
dynamico_lmdz/simple_physics/phyparam/param/paramdef.F
r4190 r4202 1 1 SUBROUTINE paramdef(ngrid,rnatur,albedo,inertie,emissiv,z0) 2 2 USE comgeomfi, ONLY : lati,sinlat,coslat 3 USE vdif_mod, ONLY : lmixmin, emin_turb3 USE turbulence, ONLY : lmixmin, emin_turb 4 4 USE planet, ONLY : coefir, coefvis 5 5 USE astronomy -
dynamico_lmdz/simple_physics/phyparam/param/phyparam.F
r4201 r4202 13 13 USE planet 14 14 USE astronomy 15 USE vdif_mod, ONLY : vdif15 USE turbulence, ONLY : vdif 16 16 USE solar, ONLY : solang, zenang, mucorr 17 17 USE radiative_sw, ONLY : sw -
dynamico_lmdz/simple_physics/phyparam/physics/turbulence.F90
r4199 r4202 1 MODULE vdif_mod1 MODULE turbulence 2 2 3 3 #include "use_logging.h" … … 444 444 END SUBROUTINE vdif 445 445 446 END MODULE vdif_mod446 END MODULE turbulence
Note: See TracChangeset
for help on using the changeset viewer.
