Changeset 1359 for trunk/NEMO/OPA_SRC/DIA/diawri.F90

Timestamp:

2009-03-31T14:36:28+02:00 (15 years ago)

Author:

smasson

Message:

first implementation of iom_put, see ticket:387

File:

• trunk/NEMO/OPA_SRC/DIA/diawri.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DIA/diawri.F90

@@ -27 +27 @@
    USE in_out_manager  ! I/O manager
    USE diadimg         ! dimg direct access file format output
+   USE iom
    USE ioipsl
 
@@ -411 +412 @@
          WRITE(numout,*) '~~~~~~ '
       ENDIF
+      
+!---------------------------------------------------------------------------
+
+      CALL iom_put("votemper",tn)
+      CALL iom_put("vosaline",sn)
+      CALL iom_put("sosstsst",tn(:,:,1))   ! sea surface temperature
+      CALL iom_put("sosaline",sn(:,:,1))   ! sea surface salinity
+#if defined key_dynspg_rl
+      CALL iom_put("sobarstf",bsfn)   ! barotropic streamfunction
+#else
+      CALL iom_put("sossheig",sshn)   ! sea surface height
+#endif
+
+      CALL iom_put("sowaflup",emp )   ! upward water flux
+      CALL iom_put("sowaflcd",emps)   ! c/d water flux
+      zw2d(:,:) = emps(:,:) * sn(:,:,1) * tmask(:,:,1)
+      CALL iom_put("sosalflx",zw2d)   ! c/d salt flux
+      CALL iom_put("sohefldo",qns + qsr )   ! total heat flux
+      CALL iom_put("soshfldo",qsr)   ! solar heat flux
+      CALL iom_put("somxl010",hmlp)   ! mixed layer depth
+      CALL iom_put("somixhgt",hmld)   ! turbocline depth
+      CALL iom_put("soicecov",fr_i)   ! ice fraction 
+#if ! defined key_coupled
+      CALL iom_put("sohefldp",qrp)   ! heat flux damping
+      CALL iom_put("sowafldp",erp)   ! freshwater flux damping
+      zw2d(:,:) = erp(:,:) * sn(:,:,1) * tmask(:,:,1)
+      CALL iom_put("sosafldp", zw2d)   ! salt flux damping
+#endif
+
+#if ( defined key_coupled && ! defined key_lim3 && ! defined key_lim2 ) 
+      CALL iom_put("sohefldp",qrp)   ! heat flux damping
+      CALL iom_put("sowafldp",erp)   ! freshwater flux damping
+         zw2d(:,:) = erp(:,:) * sn(:,:,1) * tmask(:,:,1)
+      CALL iom_put("sosafldp",zw2d)   ! salt flux damping
+#endif
+#if defined key_diaspr
+      CALL iom_put("sosurfps",gps)   ! surface pressure
+#endif
+         zw2d(:,:) = FLOAT( nmln(:,:) ) * tmask(:,:,1)
+      CALL iom_put("sobowlin",zw2d)   ! ???
+
+#if defined key_diahth
+      CALL iom_put("sothedep",hth)   ! depth of the thermocline
+      CALL iom_put("so20chgt",hd20)   ! depth of the 20 isotherm
+      CALL iom_put("so28chgt",hd28)   ! depth of the 28 isotherm
+      CALL iom_put("sohtc300",htc3)   ! first 300m heaat content
+#endif
+
+#if defined key_coupled 
+#  if defined key_lim3
+      Must be adapted for LIM3
+#  else
+      CALL iom_put("soicetem",tn_ice)   ! surf. ice temperature
+      CALL iom_put("soicealb",alb_ice)   ! ice albedo
+#  endif
+#endif
+         ! Write fields on U grid
+      CALL iom_put("vozocrtx",un)    ! i-current
+#if defined key_diaeiv
+      CALL iom_put("vozoeivu",u_eiv)    ! i-eiv current
+#endif
+      CALL iom_put("sozotaux",utau)   ! i-wind stress
+#if defined key_dynspg_rl
+      CALL lbc_lnk( spgu, 'U', -1. )
+      CALL iom_put("sozospgx",spgu)   ! i-surf. press. grad.
+#endif
+
+         ! Write fields on V grid
+      CALL iom_put("vomecrty",vn)   ! j-current
+#if defined key_diaeiv
+      CALL iom_put("vomeeivv",v_eiv)   ! j-eiv current
+#endif
+      CALL iom_put("sometauy", vtau)   ! j-wind stress
+#if defined key_dynspg_rl
+      CALL lbc_lnk( spgv, 'V', -1. )
+      CALL iom_put("somespgy",spgv)   ! j-surf. pressure grad.
+#endif
+
+         ! Write fields on W grid
+      CALL iom_put("vovecrtz",wn)    ! vert. current
+#   if defined key_diaeiv
+      CALL iom_put("voveeivw",w_eiv)    ! vert. eiv current
+#   endif
+      CALL iom_put("votkeavt",avt)    ! T vert. eddy diff. coef.
+      CALL iom_put("votkeevd",avt_evd)    ! T enhan. vert. eddy diff. coef.
+      CALL iom_put("votkeavm",avmu)    ! T vert. eddy visc. coef.
+      CALL iom_put("votkeevm",avmu_evd)    ! T enhan. vert. eddy visc. coef.
+      IF( lk_zdfddm ) THEN
+         CALL iom_put("voddmavs",fsavs(:,:,:) )    ! S vert. eddy diff. coef.
+      ENDIF
+#if defined key_traldf_c2d
+      CALL iom_put("soleahtw",ahtw)   ! lateral eddy diff. coef.
+# if defined key_traldf_eiv
+      CALL iom_put("soleaeiw",aeiw)   ! EIV coefficient at w-point
+# endif
+#endif
+!---------------------------------------------------------------------------
 
       ! Write fields on T grid