Changeset 2294

Timestamp:

Apr 22, 2020, 3:39:42 PM (5 years ago)

Author:

adelavois

Message:

MARS GCM:
Minor update of revision 2281 about initialisation of emissivity and albedo when no startfi
Flags to change these fields are now more explicit: surfemis_without_startfi instead of surfemis (idem for surfalbedo)

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/phyetat0_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -2959 +2959 @@
 == 20/04/2020 == EM
 Fix for the 1D, follow up from the recent addition of watercap.
+
+== 22/04/2020 == AD
+Minor update of revision 2281 about initialisation of emissivity and albedo when no startfi
+Now default emissivity/albedo is 0.95/0.1
+Flags to change these fields are now more explicit: surfemis_without_startfi instead of surfemis (idem for surfalbedo)

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

@@ -107 +107 @@
 write(*,*) "phyetat0: startphy_file", startphy_file
 
-if(.not.startphy_file) then
-   surfemis = -999
-   surfalbedo = -999
-   call getin_p("surfemis",surfemis)
-   call getin_p("surfalbedo",surfalbedo)
-   if (surfemis.eq.-999) then
-      call abort_physic(modname, &
-               "Missing value for surfemis in def files!",1)
-   endif
-   if (surfalbedo.eq.-999) then
-      call abort_physic(modname, &
-               "Missing value for surfalbedo in def files!",1)
-   endif
-   write(*,*) "phyetat0: surfemis: ", surfemis
-   write(*,*) "phyetat0: surfalbedo: ", surfalbedo
-endif
-
 if (startphy_file) then
    ! open physics initial state file:
@@ -159 +142 @@
    endif
 else ! If no startfi file, use parameter surfalbedo in def file
-  surfalbedo=0.5
-  call getin_p("surfalbedo",surfalbedo)
-  print*,"surfalbedo",surfalbedo
+  surfalbedo=0.1
+  call getin_p("surfalbedo_without_startfi",surfalbedo)
+  print*,"surfalbedo_without_startfi",surfalbedo
   albedodat(:)=surfalbedo
 endif ! of if (startphy_file)
@@ -508 +491 @@
 else
   ! If no startfi file, use parameter surfemis in def file
-  surfemis=1.0
-  call getin_p("surfemis",surfemis)
-  print*,"surfemis",surfemis
+  surfemis=0.95
+  call getin_p("surfemis_without_startfi",surfemis)
+  print*,"surfemis_without_startfi",surfemis
   emis(:)=surfemis
 endif ! of if (startphy_file)