Changeset 1303 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jul 1, 2014, 2:44:17 AM (11 years ago)

Author:

bclmd

Message:

LMDZ Generic: slab ocean added for rcm1d

File:

• trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F

@@ -15 +15 @@
       use phyredem, only: physdem0,physdem1
       use comgeomphy, only: initcomgeomphy
+      use slab_ice_h, only: noceanmx
+
       implicit none
 
@@ -115 +117 @@
       REAL hice(1),totcloudfrac(1)
       REAL qzero1D   !initial water amount on the ground
+
+!     added by BC for ocean
+      real rnat(1)
+      REAL tslab(1,noceanmx),tsea_ice(1),sea_ice(1)
+      real pctsrf_sic(1)
+
+
 
 !     added by AS to avoid the use of adv trac common
@@ -718 +727 @@
 
 
+! Initialize slab ocean
+! -----------------
+      rnat=1. ! default value for rnat
+      if(inertiedat(1,1).GE.10000.)then
+         rnat=0.
+      endif
+      if(ok_slab_ocean)then
+      rnat=0.
+      tslab(1,1)=tsurf(1)
+      tslab(1,2)=tsurf(1)
+      tsea_ice=tsurf
+      pctsrf_sic=0.
+      sea_ice=0.
+      endif
+
+
+
 c  Write a "startfi" file
 c  --------------------
@@ -729 +755 @@
      &                dtphys,time,
      &                tsurf,tsoil,emis,q2,qsurf,
-     &                cloudfrac,totcloudfrac,hice)
+     &                cloudfrac,totcloudfrac,hice, 
+     &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
 
 !      call physdem1(1,"startfi.nc",long,lati,nsoilmx,nq,
@@ -891 +918 @@
 #include "../dyn3d/disvert.F"
 #include "../dyn3d/abort_gcm.F"
+#include "../dyn3d/diverg.F"
+#include "../dyn3d/grad.F"
+#include "../dyn3d/gr_u_scal.F"
+#include "../dyn3d/gr_v_scal.F"
+#include "../dyn3d/gr_dyn_fi.F"
+