Changeset 2990 for trunk/LMDZ.MARS/util/expandstartfi.F90

Timestamp:

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

Author:

emillour

Message:

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

File:

• trunk/LMDZ.MARS/util/expandstartfi.F90 (modified) (1 diff)

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)