Changeset 2239 for LMDZ5/trunk/libf/dynlmdz_phylmd/iniphysiq.F90

Timestamp:

Mar 23, 2015, 8:27:30 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Reorganizing physics/dynamics interface:

• what is related to dynamics-physics interface is now in a seperate directory: dynlmdz_phy* for physics in phy*
• 1d model and related dependencies (including a couple from "dynamics", set up as symbolic links) is now in subdirectory "dyn1d" of phy*.
• "bibio" directory is now "misc" and should only contain autonomous utilities.
• "cosp" is now a subdirectory of phylmd.

EM

Location:

LMDZ5/trunk/libf/dynlmdz_phylmd

Files:

• . (added)
• iniphysiq.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/iniphysiq.F90) (5 diffs)

LMDZ5/trunk/libf/dynlmdz_phylmd/iniphysiq.F90

@@ -148 +148 @@
     WRITE (lunout, *) '  in the dynamics punjours=', punjours
     WRITE (lunout, *) '   but in the physics RDAY=', rday
-    IF (abs(rday-punjours)>0.01) THEN
+    IF (abs(rday-punjours)>0.01*punjours) THEN
         ! stop here if the relative difference is more than 1%
       abort_message = 'length of day discrepancy'
@@ -158 +158 @@
     WRITE (lunout, *) '     in the dynamics pg=', pg
     WRITE (lunout, *) '  but in the physics RG=', rg
-    IF (abs(rg-pg)>0.01) THEN
+    IF (abs(rg-pg)>0.01*pg) THEN
         ! stop here if the relative difference is more than 1%
       abort_message = 'gravity discrepancy'
@@ -168 +168 @@
     WRITE (lunout, *) '   in the dynamics prad=', prad
     WRITE (lunout, *) '  but in the physics RA=', ra
-    IF (abs(ra-prad)>0.01) THEN
+    IF (abs(ra-prad)>0.01*prad) THEN
         ! stop here if the relative difference is more than 1%
       abort_message = 'planet radius discrepancy'
@@ -178 +178 @@
     WRITE (lunout, *) '     in the dynamics pr=', pr
     WRITE (lunout, *) '  but in the physics RD=', rd
-    IF (abs(rd-pr)>0.01) THEN
+    IF (abs(rd-pr)>0.01*pr) THEN
         ! stop here if the relative difference is more than 1%
       abort_message = 'reduced gas constant discrepancy'
@@ -188 +188 @@
     WRITE (lunout, *) '     in the dynamics pcpp=', pcpp
     WRITE (lunout, *) '  but in the physics RCPD=', rcpd
-    IF (abs(rcpd-pcpp)>0.01) THEN
+    IF (abs(rcpd-pcpp)>0.01*pcpp) THEN
         ! stop here if the relative difference is more than 1%
       abort_message = 'specific heat discrepancy'