Changeset 707 for trunk/LMDZ.MARS/libf

Timestamp:

Jun 22, 2012, 1:06:09 PM (13 years ago)

Author:

emillour

Message:

Mars GCM:

• Enforce that imposed TES ice albedo can never be greater than 0.9

EM

File:

• trunk/LMDZ.MARS/libf/phymars/albedocaps.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/albedocaps.F90

@@ -159 +159 @@
     write(*,*)'Problem opening npsc_albedo.nc (phymars/albedocaps.F90)'
     write(*,*)'It should be in :',trim(datafile),'/'
-    write(*,*)'1) You can change this directory address in '
-    write(*,*)'   file phymars/datafile.h'
+    write(*,*)'1) You can change this directory address in callfis.def with'
+    write(*,*)'   datadir=/path/to/datafiles'
     write(*,*)'2) If necessary, npsc_albedo.nc (and other datafiles)'
     write(*,*)'   can be obtained online on:'
@@ -200 +200 @@
     write(*,*)'Problem opening spsc_albedo.nc (phymars/albedocaps.F90)'
     write(*,*)'It should be in :',trim(datafile),'/'
-    write(*,*)'1) You can change this directory address in '
-    write(*,*)'   file phymars/datafile.h'
+    write(*,*)'1) You can change this directory address in callfis.def with'
+    write(*,*)'   datadir=/path/to/datafiles'
     write(*,*)'2) If necessary, spsc_albedo.nc (and other datafiles)'
     write(*,*)'   can be obtained online on:'
@@ -379 +379 @@
 endif ! of if (icap.eq.1)
 
+! Make sure that returned albedo is never greater than 0.90
+if (alb.gt.0.90) alb=0.90
+
 end subroutine TES_icecap_albedo