if (err .ne. 0)
+ call errore('nf_open: ', err)
canConvert = (var_type(i) .eq. NF_CHAR) .eqv.
+ (NFT_TEXT .eq. NFT_TEXT)
if (err .ne. NF_EBADID)
+ call errore('bad ncid: ', err)
err = nf_get_var1_text(ncid, BAD_VARID,
+ index, value)
if (err .ne. NF_ENOTVAR)
+ call errore('bad var id: ', err)
do 3, j = 1, var_rank(i)
index(j) = var_shape(j,i) + 1
if (err .ne. NF_ECHAR)
+ call errore('conversion: ', err)
if (err .ne. NF_EINVALCOORDS)
+ call errore('bad index: ', err)
common /dims/timedim, latdim, londim, leveldim, lendim,
+ dimsiz
ncid = nccre(name, NCCLOB, iret)
if (ncid .eq. -1) then nfails = nfails + 1
incdf = ncopn(iname, NCNOWRIT, iret)
if (iret .ne. 0) nfails = nfails + 1
outcdf = nccre(oname, NCCLOB, iret)
if (iret .ne. 0) nfails = nfails + 1
call ncinq (incdf, ndims, nvars, natts, recdim, iret)
if (iret .ne. 0) nfails = nfails + 1
call ncanam (incdf, NCGLOBAL, j, attnam, iret)
if (iret .ne. 0) nfails = nfails + 1
call ncacpy (incdf, NCGLOBAL, attnam, outcdf, NCGLOBAL, iret)
call ncvinq (incdf, i, vname, vartyp, nvdims,
+ vdims, nvatts, iret)
if (iret .ne. 0) nfails = nfails + 1
if (err .ne. 0)
+ call errore('nf_open: ', err)
canConvert = (var_type(i) .eq. NF_CHAR) .eqv.
+ (NFT_ITYPE($1) .eq. NFT_TEXT)
if (err .ne. NF_EBADID)
+ call errore('bad ncid: ', err)
err = nf_get_var1_$1(ncid, BAD_VARID,
+ index, value)
if (err .ne. NF_ENOTVAR)
+ call errore('bad var id: ', err)
do 3, j = 1, var_rank(i)
index(j) = var_shape(j,i) + 1
if (err .ne. NF_ECHAR)
+ call errore('conversion: ', err)
if (err .ne. NF_EINVALCOORDS)
+ call errore('bad index: ', err)
hash_text = max(minimum, min(maximum, hash4( type, rank,
+ index, itype)))
hash_int1 = max(minimum, min(maximum, hash4( type, rank,
+ index, itype)))
hash_double = max(minimum, min(maximum, hash4( type, rank,
+ index, itype)))
if (err .ne. 0)
+ call errore('nf_open: ', err)
canConvert = (var_type(i) .eq. NF_CHAR) .eqv.
+ (NFT_TEXT .eq. NFT_TEXT)
err = nf_inq_var(ncid, i, name, datatype, ndims, dimids,
+ ngatts)
if (err .ne. 0)
if (name .ne. var_name(i))
+ call error('Unexpected var_name')
if (datatype .ne. var_type(i))
if (name .ne. var_name(i))
+ call error('Unexpected var_name')
if (datatype .ne. var_type(i))
+ call error('Unexpected type')
if (ndims .ne. var_rank(i))
+ call error('Unexpected rank')
if (err .ne. 0)
+ call errore('nf_inq_dim: ', err)
if (length .ne. var_shape(j,i))
+ call error('Unexpected shape')
err = index2indexes(j, var_rank(i), var_shape(1,i),
+ index)
if (err .ne. 0)
+ call error('error in index2indexes()')
expect = hash4( var_type(i), var_rank(i), index,
data msg(13)
+ / 'Not a netCDF data type or _FillValue type mismatch' /
data msg(24) / 'Attempt to convert between text & numbers' /
data msg(25) / 'Start+count exceeds dimension bound' /
data msg(27)
+ / 'Attribute or variable name contains illegal characters' /
if (message .ne. 'Unknown Error')
+ call errorc('nf_strerror on bad error status returned: ',
+ message)
if (err .ne. 0)
+ call errore('nf_open: ', err)
if (err .ne. NF_EPERM)
+ call error('nf_redef of read-only file should fail')
if (ncid2 .eq. ncid)
+ call error(
+ 'netCDF IDs for first and second nf_open calls should differ')
static schar schars[] = {
SCHAR_MIN, SCHAR_MIN +1,
#if SHORT_MAX <= X_SHORT_MAX
SHORT_MIN, SHORT_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MIN - 1, SCHAR_MIN, SCHAR_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MAX - 1, SCHAR_MAX, SCHAR_MAX + 1,
#else
X_SHORT_MIN, X_SHORT_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MIN - 1, SCHAR_MIN, SCHAR_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MAX - 1, SCHAR_MAX, SCHAR_MAX + 1,
#if INT_MAX <= X_INT_MAX
INT_MIN, INT_MIN +1,
# if SHORT_MAX < X_INT_MAX
static schar schars[] = {
SCHAR_MIN, SCHAR_MIN +1,
#if SHORT_MAX <= X_SHORT_MAX
SHORT_MIN, SHORT_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MIN - 1, SCHAR_MIN, SCHAR_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MAX - 1, SCHAR_MAX, SCHAR_MAX + 1,
#else
X_SHORT_MIN, X_SHORT_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MIN - 1, SCHAR_MIN, SCHAR_MIN + 1,
# if SCHAR_MAX < X_SHORT_MAX
SCHAR_MAX - 1, SCHAR_MAX, SCHAR_MAX + 1,
#if INT_MAX <= X_INT_MAX
INT_MIN, INT_MIN +1,
# if SHORT_MAX < X_INT_MAX
error("nc_open: %s", nc_strerror(err));
for (i = 0; i < NVARS; i++) {
canConvert = (var_type[i] == NC_CHAR) == (NCT_TEXT == NCT_TEXT);
for (j = 0; j < var_rank[i]; j++)
error("bad var id: status = %d", err);
for (j = 0; j < var_rank[i]; j++) {
}
for (j = 0; j < var_nels[i]; j++) {
} else {
nok++;
error("nc_open: %s", nc_strerror(err));
for (i = 0; i < NVARS; i++) {
canConvert = (var_type[i] == NC_CHAR) == (NCT_SCHAR == NCT_TEXT);
for (j = 0; j < var_rank[i]; j++)
error("bad var id: status = %d", err);
for (j = 0; j < var_rank[i]; j++) {
<font face="Arial">call
h5pset_szip_f(crp_list,H5_SZIP_AK13_OM_F+H5_SZIP_NN_OM_F,8,hdf5err)</font><br>
<br>
+ if [ -n "$HDF5" ] ; then<br>
+ echo "Will use HDF5 in dir: $HDF5"<br>
+ else<br>
+ echo "Will configure for use without HDF5"<br>
+ fi<br>
+ <br>
+ if [ -n "$PHDF5" ] ; then<br>
+ echo "Will use PHDF5 in dir: $PHDF5"<br>
$sw_netcdf_path = "" ;<br>
+ $sw_hdf5_path ="";<br>
+ $sw_phdf5_path=""; <br>
$sw_os = "ARCH" ; # ARCH will match any<br>
p = p1 - 1
ranks(p1) = mod( p , ntasks_x ) * ntasks_y + p / ntasks_x
CASE ( DATA_ORDER_ZXY )
idim = ed2-sd2+1
jdim = ed3-sd3+1
kdim = ed1-sd1+1
c_kds = sd1 ; c_kde = ed1
CASE ( DATA_ORDER_ZYX )
idim = ed3-sd3+1
jdim = ed2-sd2+1
kdim = ed1-sd1+1
c_kds = sd1 ; c_kde = ed1
CASE ( DATA_ORDER_XYZ )
idim = ed1-sd1+1
jdim = ed2-sd2+1
INTEGER, PARAMETER :: int_num_handles = 99
INTEGER, PARAMETER :: max_servers = int_num_handles+1 ! why +1?
! Server polling is disabled, so cycle through servers:
prev_server_for_handle = mod ( prev_server_for_handle + 1 , nio_groups )
server_for_handle( dhandle ) = prev_server_for_handle+1
ENDIF
! and designating the groups of i/o tasks
DO i = ncompute_tasks+1, ntasks, nio
DO j = i, i+nio-1
icolor(j) = ii
ENDDO
ii = ii+1
do j=1,n_groups
i=ncompute_tasks+j*nio-1
if(me+1==i) then
in_avail=.true. ! I/O server masters are in mpi_comm_avail
CASE ( DATA_ORDER_XYZ )
grid%sm31 = grid%sm31 + dx
grid%em31 = grid%em31 + dx
grid%sm32 = grid%sm32 + dy
grid%em32 = grid%em32 + dy
grid%sp31 = grid%sp31 + dx
grid%ep31 = grid%ep31 + dx
grid%sp32 = grid%sp32 + dy
grid%ep32 = grid%ep32 + dy
grid%sd31 = grid%sd31 + dx
CASE ( DATA_ORDER_YXZ )
grid%sm31 = grid%sm31 + dy
grid%em31 = grid%em31 + dy
grid%sm32 = grid%sm32 + dx
grid%em32 = grid%em32 + dx
! hdrbuf(5) = LEN(TRIM(DateStr))
! hdrbuf(6:5+n1) = DateStr ! n1 = LEN(TRIM(DateStr)) + 1
! hdrbuf(6+n1) = LEN(TRIM(VarName))
! hdrbuf(7+n1:6+n1+n2) = VarName ! n2 = LEN(TRIM(VarName)) + 1
! hdrbuf(7+n1+n2) = FieldType
! hdrbuf(8+n1+n2) = LEN(TRIM(MemoryOrder))
! hdrbuf(9+n1+n2:8+n1+n2+n3) = MemoryOrder ! n3 = LEN(TRIM(MemoryOrder)) + 1
! hdrbuf(9+n1+n2+n3) = LEN(TRIM(Stagger))
! hdrbuf(9+n1+n2+n3:8+n1+n2+n3+n4) = Stagger ! n4 = LEN(TRIM(Stagger)) + 1
! hdrbuf(9+n1+n2+n3+n4) = LEN(TRIM(DimNames(1)))
! hdrbuf(9+n1+n2+n3+n4:8+n1+n2+n3+n4+n5) = DimNames(1) ! n5 = LEN(TRIM(DimNames(1))) + 1
! hdrbuf(9+n1+n2+n3+n4+n5) = LEN(TRIM(DimNames(2)))
! hdrbuf(9+n1+n2+n3+n4+n5:8+n1+n2+n3+n4+n5+n6) = DimNames(2) ! n6 = LEN(TRIM(DimNames(2))) + 1
! hdrbuf(9+n1+n2+n3+n4+n5+n6) = LEN(TRIM(DimNames(3)))
! hdrbuf(9+n1+n2+n3+n4+n5+n6:8+n1+n2+n3+n4+n5+n6+n7) = DimNames(3) ! n7 = LEN(TRIM(DimNames(3))) + 1
IF ( ( config_flags%polar ) .AND. ( flag_mf_xy .EQ. 1 ) ) THEN
DO j=max(jds+1,jts),min(jde-1,jte)
grid%ght_gc(i,k,j) = grid%ght_gc(i,k-1,j) - &
R_d / g * 0.5 * ( grid%t_gc(i,k ,j) * ( 1 + 0.608 * grid%qv_gc(i,k ,j) ) + &
grid%t_gc(i,k-1,j) * ( 1 + 0.608 * grid%qv_gc(i,k-1,j) ) ) * &
LOG ( grid%p_gc(i,k,j) / grid%p_gc(i,k-1,j) )
IF ( grid%p_gc(i_valid,num_metgrid_levels/2,j_valid) .LT. grid%p_gc(i_valid,num_metgrid_levels/2+1,j_valid) ) THEN
IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE
ptemp = ((grid%p_gc(i,k,j) - grid%pd_gc(i,k,j)) + (grid%p_gc(i,k+1,j) - grid%pd_gc(i,k+1,j)))/2
IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE
grid%landusef(i,grid%iswater,j) = grid%landusef(i,grid%iswater,j) + &
print *,'Noah -> Noah: bad soil moisture at i,j = ',i,j,grid%smois(i,:,j)
iicount = iicount + 1
grid%smois(i,:,j) = MAX ( grid%smois(i,:,j) , 0.005 )
! grid%smois(i,:,j) = MAX ( grid%smois(i,:,j) + lqmi(grid%isltyp(i,j)) , 0.005 )
DO i=its,itf
muu(i,j) = 0.5*(mu(i,j)+mu(i-1,j)+mub(i,j)+mub(i-1,j))
DO j=jts,jtf
DO i=its+1,itf
muu(i,j) = 0.5*(mu(i,j)+mu(i-1,j)+mub(i,j)+mub(i-1,j))
ENDDO
DO j=jts,jtf
! muu(i,j) = mu(i,j) +mub(i,j)
! fix for periodic b.c., 13 march 2004, wcs
muu(i,j) = 0.5*(mu(i,j)+mu(im,j)+mub(i,j)+mub(im,j))
DO i=its,itf-1
muu(i,j) = 0.5*(mu(i,j)+mu(i-1,j)+mub(i,j)+mub(i-1,j))
DO j=jts,jtf
! muu(i,j) = mu(i-1,j) +mub(i-1,j)
! fix for periodic b.c., 13 march 2004, wcs
DO i=i_start, i_end
muts(i,j)=mub(i,j)+mu_2(i,j)
! u coupled with my
u_2(i,k,j) = (msfuy(i,j)*u_2(i,k,j) + u_save(i,k,j)*muu(i,j))/muus(i,j)
! w coupled with my
w_2(i,k,j) = (msfty(i,j)*w_2(i,k,j) + w_save(i,k,j)*mut(i,j))/muts(i,j)
ph_2(i,k,j) = ph_2(i,k,j) + ph_save(i,k,j)
ww(i,k,j) = ww(i,k,j) + ww1(i,k,j)
ENDDO
DO i = i_start, i_end
t_2(i,k,j) = (t_2(i,k,j) + t_save(i,k,j)*mut(i,j))/muts(i,j)
DO i = i_start, i_end
t_2(i,k,j) = (t_2(i,k,j) + t_save(i,k,j)*mut(i,j))/muts(i,j)
t_2(i,k,j) = (t_2(i,k,j) - dts*number_of_small_timesteps*mut(i,j)*h_diabatic(i,k,j) &
+ t_save(i,k,j)*mut(i,j))/muts(i,j)