Changeset 1052 for codes/icosagcm/devel/src/diagnostics/observable.f90

Timestamp:

09/23/20 16:00:58 (4 years ago)

Author:

dubos

Message:

devel : diagnose divergence and vorticity

File:

• codes/icosagcm/devel/src/diagnostics/observable.f90 (modified) (3 diffs)

codes/icosagcm/devel/src/diagnostics/observable.f90

@@ -46 +46 @@
     USE compute_temperature_mod
     USE compute_velocity_mod
+    USE compute_vorticity_mod
+    USE compute_divergence_mod
     USE vertical_interp_mod
     USE theta2theta_rhodz_mod
@@ -109 +111 @@
 
        CALL output_field("temp_init",f_buf_i)
+       CALL output_field("vort_init",f_buf_v)
     ELSE
        CALL output_field("p",f_pmid)
@@ -123 +126 @@
        CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,preff)
        CALL output_field("SST",f_buf_s)       
+    END IF
+
+    CALL vorticity(f_u, f_buf_v)
+    CALL divergence(f_u, f_buf_i)
+    IF(init) THEN
+       CALL output_field("vort_init",f_buf_v)
+       CALL output_field("div_init",f_buf_i)
+    ELSE
+       CALL output_field("vort",f_buf_v)
+       CALL output_field("div",f_buf_i)
     END IF