Changeset 4240

Timestamp:

May 21, 2020, 2:22:24 AM (5 years ago)

Author:

dubos

Message:

simple_physics : physics output with DYNAMICO, tested on Camelot

File:

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

dynamico_lmdz/simple_physics/phyparam/DYNAMICO/icosa_phyparam_mod.F90

@@ -4 +4 @@
 ! FCM gets confused when external modules are USEd at module level
 ! => USE statements to DYNAMICO modules go into subroutines
+  
+  USE icosa, ONLY : t_field
 
   IMPLICIT NONE
@@ -14 +16 @@
   REAL, PARAMETER :: oneday = 86400. ! hard-coded                                                                                                                                                         
   INTEGER, PARAMETER :: log_unit = 15
+
+  TYPE(t_field),POINTER :: f_write2d(:), f_write_llm(:), f_write_llmp1(:)
 
   PUBLIC :: init_physics, physics
@@ -40 +44 @@
       flush_plugin => flush_log_
 
+      CALL init_plugin_writefield
+
       physics_log_level='INF'
       CALL getin('physics_log_level', physics_log_level)
@@ -94 +100 @@
       END DO
 
-      IF(is_mpi_master) PRINT *, 'phyparam phi :', pphi(inout%ngrid/2+1, :)
+!      IF(is_mpi_master) PRINT *, 'phyparam phi :', pphi(inout%ngrid/2+1, :)
 
       CALL check_NaN('physics', 'ulon', inout%ulon)
@@ -108 +114 @@
            &        inout%dulon, inout%dulat, inout%dtemp, dps)
 
-      IF(is_mpi_master) PRINT *, 'phyparam dT :', inout%dtemp(inout%ngrid/2+1, :)
+!      IF(is_mpi_master) PRINT *, 'phyparam dT :', inout%dtemp(inout%ngrid/2+1, :)
 
       CALL check_NaN('physics', 'dulon', inout%dulon)
@@ -140 +146 @@
       INTEGER(c_int), INTENT(IN), VALUE :: lev, taglen, buflen, bufsize
       CHARACTER(KIND=c_char), INTENT(IN) :: tag(taglen), buf(buflen, bufsize)
+
       CHARACTER(buflen+1) :: line
-      !    SUBROUTINE flush_log_(lev, tag, buf)
-      !      INTEGER, INTENT(IN) :: lev
-      !      CHARACTER(*), INTENT(IN) :: tag, buf(:)
       CHARACTER(100) :: prefix
       INTEGER :: i
@@ -151 +155 @@
          DO i=1, bufsize
             WRITE(line,*) buf(:,i)
-            WRITE(log_unit,*) TRIM(prefix), TRIM(line)
+            WRITE(log_unit,*) TRIM(prefix) // TRIM(line)
          END DO
          WRITE(log_unit, *) ''
@@ -186 +190 @@
     END SUBROUTINE read_paramb
 
+!--------------------------------------- writefield ------------------------------------
+
+    SUBROUTINE init_plugin_writefield
+      USE icosa, ONLY : t_field, field_t, type_real, allocate_field, llm
+      USE writefield_mod, ONLY : writefield1_plugin, writefield2_plugin
+      CALL allocate_field(f_write2d,     field_t, type_real,        name='phyparam_write2d')
+      CALL allocate_field(f_write_llm,   field_t, type_real, llm,   name='phyparam_write_llm')
+      CALL allocate_field(f_write_llmp1, field_t, type_real, llm+1, name='phyparam_write_llmp1')
+      writefield1_plugin => plugin_writefield1
+    END SUBROUTINE init_plugin_writefield
+
+    SUBROUTINE plugin_writefield1(name,longname,unit, var)
+      USE physics_interface_mod, ONLY : unpack_field, inout => physics_inout
+      USE output_field_mod, ONLY : output_field
+      CHARACTER(*), INTENT(IN) :: name, longname, unit
+      REAL, INTENT(IN)         :: var(:)
+      WRITELOG(*,*) TRIM(name), ' : ', TRIM(longname), MINVAL(var), MAXVAL(var), inout%it
+      LOG_INFO('writefield1')
+      CALL unpack_field(f_write2d, var)
+      CALL output_field('phyparam_'//TRIM(name), f_write2d)
+    END SUBROUTINE plugin_writefield1
+
+    SUBROUTINE plugin_writefield2(name,longname,unit, var)
+      CHARACTER(*), INTENT(IN) :: name, longname, unit
+      REAL, INTENT(IN)         :: var(:,:)
+    END SUBROUTINE plugin_writefield2
+
 END MODULE icosa_phyparam_mod