Re: use (/ .. /) if this is not the desired result

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah> Date: Thu, 30 Mar 2006 11:58:13 -0700 (MST) · >

---
So, you are reading from the output file? or did u mean
                                          grd=intide->tide_Ephase
If u do a                                     ^^^^^^
 printVarSummary(grd)
 
You will see that "grid" has meta data associated with it.
meta data: attributes and named dimensions.
values:    all 0.0
>
>
>; Fill the scratch array and write it back out to pre-defined variable
>; Cangle
>print("Cangle")
>nodevalues=intide->tide_Cangle
Again the variable "nodevalues" will (presumably) have 
meta data associated with it.
  printVarSummary(nodevalues)
>print("After nodevalues=")
>do iperiod = 0,tide_period-1
>  kperiod = periodpointer(iperiod)
>  Fperiod = tidefacpointer(iperiod)
>   do k = 0,node-1
>     grd(iperiod,i_rho(k),j_rho(k))=nodevalues(k,kperiod)
Here is where the error messages occur. 
The above is what NCL calls a "variable-to-variable" assignment.
[No operation is being performed on the right hand
 side of the = symbol.] 
 
When performing variable-to-variable  assignment NCL will 
transfer [copy] meta data from the rhs variable to the lhs
variable. When there are conflicts, NCL will use meta data 
associated with the rhs variable to overwrite the corresponding
meta data associated with the lhs variable.
http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml#VariableToVariable
      warning:VarVarWrite: 
      Dimension names for dimension number (0) don't match,
      assigning name of rhs dimension to lhs, 
      use "(/../)" if this change is not desired
      
      lhs - left hand side   [included so that non-native
      rhs - right hand side   English speakers know what
                              the abreviations mean]
I'll bet that the 0th [leftmost] dimension  name of
variable "grd" [lhs variable] does not match the 0th  
[leftmost] dimension name of "nodevalues" [rhs variable].
The warning message is informing you that it will
overwrite variable "grd" 0th dimension name with 
"nodevalues" 0th dimension name.
Same for the 1-th dimension name.
   grd(iperiod,i_rho(k),j_rho(k)) = (/ nodevalues(k,kperiod) /)
means that the values will be transferred but 
no meta data will be transferred. Thus, the "grd" variable
will retain its original meta data.
If you want "grd" to have the same dimension name as "nodevalues"
you can do the following after the loop
   grd!0 = nodevalues!0    ; the leftmost dimension 
                           ; is the 0th dimension
Note: I am attaching a code snippet that creates a netCDF file.
      It uses the "setfileoption" procedure [available
      in v033] to set "DefineMode" for netCDF.
      
      This can *substantially* reduce the time
      it takes to write netCDF files.
      
Regards
D