Changeset 5157 for TOOLS

Timestamp:

09/18/20 17:02:09 (4 years ago)

Author:

omamce

Message:

O.M. CreateRestartAtm4Oasis.bash correct problem on dimensions in nco script

File:

• TOOLS/CPLRESTART/CreateRestartAtm4Oasis.bash (modified) (3 diffs)

TOOLS/CPLRESTART/CreateRestartAtm4Oasis.bash

@@ -112 +112 @@
     IsUnstructured=yes
     dim_cell=$(ncdump -h flxat_fields_notime.nc |& grep "cell *=" | awk '{print $3}' ) # Get nominal DYNAMICO resolution (nbp parameter)
-    ico_nbp=$( echo "sqrt(($dim_cell-2)/10)+1" | bc -l | sed 's/\..*//' )
+    ico_nbp=$( echo "sqrt(($dim_cell-2)/10)" | bc -l | sed 's/\..*//' )
     ATM=ICO${ico_nbp}
     ncrename --history --dimension cell,y flxat_fields_notime.nc
@@ -188 +188 @@
 COTOTSNO[y,x] = double(snow)           ;  //  Solid precipitation
 COTOTEVA[y,x] = double(evap_oce*fract_oce_norm + evap_sic*fract_sic_norm) ;
-COICEVAP[y,x] = double(evap_sic(:,:))  ; // Evaporation on sea ice
-COQSRMIX[y,x] = double(sols(:,:))      ; // Heat flux short wave
-COQNSMIX[y,x] = double(soll(:,:))      ; // Heat flux minus short wave
-COSHFICE[y,x] = double(sols(:,:))      ; // Heat flux short wave over sea ice
-CONSFICE[y,x] = double(soll(:,:))      ; // Heat flux minus short wave over sea ice
+COICEVAP[y,x] = double(evap_sic)       ; // Evaporation on sea ice
+COQSRMIX[y,x] = double(sols)           ; // Heat flux short wave
+COQNSMIX[y,x] = double(soll)           ; // Heat flux minus short wave
+COSHFICE[y,x] = double(sols)           ; // Heat flux short wave over sea ice
+CONSFICE[y,x] = double(soll)           ; // Heat flux minus short wave over sea ice
 CODFLXDT[y,x] = -20.0d       ; // W/m2 - dQ/dt
 COCALVIN[y,x] =   0.0d       ; // Calving of icebergs, solid
 COLIQRUN[y,x] =   0.0d       ; // River run-off , liquid
-COWINDSP[y,x] = double( wind10m(:,:) ) ; // Wind speed at 10m high
+COWINDSP[y,x] = double( wind10m )      ; // Wind speed at 10m high
 COTAUMOD[y,x] =   0.0d                 ; // Wind stress modulus
 // Wind stress
@@ -204 +204 @@
 // Wind stress, from east/north components to geocentric
 *rpi = acos (-1.0d) ;  *rad = rpi / 180.0d  ; 
-COTAUXXU[y,x] = double( -tau_x(:,:) * sin(rad * lon(:,:)) - tau_y(:,:) * sin(rad * lat(:,:)) * cos(rad * lon(:,:)) ); 
-COTAUYYU[y,x] = double(  tau_x(:,:) * cos(rad * lon(:,:)) - tau_y(:,:) * sin(rad * lat(:,:)) * sin(rad * lon(:,:)) );
-COTAUZZU[y,x] = double(  tau_y(:,:) * cos(rad * lat(:,:)) ) ;
+COTAUXXU[y,x] = double( -tau_x * sin(rad * lon) - tau_y * sin(rad * lat) * cos(rad * lon) ); 
+COTAUYYU[y,x] = double(  tau_x * cos(rad * lon) - tau_y * sin(rad * lat) * sin(rad * lon) );
+COTAUZZU[y,x] = double(  tau_y * cos(rad * lat) ) ;
 // Value at North Pole
 #include "lonlat2xyz.nco"
 //
-COTAUXXV[y,x] = COTAUXXU(:,:) ;
-COTAUYYV[y,x] = COTAUYYU(:,:) ;
-COTAUZZV[y,x] = COTAUZZU(:,:) ;
+COTAUXXV[y,x] = COTAUXXU ;
+COTAUYYV[y,x] = COTAUYYU ;
+COTAUZZV[y,x] = COTAUZZU ;
 //
 EOF