Changeset 5101 for LMDZ6/branches/Amaury_dev/tools
- Timestamp:
- Jul 23, 2024, 8:22:55 AM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/tools
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/tools/fcm/lib/Ecmwf/Fortran90_stuff.pm
r1578 r5101 601 601 next if(/^ *$/ | /^ *!/); 602 602 # The following two substitutions should be restored at end of processing 603 s/(\'[^!]*)!+(.*\')/$1\ £$2/; # Protect against mischief604 s/(["][^!]*)!+(.*["])/$1\ £$2/; # Protect against mischief603 s/(\'[^!]*)!+(.*\')/$1\�$2/; # Protect against mischief 604 s/(["][^!]*)!+(.*["])/$1\�$2/; # Protect against mischief 605 605 $im=$im+/JPIM\b/i unless($im); 606 606 $rb=$rb+/JPRB\b/i unless($rb); … … 629 629 next unless (/^\#/); 630 630 if(@pars) { 631 s/^#include +"tsmbkind.h"/USE PARKIND1 ,ONLY 631 s/^#include +"tsmbkind.h"/USE PARKIND1 ,ONLY: $string/ ; 632 632 } 633 633 else { … … 635 635 } 636 636 # if($rh) { 637 s/^#include +"hugekind.h"/USE PARKIND2 ,ONLY 637 s/^#include +"hugekind.h"/USE PARKIND2 ,ONLY: JPRH/ ; 638 638 # } 639 639 # else { … … 1418 1418 if($hook_status == 0) { # $hook_status == 0 means we have not done anything yet 1419 1419 if($cont eq 'USE') { # Add USE YOMHOOK as second use statement 1420 $href->{post_insert}="USE YOMHOOK ,ONLY 1420 $href->{post_insert}="USE YOMHOOK ,ONLY: LHOOK, DR_HOOK\n"; 1421 1421 $hook_status=1; 1422 1422 } 1423 1423 elsif($cont eq 'IMPLICIT NONE') { # No previous USE, add USE YOMHOOK before IMPLICIT NONE 1424 $href->{pre_insert} ="USE PARKIND1 ,ONLY 1425 "USE YOMHOOK ,ONLY 1424 $href->{pre_insert} ="USE PARKIND1 ,ONLY: JPRB\n". 1425 "USE YOMHOOK ,ONLY: LHOOK, DR_HOOK\n"; 1426 1426 $hook_status=1; 1427 1427 } … … 2471 2471 s/($name)%$name/$1/g; 2472 2472 $_=uc($_); 2473 # print "PROPAG $called $_ ££@inpars \n";2473 # print "PROPAG $called $_ �� @inpars \n"; 2474 2474 my @call_args=split ',' , $_; 2475 2475 my $i=0; -
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.