Changeset 321 for codes

Timestamp:

01/22/15 01:26:02 (9 years ago)

Author:

dubos

Message:

Partial fix for etat0 : W91 works again

File:

• codes/icosagcm/trunk/src/etat0.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/etat0.f90

@@ -112 +112 @@
     TYPE(t_field),POINTER :: f_phis(:)
     TYPE(t_field),POINTER :: f_theta_rhodz(:)
-    TYPE(t_field),POINTER :: f_temp(:)
     TYPE(t_field),POINTER :: f_u(:)
     TYPE(t_field),POINTER :: f_q(:)
   
+    TYPE(t_field),POINTER,SAVE :: f_temp(:)
     REAL(rstd),POINTER :: ps(:)
     REAL(rstd),POINTER :: mass(:,:)
@@ -124 +124 @@
     REAL(rstd),POINTER :: q(:,:,:)
     INTEGER :: ind
+
+    CALL allocate_field(f_temp,field_t,type_real,llm,name='temp')
 
     DO ind=1,ndomain
@@ -146 +148 @@
     IF( TRIM(etat0_type)/='williamson91.6' ) CALL temperature2theta_rhodz(f_ps,f_temp,f_theta_rhodz)
     
+    CALL deallocate_field(f_temp)
     
   END SUBROUTINE etat0_collocated
@@ -162 +165 @@
     REAL(rstd),INTENT(OUT) :: q(iim*jjm,llm,nqtot)
 
-    REAL(rstd) :: lon_i(iim*jjm)
-    REAL(rstd) :: lat_i(iim*jjm)
     REAL(rstd) :: ulon_i(iim*jjm,llm)
     REAL(rstd) :: ulat_i(iim*jjm,llm)
 
-    REAL(rstd) :: lon_e(3*iim*jjm)
-    REAL(rstd) :: lat_e(3*iim*jjm)
     REAL(rstd) :: ps_e(3*iim*jjm)
     REAL(rstd) :: mass_e(3*iim*jjm,llm)