Changeset 1181 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90

Timestamp:

Jun 15, 2009, 6:49:35 PM (15 years ago)

Author:

jghattas

Message:

• Ajoute le possibilite de forcer avec les aerosols poussiere et sel de mer.
• Ajoute le module aerosol_mod qui contient des parmetres pour les aerosols
• flag_aerosol a legerement change signification

A Cozic

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90 (modified) (7 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90

@@ -16 +16 @@
   USE mod_phys_lmdz_para, ONLY : mpi_rank  
   USE readaerosol_mod
+  USE aerosol_mod, ONLY : naero_spc, name_aero
   USE write_field_phy
 
@@ -49 +50 @@
   INTEGER, SAVE                   :: klev_src    ! Number of vertical levles in source field
 !$OMP THREADPRIVATE(klev_src)
-  INTEGER, PARAMETER              :: nb_aero=8
-
-  CHARACTER(len=7),DIMENSION(nb_aero)   :: name_aero
 
   REAL, DIMENSION(klon)             :: psurf_day, pi_psurf_day
@@ -77 +75 @@
   LOGICAL            :: lnewday      ! Indicates if first time step at a new day
   LOGICAL,SAVE       :: vert_interp  ! Indicates if vertical interpolation will be done
-  LOGICAL,SAVE       :: debug=.TRUE. ! Debugging in this subroutine
+  LOGICAL,SAVE       :: debug=.FALSE.! Debugging in this subroutine
 !$OMP THREADPRIVATE(vert_interp, debug)
 
@@ -85 +83 @@
 !
 !****************************************************************************************
-! Variable containing aerosols name
-  name_aero(1) = "SSSSM  "
-  name_aero(2) = "ASSSM  "
-  name_aero(3) = "ASBCM  "
-  name_aero(4) = "ASPOMM "
-  name_aero(5) = "SO4    "
-  name_aero(6) = "CIDUSTM"
-  name_aero(7) = "AIBCM  "
-  name_aero(8) = "AIPOMM "
 
 ! Calculation to find if it is a new day
@@ -106 +95 @@
 
   IF (.NOT. ALLOCATED(var_day)) THEN
-     ALLOCATE( var_day(klon, klev, nb_aero), stat=ierr)
+     ALLOCATE( var_day(klon, klev, naero_spc), stat=ierr)
      IF (ierr /= 0) CALL abort_gcm('readaerosol_interp', 'pb in allocation 1',1)
-     ALLOCATE( pi_var_day(klon, klev, nb_aero), stat=ierr)
+     ALLOCATE( pi_var_day(klon, klev, naero_spc), stat=ierr)
      IF (ierr /= 0) CALL abort_gcm('readaerosol_interp', 'pb in allocation 2',1)
 
-     ALLOCATE( psurf_year(klon, 12, nb_aero), pi_psurf_year(klon, 12, nb_aero), stat=ierr)
+     ALLOCATE( psurf_year(klon, 12, naero_spc), pi_psurf_year(klon, 12, naero_spc), stat=ierr)
      IF (ierr /= 0) CALL abort_gcm('readaerosol_interp', 'pb in allocation 3',1)
 
-     ALLOCATE( load_year(klon, 12, nb_aero), pi_load_year(klon, 12, nb_aero), stat=ierr)
+     ALLOCATE( load_year(klon, 12, naero_spc), pi_load_year(klon, 12, naero_spc), stat=ierr)
      IF (ierr /= 0) CALL abort_gcm('readaerosol_interp', 'pb in allocation 4',1)
 
@@ -136 +125 @@
           psurf_year(:,:,id_aero), load_year(:,:,id_aero))
      IF (.NOT. ALLOCATED(var_year)) THEN
-        ALLOCATE(var_year(klon, klev_src, 12, nb_aero), stat=ierr)
+        ALLOCATE(var_year(klon, klev_src, 12, naero_spc), stat=ierr)
         IF (ierr /= 0) CALL abort_gcm('readaerosol_interp', 'pb in allocation 5',1)
      END IF
@@ -154 +143 @@
 
      IF (.NOT. ALLOCATED(pi_var_year)) THEN
-        ALLOCATE(pi_var_year(klon, klev_src, 12, nb_aero), stat=ierr)
+        ALLOCATE(pi_var_year(klon, klev_src, 12, naero_spc), stat=ierr)
         IF (ierr /= 0) CALL abort_gcm('readaerosol_interp', 'pb in allocation 6',1)
      END IF