Changeset 801 for trunk/LMDZ.MARS/libf

Timestamp:

Oct 5, 2012, 2:25:36 AM (12 years ago)

Author:

aslmd

Message:

LMDZ.MARS adapted icap trick to mesoscale modeling

File:

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

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

@@ -13 +13 @@
 #include"surfdat.h"
 #include"callkeys.h"
+#ifdef MESOSCALE
+#include"comgeomfi.h"
+#endif
 
 ! arguments:
@@ -51 +54 @@
 
 do ig=1,ngrid
+#ifndef MESOSCALE
   if (ig.gt.ngrid/2+1) then
+#else
+  if (lati(ig)*180./acos(-1.).lt.0.) then
+#endif
     icap=2 ! Southern hemisphere
   else
@@ -64 +71 @@
 !      write(*,*) "albedocaps: call TES_icecap_albedo"
 !      write(*,*) "albedocaps: zls=",zls," ig=",ig
-      call TES_icecap_albedo(zls,ig,psolaralb(ig,1))
+      call TES_icecap_albedo(zls,ig,psolaralb(ig,1),icap)
 !      write(*,*) "albedocaps: psolaralb(ig,1)=",psolaralb(ig,1)
       psolaralb(ig,2)=psolaralb(ig,1)
@@ -89 +96 @@
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-subroutine TES_icecap_albedo(zls,ig,alb)
+subroutine TES_icecap_albedo(zls,ig,alb,icap)
 
 implicit none
@@ -103 +110 @@
 integer,intent(in) :: ig ! grid point index
 real,intent(out) :: alb ! (interpolated) TES ice albedo at that grid point
+integer :: icap ! =1: Northern hemisphere =2: Southern hemisphere
 
 ! local variables:
@@ -118 +126 @@
 real :: lond ! longitude, in degrees
 integer :: i
-integer :: icap ! =1: Northern hemisphere =2: Southern hemisphere
 
 ! TES datasets: (hard coded fixed length/sizes; for now)
@@ -231 +238 @@
 endif ! of if firstcall
 
-! 1. Identify hemisphere and encompassing latitudes
-if (ig.GT.ngridmx/2+1) then
-  icap=2 ! Southern hemisphere
-else
-  icap=1 ! Northern hemisphere
-endif
+! 1. Identify encompassing latitudes
 
 ! Check that latitude is such that there is TES data to use