Changeset 2740 for LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

Timestamp:

Jan 3, 2017, 6:00:54 PM (7 years ago)

Author:

lguez

Message:

Added output of nudging coefficients for temperature and humidity,
along those for the wind.

File:

• LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -2130 +2130 @@
     USE comconst_mod, ONLY: pi
     USE comvert_mod, ONLY: presnivs
+    use netcdf95, only: nf95_def_var, nf95_put_var
+    use netcdf, only: nf90_float
+
     IMPLICIT NONE
 
@@ -2149 +2152 @@
     INTEGER       :: nid, id_lonu, id_lonv, id_latu, id_latv, id_tim, id_lev
     INTEGER       :: vid_lonu,vid_lonv,vid_latu,vid_latv,vid_cu,vid_cv,vid_lev
-    INTEGER       :: vid_au,vid_av
+    INTEGER       :: vid_au,vid_av, varid_alpha_t, varid_alpha_q
     INTEGER, DIMENSION (3) :: dim3
     INTEGER, DIMENSION (4) :: dim4,count,start
     INTEGER                :: ierr, varid,l
-    REAL zu(ip1jmp1),zv(ip1jm)
+    REAL zu(ip1jmp1),zv(ip1jm), zt(iip1, jjp1), zq(iip1, jjp1)
     REAL, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: field_glo
     
@@ -2171 +2174 @@
     print*,'guide_out apres gather '
     CALL Gather_field_u(alpha_u,zu,1)
+    CALL Gather_field_u(alpha_t,zt,1)
+    CALL Gather_field_u(alpha_q,zq,1)
     CALL Gather_field_v(alpha_v,zv,1)
 
@@ -2207 +2212 @@
         ierr=NF_DEF_VAR(nid,"au",NF_FLOAT,2,(/id_lonu,id_latu/),vid_au)
         ierr=NF_DEF_VAR(nid,"av",NF_FLOAT,2,(/id_lonv,id_latv/),vid_av)
-
+        call nf95_def_var(nid, "alpha_T", nf90_float, (/id_lonv, id_latu/), &
+             varid_alpha_t)
+        call nf95_def_var(nid, "alpha_q", nf90_float, (/id_lonv, id_latu/), &
+             varid_alpha_q)
+        
         ierr=NF_ENDDEF(nid)
 
@@ -2232 +2241 @@
         ierr = NF_PUT_VAR_REAL(nid,vid_av,alpha_v)
 #endif
+        call nf95_put_var(nid, varid_alpha_t, zt)
+        call nf95_put_var(nid, varid_alpha_q, zq)
 ! --------------------------------------------------------------------
 ! Cr�ation des variables sauvegard�es