Changeset 3718

Timestamp:

03/29/18 10:56:01 (7 years ago)

Author:

omamce

Message:

O.M. : add computation of runoff weights

Location:

TOOLS/MOSAIX

Files:

• CreateOasisGrids.bash (modified) (2 diffs)
• CreateWeightsMask.bash (modified) (8 diffs)
• src/MOSAIX/interpol.f90 (modified) (2 diffs)

TOOLS/MOSAIX/CreateOasisGrids.bash

@@ -361 +361 @@
             --attribute LOGNAME,global,o,c,"$(whoami)"                 \
             --attribute Program,global,o,c,"Generated by ${0}"         \
-            --attribute name,global,o,c,"${type}_${CplModel}.nc"     \
+            --attribute name,global,o,c,"${type}_${CplModel}.nc"       \
             --attribute directory,global,o,c,"$(pwd)"                  \
             --attribute description,global,o,c,"Fields needed by OASIS-MCT" \
-            --attribute title,global,o,c,"${type}_${CplModel}.nc"    \
+            --attribute title,global,o,c,"${type}_${CplModel}.nc"      \
             --attribute timeStamp,global,o,c,"$(date)"                 \
             --attribute associatedFiles,global,o,c,"grids_${CplModel}.nc areas_${CplModel}.nc masks_${CplModel}.nc" \
@@ -373 +373 @@
             --attribute hardware,global,o,c,"$(uname -i)"              \
             --attribute Comment,global,o,c,"Preliminary attempt - Do not trust !" \
-            --attribute SVN_Author,global,o,c,"$Author$"                           \
-            --attribute SVN_Date,global,o,c,"$Date$"                            \
+            --attribute SVN_Author,global,o,c,"$Author$"      \
+            --attribute SVN_Date,global,o,c,"$Date$" \
             --attribute SVN_Revision,global,o,c,"$Revision$"                    \
-            --attribute SVN_Id,global,o,c,"$Id$"                                \
+            --attribute SVN_Id,global,o,c,"$Id$" \
+            --attribute SVN_HeardURL,global,o,c,"$HeadURL $" \
             ${type}_${CplModel}.nc
 done

TOOLS/MOSAIX/CreateWeightsMask.bash

@@ -48 +48 @@
 #OCE=ORCA025
 
-ATM=ICO30
+#ATM=ICO30
 #ATM=ICO40
 #ATM=ICO450
-#ATM=LMD9695
+ATM=LMD9695
 #ATM=LMD144X142
 
@@ -180 +180 @@
 ncks --alphabetize --history --append --variable OceMask dia_t${oce}_to_t${atm}_1storder_mask.nc   ${ATM}_grid_maskFrom_${OCE}.nc
 
+##
+## Create mask of coastal ATM point
+## ===========================================================================
+cat <<EOF > coastal.nco
+AtmCoastal = OceFrac * 0.0 ;
+where (OceFrac > 0.0 && OceFrac < 1.0 )  AtmCoastal = 1.0 ;
+EOF
+ncap2 --history --overwrite --script-file coastal.nco ${ATM}_grid_maskFrom_${OCE}.nc ${ATM}_coastal_maskFrom_${OCE}.nc
+ncks --history --append --variable AtmCoastal ${ATM}_coastal_maskFrom_${OCE}.nc ${ATM}_grid_maskFrom_${OCE}.nc
+rm ${ATM}_coastal_maskFrom_${OCE}.nc
 
 ##
@@ -255 +265 @@
 
 ##
+## ATM towards NEMO points - runoff
+## ===========================================================================
+cp ${SUBMIT_DIR}/iodef_atm_to_oce.xml           iodef.xml
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_read"]/domain_definition/domain[@id="domain_src"]'                   -k type  -v ${atm_domain_type}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_read"]/domain_definition/domain[@id="domain_dst"]'                   -k type  -v ${oce_domain_type}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/file_definition/file[@id="dia"]/variable[@name="title"]'       -t "${ATM} coastal mask interpolated to ${OCE}"
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/file_definition/file[@id="dia"]/variable[@name="source_grid"]' -t ${atm_domain_type}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/file_definition/file[@id="dia"]/variable[@name="dest_grid"]'   -t ${oce_domain_type}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/domain_definition/domain[@id="domain_src"]'                    -k type  -v ${atm_domain_type}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/domain_definition/domain[@id="domain_dst"]'                    -k type  -v ${oce_domain_type}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/domain_definition/domain[@id="domain_dst"]/interpolate_domain' -k order -v 1
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/file_definition/file[@id="dia"]/variable[@name="order"]'       -t 1
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_read"]/file_definition/file[@id="file_src"]'                         -k name  -v ${ATM}_grid_maskFrom_${OCE}
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_read"]/file_definition/file[@id="file_src"]/field[@id="mask_src"]'   -k name  -v AtmCoastal
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/file_definition/file[@id="dia"]'                               -k name  -v dia_c${atm}_to_t${oce}_1storder
+python update_xml.py -i iodef.xml -n 'context[@id="interpol_run"]/domain_definition/domain[@id="domain_dst"]/interpolate_domain' -k weight_filename -v rmp_c${atm}_to_t${oce}_1storder.nc
+
+cp iodef.xml iodef_c${atm}_to_t${oce}_1storder.xml
+
+time /usr/bin/time ccc_mprun ./interpol.exe --mask_src=yes --mask_dst=yes
+
+##
 ## Change all NetCDF files to NetCDF 3 format (needed for OASIS)
 ## ===========================================================================
 for InFile in *.nc ; do
-    if [[ $(ncdump -k ${InFile}) = *netCDF-4* ]] ; then
-        mv ${InFile} tmp_${InFile}
-        ncks --alphabetize --fl_fmt=${FL_FMT} --history  tmp_${InFile}  ${InFile}
-    fi
+    mv ${InFile} tmp_${InFile}
+    ncks --alphabetize --fl_fmt=${FL_FMT} --history  tmp_${InFile}  ${InFile}
 done
 
@@ -385 +415 @@
 ## ===========================================================================
 for OGRID in ${ListOCEGRID} ; do
-    for order in ${Listorder} ; do
-        ogrid=${OGRID~}
-        agrid=t
+    for AGRID in T C ; do
+        for order in ${Listorder} ; do
+            ogrid=${OGRID~}
+            agrid=${AGRID~}
         
-        file_o_to_a=rmp_${ogrid}${oce}_to_${agrid}${atm}_${order}order.nc
-        file_a_to_o=rmp_${agrid}${atm}_to_${ogrid}${oce}_${order}order.nc
-        
-        cat <<EOF >add_varoce.nco
+            file_o_to_a=rmp_${ogrid}${oce}_to_${agrid}${atm}_${order}order.nc
+            file_a_to_o=rmp_${agrid}${atm}_to_${ogrid}${oce}_${order}order.nc
+            
+            cat <<EOF >add_varoce.nco
 defdim ("src_grid_size"   , \$x_grid_${OGRID}.size*\$y_grid_${OGRID}.size) ;
 defdim ("src_grid_corners", 4) ;
@@ -417 +448 @@
 EOF
 
-        ncap2 --history --append --script-file add_varoce.nco oce_grid.nc ${file_o_to_a}
-        sed --in-place "s/src_/dst_/g" add_varoce.nco
-        ncap2 --history --append --script-file add_varoce.nco oce_grid.nc ${file_a_to_o}
-
-        if [[ ${atm} = ico ]] ; then
-            cat <<EOF >add_varatm.nco
+            [[ -f ${file_o_to_a} ]] && ncap2 --history --append --script-file add_varoce.nco oce_grid.nc ${file_o_to_a}
+            sed --in-place "s/src_/dst_/g" add_varoce.nco
+            [[ -f ${file_a_to_o} ]] && ncap2 --history --append --script-file add_varoce.nco oce_grid.nc ${file_a_to_o}
+            
+            if [[ ${atm} = ico ]] ; then
+                cat <<EOF >add_varatm.nco
 defdim ("dst_grid_size"   , \$cell.size) ;
 defdim ("dst_grid_corners", 6) ;
@@ -446 +477 @@
 dst_grid_area  (:) = aire(:) ;
 EOF
-            ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_o_to_a}
-            sed --in-place "s/dst_/src_/g" add_varatm.nco
-            ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_a_to_o}
-            
-            ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_corner_lon,src_grid_corner_lat,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_corner_lon,dst_grid_corner_lat,dst_grid_area,dst_grid_imask ${file_o_to_a} rmp_tmp.nc
-            mv rmp_tmp.nc ${file_o_to_a}
-            ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_corner_lon,src_grid_corner_lat,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_corner_lon,dst_grid_corner_lat,dst_grid_area,dst_grid_imask ${file_a_to_o} rmp_tmp.nc
-            mv rmp_tmp.nc ${file_a_to_o}
-
-        fi
-        if [[ ${atm} = lmd ]] ; then
-            cat <<EOF >add_varatm.nco
+                if [[ -f ${file_o_to_a} ]] ; then
+                    ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_o_to_a}
+                    ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_corner_lon,src_grid_corner_lat,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_corner_lon,dst_grid_corner_lat,dst_grid_area,dst_grid_imask ${file_o_to_a} rmp_tmp.nc
+                    mv rmp_tmp.nc ${file_o_to_a}
+                fi
+                    
+                sed --in-place "s/dst_/src_/g" add_varatm.nco
+                if [[ -f ${file_a_to_o} ]] ; then
+                    ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_a_to_o}
+                    ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_corner_lon,src_grid_corner_lat,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_corner_lon,dst_grid_corner_lat,dst_grid_area,dst_grid_imask ${file_a_to_o} rmp_tmp.nc
+                    mv rmp_tmp.nc ${file_a_to_o}
+                fi
+            fi
+            if [[ ${atm} = lmd ]] ; then
+                cat <<EOF >add_varatm.nco
 defdim ("dst_grid_size"   , \$lon.size*\$lat.size) ;
 defdim ("dst_grid_corners", 4) ;
@@ -484 +518 @@
 dst_grid_frac (:)  = OceFrac(0,:,:) ;
 EOF
-            ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_o_to_a}
-            sed --in-place "s/dst/src/g" add_varatm.nco
-            ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_a_to_o}
-
-            ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_corner_lon,src_grid_corner_lat,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_area,dst_grid_imask ${file_o_to_a} rmp_tmp.nc
-            mv rmp_tmp.nc ${file_o_to_a}
-            ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_corner_lon,dst_grid_corner_lat,dst_grid_area,dst_grid_imask ${file_a_to_o} rmp_tmp.nc
-            mv rmp_tmp.nc ${file_a_to_o}
-
-        fi
+                if [[ -f ${file_o_to_a} ]] ; then
+                    ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_o_to_a}
+                    ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_corner_lon,src_grid_corner_lat,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_area,dst_grid_imask ${file_o_to_a} rmp_tmp.nc
+                    mv rmp_tmp.nc ${file_o_to_a}
+                fi
+                
+                sed --in-place "s/dst/src/g" add_varatm.nco
+                if [[ -f ${file_a_to_o} ]] ; then
+                    ncap2 --history --append --script-file add_varatm.nco ${ATM}_grid_maskFrom_${OCE}.nc ${file_a_to_o}
+                    ncks --alphabetize --history --overwrite --variable src_address,dst_address,remap_matrix,src_grid_dims,src_grid_center_lat,src_grid_center_lon,src_grid_area,src_grid_imask,dst_grid_dims,dst_grid_center_lat,dst_grid_center_lon,dst_grid_corner_lon,dst_grid_corner_lat,dst_grid_area,dst_grid_imask ${file_a_to_o} rmp_tmp.nc
+                    mv rmp_tmp.nc ${file_a_to_o}
+                fi
+            fi
+        done
     done
 done
@@ -501 +539 @@
 ## ===========================================================================
 cp ${ATM}_grid_maskFrom_${OCE}.nc ${SUBMIT_DIR}
+
+cp rmp_c${atm}_to_t${oce}_${order}order.nc ${SUBMIT_DIR}/rmp_c${ATM}_to_t${OCE}_1storder.nc
+cp dia_c${atm}_to_t${oce}_${order}order.nc ${SUBMIT_DIR}/dia_c${ATM}_to_t${OCE}_1storder.nc
+
 for order in ${Listorder}
 do
     cp dia_t${oce}_to_t${atm}_${order}order.nc ${SUBMIT_DIR}/dia_t${OCE}_to_t${ATM}_${order}order.nc
     cp rmp_t${oce}_to_t${atm}_${order}order.nc ${SUBMIT_DIR}/rmp_t${OCE}_to_t${ATM}_${order}order.nc
+  
     for grid in ${Listocegrid} ; do
         cp rmp_t${atm}_to_${grid}${oce}_${order}order.nc ${SUBMIT_DIR}/rmp_t${ATM}_to_${grid}${OCE}_${order}order.nc

TOOLS/MOSAIX/src/MOSAIX/interpol.f90

@@ -111 +111 @@
       CALL xios_recv_field ("mask_src", imask_src)
       lmask_src = .FALSE.
-      WHERE (imask_src>0) lmask_src = .TRUE.
+      WHERE (imask_src > 0.5 ) lmask_src = .TRUE.
    ELSE
       imask_src (:,:) = 1 ; lmask_src (:,:) = .TRUE.
    ENDIF
+   WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- counting mask_src : ', 1I8)" ) rank, COUNT(lmask_src)
    
    !< Read mask on the destination grid
@@ -125 +126 @@
       CALL xios_recv_field ("mask_dst", imask_dst)
       lmask_dst = .FALSE.
-      WHERE (imask_dst>0) lmask_dst = .TRUE.
+      WHERE (imask_dst > 0.5) lmask_dst = .TRUE.
    ELSE
       imask_dst (:,:) = 1 ; lmask_dst (:,:) = .TRUE.
    ENDIF
+   WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- counting mask_dst : ', 1I8)" ) rank, COUNT(lmask_dst)
 
    WRITE (UNIT=nout, FMT="('-- ', 1I4.4, ' -- fin interpol_read (1)')") rank