Changeset 2294 for trunk/LMDZ.MARS/libf/phymars

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)

File:

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

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)