Changeset 5101 for LMDZ6/branches/Amaury_dev/tools/make_sso
- Timestamp:
- Jul 23, 2024, 8:22:55 AM (17 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/tools/make_sso/make_sso_SpherePack.f90
r5100 r5101 136 136 msg='Missing or wrong "-m" option ; can be "noro", "spec" or a mask file' 137 137 CALL err(nf90_open(fmsk,nf90_nowrite,fID)/=nf90_noerr,msg) 138 CALL nc( NF90_INQ_VARID(fID,"MaskOcean",vID),"MaskOcean") !--- MASK ID139 CALL nc( NF90_INQUIRE_VARIABLE(fID,vID,dimids=dIDs)) !--- DIMS IDS138 CALL nc(nf90_inq_varid(fID,"MaskOcean",vID),"MaskOcean") !--- MASK ID 139 CALL nc(nf90_inquire_variable(fID,vID,dimids=dIDs)) !--- DIMS IDS 140 140 CALL nc(nf90_inquire_dimension(fID,dIDs(1),len=nlon_ou),'x')!--- NB LONG 141 141 CALL nc(nf90_inquire_dimension(fID,dIDs(2),len=nlat_ou),'y')!--- NB LAT … … 151 151 152 152 !--- GET HEIGHT VARIABLE AND ITS DIMENSIONS 153 CALL nc( NF90_INQ_VARID(fID,vnam,vID))154 CALL nc( NF90_INQUIRE_VARIABLE(fID,vID,dimids=dIDs))153 CALL nc(nf90_inq_varid(fID,vnam,vID)) 154 CALL nc(nf90_inquire_variable (fID,vID,dimids=dIDs)) 155 155 CALL nc(nf90_inquire_dimension(fID,dIDs(1),len=nlon_in,name=lonn)) 156 156 CALL nc(nf90_inquire_dimension(fID,dIDs(2),len=nlat_in,name=latn)) … … 159 159 nlon_in=nlon_in-1 !--- DISTINCT POINTS NUMBER 160 160 ALLOCATE(lon_in(nlon_in+1),lat_in(nlat_in),h(nlon_in+1,nlat_in)) 161 CALL nc( NF90_INQ_VARID(fID,lonn,loID) ,lonn)162 CALL nc( NF90_INQ_VARID(fID,latn,laID) ,latn)161 CALL nc(nf90_inq_varid(fID,lonn,loID) ,lonn) 162 CALL nc(nf90_inq_varid(fID,latn,laID) ,latn) 163 163 CALL nc(nf90_get_var (fID,loID,lon_in) ,lonn) 164 164 CALL nc(nf90_get_var (fID,laID,lat_in) ,latn) … … 369 369 msg='Missing or wrong "-m" option ; can be "noro", "spec" or a mask file' 370 370 CALL err(nf90_open(fmsk,nf90_nowrite,fID)/=nf90_noerr,msg) 371 CALL nc( NF90_INQ_VARID(fID,"MaskOcean",vID),"MaskOcean") !--- MASK ID371 CALL nc(nf90_inq_varid(fID,"MaskOcean",vID),"MaskOcean") !--- MASK ID 372 372 CALL nc(nf90_get_var(fID,vID,msko(:,:))) !--- MASK 373 373 CALL nc(nf90_close(fID)) … … 460 460 461 461 CALL nc(nf90_def_dim(fID,'x',nlon_ou,xID)) 462 CALL nc(nf90_def_var(fID,'x', NF90_REAL,xID,loID) ,'x')462 CALL nc(nf90_def_var(fID,'x',nf90_real,xID,loID) ,'x') 463 463 CALL nc(nf90_put_att(fID,loID,'long_name','Longitude'),'x') 464 464 CALL nc(nf90_put_att(fID,loID,'units','degrees_east') ,'x') 465 465 466 466 CALL nc(nf90_def_dim(fID,'y',nlat_ou,yID)) 467 CALL nc(nf90_def_var(fID,'y', NF90_REAL,yID,laID) ,'y')467 CALL nc(nf90_def_var(fID,'y',nf90_real,yID,laID) ,'y') 468 468 CALL nc(nf90_put_att(fID,laID,'long_name','Latitude') ,'y') 469 469 CALL nc(nf90_put_att(fID,laID,'units','degrees_north'),'y') 470 470 471 CALL nc(nf90_def_var(fID,'mask', NF90_REAL,[xID,yID],mskID),'mask')471 CALL nc(nf90_def_var(fID,'mask',nf90_real,[xID,yID],mskID),'mask') 472 472 IF(fmsk=='noro') & 473 CALL nc(nf90_def_var(fID,'Zphi', NF90_REAL,[xID,yID],phiID),'Zphi')474 CALL nc(nf90_def_var(fID,'Zmea', NF90_REAL,[xID,yID],meaID),'Zmea')475 CALL nc(nf90_def_var(fID,'mu' , NF90_REAL,[xID,yID], muID),'mu' )476 CALL nc(nf90_def_var(fID,'Zsig', NF90_REAL,[xID,yID],sigID),'Zsig')477 CALL nc(nf90_def_var(fID,'Zgam', NF90_REAL,[xID,yID],gamID),'Zgam')478 CALL nc(nf90_def_var(fID,'Zthe', NF90_REAL,[xID,yID],theID),'Zthe')479 CALL nc(nf90_def_var(fID,'Zpic', NF90_REAL,[xID,yID],picID),'Zpic')480 CALL nc(nf90_def_var(fID,'Zval', NF90_REAL,[xID,yID],valID),'Zval')473 CALL nc(nf90_def_var(fID,'Zphi',nf90_real,[xID,yID],phiID),'Zphi') 474 CALL nc(nf90_def_var(fID,'Zmea',nf90_real,[xID,yID],meaID),'Zmea') 475 CALL nc(nf90_def_var(fID,'mu' ,nf90_real,[xID,yID], muID),'mu' ) 476 CALL nc(nf90_def_var(fID,'Zsig',nf90_real,[xID,yID],sigID),'Zsig') 477 CALL nc(nf90_def_var(fID,'Zgam',nf90_real,[xID,yID],gamID),'Zgam') 478 CALL nc(nf90_def_var(fID,'Zthe',nf90_real,[xID,yID],theID),'Zthe') 479 CALL nc(nf90_def_var(fID,'Zpic',nf90_real,[xID,yID],picID),'Zpic') 480 CALL nc(nf90_def_var(fID,'Zval',nf90_real,[xID,yID],valID),'Zval') 481 481 482 482 CALL nc(nf90_put_att(fID,mskID,'long_name','Fractional land mask' ),'mask') … … 834 834 IF(fnam/='') msg=TRIM(msg)//' for file "'//TRIM(fnam)//'"' 835 835 IF(PRESENT(var)) msg=TRIM(msg)//'" and variable "'//TRIM(var)//'"' 836 WRITE(*,*)TRIM(msg)//': '// NF90_STRERROR(ncres); STOP836 WRITE(*,*)TRIM(msg)//': '//nf90_strerror(ncres); STOP 837 837 END IF 838 838
Note: See TracChangeset
for help on using the changeset viewer.
