Changeset 4242 for dynamico_lmdz/simple_physics/phyparam/DYNAMICO

Timestamp:

Jun 5, 2020, 1:59:19 PM (5 years ago)

Author:

dubos

Message:

simple_physics : output SW fluxes

File:

• dynamico_lmdz/simple_physics/phyparam/DYNAMICO/icosa_phyparam_mod.F90 (modified) (3 diffs)

dynamico_lmdz/simple_physics/phyparam/DYNAMICO/icosa_phyparam_mod.F90

@@ -199 +199 @@
       CALL allocate_field(f_write_llmp1, field_t, type_real, llm+1, name='phyparam_write_llmp1')
       writefield1_plugin => plugin_writefield1
+      writefield2_plugin => plugin_writefield2
     END SUBROUTINE init_plugin_writefield
 
@@ -206 +207 @@
       CHARACTER(*), INTENT(IN) :: name, longname, unit
       REAL, INTENT(IN)         :: var(:)
-      WRITELOG(*,*) TRIM(name), ' : ', TRIM(longname), MINVAL(var), MAXVAL(var), inout%it
+      WRITELOG(*,*) TRIM(name), ' : ', TRIM(longname), SHAPE(var), inout%it
+      WRITELOG(*,*) TRIM(name), ' : ', MINVAL(var), MAXVAL(var)
       LOG_INFO('writefield1')
       CALL unpack_field(f_write2d, var)
@@ -213 +215 @@
 
     SUBROUTINE plugin_writefield2(name,longname,unit, var)
+      USE physics_interface_mod, ONLY : unpack_field, inout => physics_inout
+      USE output_field_mod, ONLY : output_field
+      USE icosa,    ONLY : llm
       CHARACTER(*), INTENT(IN) :: name, longname, unit
       REAL, INTENT(IN)         :: var(:,:)
+      INTEGER :: nlev
+      WRITELOG(*,*) TRIM(name), ' : ', TRIM(longname), SHAPE(var), inout%it
+      WRITELOG(*,*) TRIM(name), ' : ', MINVAL(var), MAXVAL(var)
+      LOG_INFO('writefield2')
+      nlev = SIZE(var, 2)
+      IF(nlev==llm) THEN
+         CALL unpack_field(f_write_llm, var)
+         CALL output_field('phyparam_'//TRIM(name), f_write_llm)
+      ELSEIF(nlev==llm+1) THEN
+         CALL unpack_field(f_write_llmp1, var)
+         CALL output_field('phyparam_'//TRIM(name), f_write_llmp1)
+      END IF
     END SUBROUTINE plugin_writefield2