Changeset 149 for codes/icosagcm/trunk/src/physics.f90

Timestamp:

04/03/13 12:05:12 (11 years ago)

Author:

sdubey

Message:

Added few new routines to read NC files and compute diagnostics to r145.

Few routines of dry physics including radiation module, surface process and convective adjustment in new routine phyparam.f90. dynetat to read start files for dynamics.

check_conserve routine to compute conservation of quatities like mass, energy etc.etat0_heldsz.f90 for held-suarez test case initial conditions.

new Key time_style=lmd or dcmip to use day_step, ndays like in LMDZ

File:

• codes/icosagcm/trunk/src/physics.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/physics.f90

@@ -8 +8 @@
   SUBROUTINE init_physics
   USE icosa
-  USE physics_dcmip_mod, init_physics_dcmip=>init_physics
+  USE physics_dcmip_mod,init_physics_dcmip=>init_physics
+  USE physics_dry_mod
   IMPLICIT NONE
     
@@ -18 +19 @@
       CASE ('dcmip')
         CALL init_physics_dcmip
+
+      CASE ('lmd')
+        CALL init_physics_dry
       
       CASE DEFAULT
-         PRINT*, 'Bad selector for variable physics <',physics_type, &
+         PRINT*, 'Bad selector for variable physics init <',physics_type, &
               '> options are <none>, <dcmip>,'
-         STOP
+
     END SELECT
     
   END SUBROUTINE init_physics
 
-  SUBROUTINE physics(it,f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
+  SUBROUTINE physics(it,jD_cur,jH_cur,f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
   USE icosa
+  USE physics_dry_mod
   USE physics_dcmip_mod, physics_dcmip=>physics
+  USE etat0_mod
+  USE etat0_heldsz_mod
   IMPLICIT NONE
     INTEGER, INTENT(IN)   :: it
+    REAL(rstd),INTENT(IN)::jD_cur,jH_cur
     TYPE(t_field),POINTER :: f_phis(:)
     TYPE(t_field),POINTER :: f_ps(:)
@@ -37 +45 @@
     TYPE(t_field),POINTER :: f_ue(:)
     TYPE(t_field),POINTER :: f_q(:)
+    LOGICAL:: firstcall,lastcall
     
     SELECT CASE(TRIM(physics_type))
       CASE ('none')
+
+        SELECT CASE(TRIM(etat0_type))
+        CASE('heldsz') 
+!       CALL transfert_request(f_ps,req_i1)
+!       CALL transfert_request(f_theta_rhodz,req_i1)
+!       CALL transfert_request(f_ue,req_e1_vect)
+!       CALL held_saurez(f_ps,f_theta_rhodz,f_ue) 
+        CASE DEFAULT
+        PRINT*,"NO PHYSICAL PACAKAGE USED" 
+        END SELECT
     
       CASE ('dcmip')
         CALL physics_dcmip(it, f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
+
+      CASE ('dry')
+        CALL physics_dry(it,jD_cur,jH_cur,f_phis, f_ps, f_theta_rhodz, f_ue, f_q)
       
       CASE DEFAULT
          PRINT*, 'Bad selector for variable physics <',physics_type, &
               '> options are <none>, <dcmip>,'
-         STOP
+  STOP
     END SELECT