Changeset 13

Timestamp:

Jan 24, 2011, 2:02:24 PM (14 years ago)

Author:

aslmd

Message:

spiga@svn-planeto:mineur

File:

• trunk/mesoscale/LMD_MM_MARS/SRC/WPS/geogrid/src/process_tile_module.F90 (modified) (7 diffs)

trunk/mesoscale/LMD_MM_MARS/SRC/WPS/geogrid/src/process_tile_module.F90

@@ -149 +149 @@
       !
       call mprintf(.true.,STDOUT,'  Processing XLAT and XLONG')
-    
+   
       if (grid_type == 'C') then
          call get_lat_lon_fields(which_domain, xlat_array, xlon_array, start_mem_i, &
@@ -162 +162 @@
          corner_lats(3) = xlat_array(end_patch_i,end_patch_j)
          corner_lats(4) = xlat_array(end_patch_i,start_patch_j)
-     
+
          corner_lats(5) = xlat_array_u(start_patch_i,start_patch_j)
          corner_lats(6) = xlat_array_u(start_patch_i,end_patch_j)
@@ -314 +314 @@
       if (grid_type == 'C') then
          call mprintf(.true.,STDOUT,'  Processing ROTANG')
-    
+   
          allocate(sina_array(start_mem_i:end_mem_i, start_mem_j:end_mem_j))
          allocate(cosa_array(start_mem_i:end_mem_i, start_mem_j:end_mem_j))
@@ -525 +525 @@
             call write_field(start_mem_i, end_mem_i, start_mem_j, end_mem_j, 1, 1, trim(domname), &
                              datestr, dominant_field)
-          
+         
             deallocate(dominant_field)
          end if
@@ -633 +633 @@
                               if (field(i,j,k) /= msg_fill_val) then
                                  field(i,j,k) = field(i,j,k) * scale_factor
+!!!! debug: to display the interpolated values 
+!                                   call mprintf(.true.,STDOUT,'%i',i1=int(field(i,j,k)))
                               end if
                            end do
@@ -678 +680 @@
 
                   end if
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!BEURK
+!!!!!!!!!!!!!!!!!!!!!!!!!!
+!if (fieldname == 'HGT_M') then
+!!          open(unit=100,file='./custom_topo/hm12',form='formatted',status='old')
+!!          open(unit=100,file='./custom_topo/hm6',form='formatted',status='old')
+!          open(unit=100,file='./custom_topo/hm18',form='formatted',status='old')
+!          rewind(100)
+!          do j=sm2+3,em2-3  !!ok
+!          do i=sm1+3,em1-3  !!ok
+!                read(100,*) field(i,j,1)
+!          enddo
+!          enddo
+!          close(100)    
+!endif
+!if (fieldname == 'HGT_U') then
+!!          open(unit=101,file='./custom_topo/hu12',form='formatted',status='old')
+!!          open(unit=101,file='./custom_topo/hu6',form='formatted',status='old')
+!          open(unit=101,file='./custom_topo/hu18',form='formatted',status='old')
+!          rewind(101)
+!          do j=sm2+3,em2-3
+!          do i=sm1+3,em1-3!4
+!                read(101,*) field(i,j,1)
+!          enddo
+!          enddo
+!          close(101)
+!endif
+!if (fieldname == 'HGT_V') then
+!!          open(unit=102,file='./custom_topo/hv12',form='formatted',status='old')
+!!          open(unit=102,file='./custom_topo/hv6',form='formatted',status='old')
+!          open(unit=102,file='./custom_topo/hv18',form='formatted',status='old')
+!          rewind(102)
+!          do j=sm2+3,em2-3!4
+!          do i=sm1+3,em1-3
+!                read(102,*) field(i,j,1)
+!          enddo
+!         enddo
+!          close(102)
+!endif
+!!!!!!!!!!!!!!!!!!!!!!!!!!
+!!!BEURK
+!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+
+
 
                   call write_field(sm1, em1, sm2, em2, &
@@ -1293 +1341 @@
       end if
 
+!**** MARS STUFF
+! the 9000. translation was done in the Mars GEOG file to avoid negative numbers
+        if (fieldname == 'HGT_M' .or. fieldname == 'HGT_U' .or. fieldname == 'HGT_V') then 
+
+                field(:,:,:)=field(:,:,:)-9000.         
+        
+        !print *,'OK'           
+        end if
+!**** MARS STUFF
+
+
       deallocate(interp_type)