adding variables in netcdf4 fortran 77 -

-

i'm trying save outputs of model (written in fortran 77) in netcdf4 format. i'm using gfortran. since i'm fresher in fortran, wrote simple code in create dummy variables , add test netcdf file emulating dimensions,variables , attributes of model. works perfectly.

i proceeded next step knowing code works. well, thing same code not work , cannot find out why. model code old , "organized" in 1 single .for file. interestingly enough, main code written subroutine. basic structure of program is:

  implicit real*8(a-h,o-z)    character dummy*80   ...   call main(delt,npl)         end        subroutine main(delt, npl)   *variable declarations*   ...   end

my approach create netcdf file outside "main" subroutine. works , netcdf file created. inside "main" define dimensions , variable (since they're based on data generated inside subroutine). code works charm when add dimensions:

  setting dimensions   nc_status = nf_def_dim(ncid,'parcel_id', ntot, prcl_dimid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)   nc_status = nf_def_dim(ncid, 'time', icomp, time_dimid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)   end define mode.   nc_status = nf_enddef(ncid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)   close file   nc_status = nf_close(ncid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)

the nc file generated has right dimensions right size. problem comes when try add variables:

  dimensions   nc_status = nf_def_dim(ncid,'parcel_id', ntot, prcl_dimid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)   nc_status = nf_def_dim(ncid, 'time', icomp, time_dimid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)   dimids(1) = prcl_dimid   dimids(2) = time_dimid    variables   **nc_status = nf_def_var(ncid, 'latitude', nf_float, 2, dimids,   +     lat_varid)**   if (nc_status .ne. nf_noerr) call handle_err(nc_status)    end define mode.   nc_status = nf_enddef(ncid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)    close file   nc_status = nf_close(ncid)   if (nc_status .ne. nf_noerr) call handle_err(nc_status)

i following error bold-font part of code.

error: netcdf: not valid data type or _fillvalue type mismatch

first, code (exactly that) works when create simple fortran file dummy values. second, no fillvalue defined me. , third, nf_float valid data type according unidata website.

i've tried create whole netcdf4 inside "main" subroutine same error. i've used imagination trying figure out happening and, apparently, no 1 else had same problem (or @ least did not mention online...). i'm beginner in fortran , possible.

all best , sorry if stupid question you. hasnt been stupid me @ all!


Comments

Post a Comment

Popular posts from this blog

java - Oracle EBS .ClassNotFoundException: oracle.apps.fnd.formsClient.FormsLauncher.class ERROR -

-
i need in configuring java on machine because, when try loading forms, loader goes in circles , gives error: java.lang.classnotfoundexception: oracle.apps.fnd.formsclient.formslauncher.class @ sun.plugin2.applet.applet2classloader.findclass(unknown source) @ sun.plugin2.applet.plugin2classloader.loadclass0(unknown source) @ sun.plugin2.applet.plugin2classloader.loadclass(unknown source) @ sun.plugin2.applet.plugin2classloader.loadclass(unknown source) @ java.lang.classloader.loadclass(unknown source) @ sun.plugin2.applet.plugin2classloader.loadcode(unknown source) @ sun.plugin2.applet.plugin2manager.createapplet(unknown source) @ sun.plugin2.applet.plugin2manager$appletexecutionrunnable.run (unknown source) @ java.lang.thread.run(unknown source) i have tried re-installing firefox, older versions of java runtime environments, disabling proxy. what missing? application working on other machines. i spent around 100 hours finding s...
Read more

c# - how to use buttonedit in devexpress gridcontrol -

-
Image
i trying create gridcontrol buttonedit in 1 of columns. when user clicks button edit, popups form select product. when selection done on popup, returns datarow selection main grid below. but when column loses focus, value wrote column disappear. here code creates gridcontrol's data , buttonedit's click event. private void frmsiparisnew_load(object sender, eventargs e) { dt = new datatable(); dt.columns.add("malzeme_kodu",typeof(string)); dt.columns.add("malzeme_aciklama", typeof(string)); dt.columns.add("adet", typeof(decimal)); dt.columns.add("birim", typeof(string)); dt.columns.add("fiyat", typeof(decimal)); dt.columns.add("kur", typeof(string)); dt.columns.add("tutar", typeof(decimal)); datarow dr = dt.newrow(); dt.rows.add(dr); gc.datasource = dt; } private void repositoryitembuttoneditmal...
Read more

How do you convert a timestamp into a datetime in python with the correct timezone? -

-
you naively expect work: import pytz datetime import datetime def tz_datetime_from_timestamp(timestamp): """convert timestamp datetime """ tz = pytz.timezone('australia/perth') server_time = datetime.utcnow() delta = tz.utcoffset(server_time) value = datetime.utcfromtimestamp(timestamp).replace(tzinfo=tz) return value + delta print tz_datetime_from_timestamp(1416387060) and convert timestamp 1416387060 wed nov 19 16:51:00 2014 gmt+8:00. ...but not. prints: 2014-11-19 16:51:00+07:43 the timezone of australia/perth not gmt+7:43. it gmt+8:00. unequivocally stated on australian government website http://www.australia.gov.au/about-australia/our-country/time : awst equal coordinated universal time plus 8 hours (utc +8). so, heck pytz pulling 7.43 from? well, turns out pytz pulls data http://www.iana.org/time-zones time zone database, , database states: # western australia # # rule name type in on @...
Read more