source: codes/icosagcm/trunk/src/physics.f90 @ 155

MODULE physics_mod

  CHARACTER(LEN=255) :: physics_type="none"


CONTAINS

  SUBROUTINE init_physics
  USE icosa
  USE physics_dcmip_mod,init_physics_dcmip=>init_physics
  USE physics_dry_mod
  IMPLICIT NONE
    
    CALL getin("physics",physics_type)
    
    SELECT CASE(TRIM(physics_type))
      CASE ('none')
    
      CASE ('dcmip')
        CALL init_physics_dcmip

      CASE ('lmd')
        CALL init_physics_dry
      
      CASE DEFAULT
         PRINT*, 'Bad selector for variable physics init <',physics_type, &
              '> options are <none>, <dcmip>,'

    END SELECT
    
  END SUBROUTINE init_physics

  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(:)
    TYPE(t_field),POINTER :: f_theta_rhodz(:)
    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
    END SELECT
    
  END SUBROUTINE physics

END MODULE physics_mod