Changeset 1150 for trunk/LMDZ.GENERIC

Timestamp:

Jan 6, 2014, 10:27:56 PM (12 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. put surface temperature update under .not.nosurf flag. commented and added prints for ichoice=8. added a modification of iradia in rcm1d.def to be able to still use same .def as subsequent 3D runs for consistency

Location:

trunk/LMDZ.GENERIC

Files:

• deftank/saturn1d/compile (modified) (1 diff)
• deftank/saturn1d/profile (added)
• deftank/saturn1d/rcm1d.def (modified) (2 diffs)
• libf/phystd/physiq.F90 (modified) (2 diffs)
• libf/phystd/profile.F (modified) (1 diff)
• libf/phystd/rcm1d.F (modified) (1 diff)

trunk/LMDZ.GENERIC/deftank/saturn1d/compile

@@ -1 +1 @@
 #! /bin/bash
 
-#/donnees/aslmd/MODELES/LMDZ.GENERIC/makegcm -d 32 -t 0 -b 38x36 rcm1d
 /donnees/aslmd/MODELES/LMDZ.GENERIC/makegcm -d 64 -t 0 -b 20x30 rcm1d
+#/donnees/aslmd/MODELES/LMDZ.GENERIC/makegcm -d 32 -t 0 -b 20x30 rcm1d

trunk/LMDZ.GENERIC/deftank/saturn1d/rcm1d.def

@@ -41 +41 @@
 # Number of steps between each writing in diagfi
 ecritphy=5000
+
+# OVERRIDE callphys.def settings
+# ... this to ensure compliancy with the low day_step value above
+iradia = 25
 
 
@@ -89 +93 @@
 #         ichoice=7   constant T  + gaussian perturbation (km)
 #         ichoice=8   Read in an ascii file "profile" 
-ichoice        = 1
+ichoice        = 8
 # Reference temperature tref (K)
 tref           = 200.

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -1350 +1350 @@
 !-----------------------------------------------------------------------
 
+      if (.not.nosurf) then
+
 !     Increment surface temperature
       tsurf(1:ngrid)=tsurf(1:ngrid)+ptimestep*zdtsurf(1:ngrid) 
@@ -1358 +1360 @@
                 ptimestep,tsurf,tsoil,capcal,fluxgrd)
       endif
+
+      endif ! if (.not.nosurf)
 
 !-------------------------

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

@@ -175 +175 @@
 c   ---------
 
+      ! first value is surface temperature
+      ! then profile of atmospheric temperature
       ELSE IF(ichoice.GE.8) THEN
       OPEN(11,file='profile',status='old',form='formatted',err=101)

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

@@ -677 +677 @@
         temp(ilayer)=tmp2(ilayer)
       ENDDO
-      
+      print*,"check"
+      PRINT*,"INPUT SURFACE TEMPERATURE",tsurf
+      PRINT*,"INPUT TEMPERATURE PROFILE",temp
 
 c  Initialisation albedo / inertie du sol