Changeset 2990

Timestamp:

Jul 6, 2023, 10:28:41 AM (18 months ago)

Author:

emillour

Message:

Mars PCM:
Minor fix in utility "expandstartfi": correctly handle the case of
polar mesh areas.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• util/expandstartfi.F90 (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4086 +4086 @@
 (ngrid,nlay) where conducted within big (ngrid,nlay) loop. Two occurrences :
 initialization of zq0 and preventing negative tracers. 
+
+== 06/07/2023 == EM
+Minor fix in utility "expandstartfi": correctly handle the case of
+polar mesh areas.

trunk/LMDZ.MARS/util/expandstartfi.F90

@@ -486 +486 @@
     
     ! "convert" from surf_field(:) to out_surf_field(:,:)
-    do i=1,lonlen
+    if (trim(varname)=="area") then 
+     ! Very specific handling of mesh area at the poles
+     ! which need to be split along all longitudes
+     do i=1,lonlen
+      out_surf_field(i,1)=surf_field(1)/(lonlen-1)
+      out_surf_field(i,latlen)=surf_field(physical_points)/(lonlen-1)
+     enddo    
+    else ! All other variables
+     do i=1,lonlen
       out_surf_field(i,1)=surf_field(1) ! North Pole
       out_surf_field(i,latlen)=surf_field(physical_points) ! South Pole
-    enddo
+     enddo
+    endif ! of if (trim(varname)=="area")
     do j=2,latlen-1
       ig0=1+(j-2)*(lonlen-1)