// test for julian()
month     = 1       ; // <integer> month
day       = 9       ; // <integer> day
year      = 1974    ; // <integer> year
hour      = 9       ; // <integer> hour
minute    = 0       ; // <integer> minute
second    = 0       ; // <integer> second
ier       = 0       ; // <integer> error flag (0 = OK)
jdate     = 0.0     ; // <real> julian date             

[ier,jdate] = julian(day,month,year,hour,minute,second);

// output results
ier
jdate

// test for call_mcd()
// memory allocation for output 
//
ier = 0;
i=1:5  ; instvar(i)=0.0;
i=1:5  ; meanvar(i)=0.0;
i=1:100; extvar(i) =0.0;
seedout = 0.0;
//
// definition of inputs
//
i=1:4; keys(i)=0;
keys(1) = 3; // zkey
keys(2) = 1; // hireskey
keys(3) = 1; // datakey
keys(4) = 1; // perturkey
scena   = 1; // scenario
i=1:100; extvarkeys(i)=0; // array output type 
extvarkeys(10)=0;
i=1:5; point(i)=0.0;
point(1) = 100.0; // xz
point(2) = 6.0;   // xlon
point(3) = 5.0;   // xlat
point(4) = 120.0; // xdate
point(5) = 2.0;   // localtime
dset ='MCD_DATA/'; // <char*50> data set
seedin   = 7.0;
gwlength = 0.0;
[ier,instvar,meanvar,extvar,seedout]=call_mcd(keys,scena,extvarkeys,point,dset,seedin,gwlength);

// output results
ier
instvar
meanvar
extvar
seedout


