text:
HANSARD REVISE * NUMERO 35
Le vendredi 21 novembre 1997
LES COMITES DE LA CHAMBRE
M. Dennis J. Mills
Mme Jocelyne Girard-Bujold
Mme Caroline St-Hilaire
M. Guy Saint-Julien
LE WYE MARSH WILDLIFE CENTRE
LA COMMISSION ROYALES SUR LES PEUPLES AUTOCHTONES
Mme Nancy Karetak-Lindell
Mme Madeleine Dalphond-Guiral
L'INSTITUT DES TECHNOLOGIES DE FABRICATION INTEGREE
Mme Rose-Marie Ur
L'hon. Pierre S. Pettigrew
L'hon. Pierre S. Pettigrew
Mme Jocelyne Girard-Bujold
Mme Jocelyne Girard-Bujold
L'hon. Pierre S. Pettigrew
L'hon. Pierre S. Pettigrew
L'EXPORTATION DE MARCHANDISES MILITAIRES
REPONSE DU GOUVERNEMENT A DES PETITIONS
LES COMITES DE LA CHAMBRE
Presentation et premiere lecture
COMITES DE LA CHAMBRE
Procedure et affaires de la Chambre
Mme Rose-Marie Ur
Mme Jocelyne Girard-Bujold
M. Guy Saint-Julien
HANSARD REVISE * NUMERO 35
Le vendredi 21 novembre 1997
La seance est ouverte a 10 heures.
M. Bob Kilger (Stormont-Dundas, Lib.):
Monsieur le President, j'invoque le Reglement.
Le president suppleant (M. McClelland):
La Chambre a entendu la requete.
Y a-t-il consentement unanime?
LES COMITES DE LA CHAMBRE
M. Dennis J. Mills (Broadview-Greenwood, Lib.):
M. Julian Reed (secretaire parlementaire du ministre du Commerce international, Lib.):
Ces entreprises et les fabricants canadiens de materiel de telecommunications sont les meilleurs au monde.
Apres deux ans de negociations, un accord a ete conclu le 15 fevrier 1997.
Mme Elinor Caplan (Thornhill, Lib.):
Monsieur le President, je suis heureuse de prendre la parole sur ce projet de loi.
En raison de ces objectifs valables, le projet de loi merite d'etre appuye.
Plusieurs mesures ont deja ete prises.
Il prevoit la mise en oeuvre des dispositions de ces accords.
Recemment, j'ai rencontre l'agent de developpement economique de la ville de Markham.
On compte environ 700 entreprises dans ce secteur.
On y a cree un centre ou les vieux ordinateurs sont remis en etat.
Le Canada est un pionnier.
Nous sommes les leaders dans ce domaine.
La facon d'y parvenir est de laisser jouer les forces du marche.
M. Dick Harris (Prince George-Bulkley Valley, Ref.):
Il ameliorera peut-etre le systeme globalement.
La deputee pourrait-elle commenter ces observations s'il-vous-plait?
M. Peter Stoffer (Sackville-Eastern Shore, NPD):
Monsieur le President, j'ai ecoute les commentaires de la deputee liberale.
Monsieur le President, je suis heureuse de repondre au depute.
Je suis persuadee que, en fin de compte, le Canada et les Canadiens seront beneficiaires.
Nous profiterons de cela de deux facons.
Cela signifie plus d'emplois et de meilleurs emplois pour tout le monde.
Je pense que ce sera le resultat de la concurrence.
Mme Helene Alarie (Louis-Hebert, BQ):
Les telecommunications sont en constante evolution au Canada.
Les reseaux de communications classiques cedent leur place a de nouvelles technologies.
En 1993, le trafic international des communications a atteint 47,7 milliards de minutes.
On remarque que ce projet de loi poursuit donc deux objectifs.
Quant a l'article 46.6,
Que signifie le terme abordabilite du service de base?
Cela devient donc un controle sur les individus et la societe.
M. Dick Harris (Prince George-Bulkley Valley, Ref.):
Notre parti est partisan de la libre entreprise.
Mon autre reserve concerne le financement de l'operation.
Le gouvernement y voit la une occasion d'elargir son cercle d'amis.
Le president suppleant (M. McClelland):
Mme Jocelyne Girard-Bujold (Jonquiere, BQ):
Toutes les questions qu'elle a posees concernant ce projet de loi sont courantes.
A ce niveau, nous sommes donc bien places.
M. Rob Anders (Calgary-Ouest, Ref.):
J'imagine tout de suite un CRTC encore plus imposant.
Monsieur le President, je suis heureux de pouvoir intervenir.
Le president suppleant (M. McClelland):
La Chambre est-elle prete a se prononcer?
Le president suppleant (M. McClelland):
Plait-il a la Chambre d'adopter la motion?
Le president suppleant (M. McClelland):
Que tous ceux qui sont en faveur de la motion veuillent bien dire oui.
Le president suppleant (M. McClelland):
Que tous ceux qui sont contre veuillent bien dire non.
Le president suppleant (M. McClelland):
A mon avis, les oui l'emportent.
M. Bob Kilger (Stormont-Dundas, Lib.):
Monsieur le President, j'invoque le Reglement.
Le president suppleant (M. McClelland):
La Chambre a entendu la proposition du whip en chef du gouvernement.
Y a-t-il unanimite?
Monsieur le President, j'invoque le Reglement.
Le president suppleant (M. McClelland):
La Chambre a entendu la motion.
Y a-t-il unanimite?
Le president suppleant (M. McClelland):
Il n'y a pas unanimite.
M. Peter Stoffer (Sackville-Eastern Shore, NPD):
Monsieur le President, j'invoque le Reglement.
Je donne mon consentement.
Le president suppleant (M. McClelland):
M. Bob Kilger (Stormont-Dundas, Lib.):
Le president suppleant (M. McClelland):
La Chambre a entendu le whip en chef du gouvernement.
Y a-t-il unanimite?
Le president suppleant (M. McClelland):
Le whip en chef du gouvernement invoque le Reglement.
M. Bob Kilger (Stormont-Dundas, Lib.):
Le president suppleant (M. McClelland):
Y a-t-il consentement unanime?
L'hon. Andy Mitchell (secretaire d'Etat (Parcs), Lib.):
Monsieur le President, je voudrais proposer un amendement de forme.
Au fond, nous voulons garantir que les deux versions, anglaise et francaise, correspondent.
flore, ainsi que les ressources naturelles exceptionnel-.
M. Paul Crete (Kamouraska-Riviere-du-Loup-Temiscouata-Les Basques, BQ):
Pour la premiere fois, le Quebec garde la propriete des fonds marins.
Je pense que c'est une belle etape et un bon resultat.
Dans ce sens-la, c'est tres souhaitable.
Le president suppleant (M. McClelland):
M. Dick Harris (Prince George-Bulkley Valley, Ref.):
Cette recompense est accordee a ceux qui ont endure de grande souffrances avec heroisme.
Tu as vraiment vecu en accord avec la promesse et la loi scoutes.
M. Eugene Bellemare (Carleton-Gloucester, Lib.):
J'ai egalement eu le plaisir de les accueillir a la Chambre des communes.
Votre compassion et votre devouement sont d'un grand secours dans des moments difficiles.
Mme Caroline St-Hilaire (Longueuil, BQ):
L'avenir appartient aux audacieuses et aux audacieux.
Les gens de la vallee de la Matapedia font le pari de l'audace.
Mme Elinor Caplan (Thornhill, Lib.):
Je prends aujourd'hui la parole a la Chambre pour signaler une reussite.
Il s'agit d'une reduction de 6 p. 100.
Le depute d'Abitibi.
M. Guy Saint-Julien (Abitibi, Lib.):
Pour reservation, telephonez au 996-6740.
Venez voir la nouvelle reforme du bingo.
M. Jay Hill (Prince George-Peace River, Ref.):
Monsieur le President, le gouvernement se defait des petits aeroports regionaux dans tout le Canada.
Cette petite localite au grand coeur ne compte que 4 500 habitants.
LE WYE MARSH WILDLIFE CENTRE
M. Paul DeVillers (Simcoe-Nord, Lib.):
Mme Marlene Catterall (Ottawa-Ouest-Nepean, Lib.):
M. Grant McNally (Dewdney-Alouette, Ref.):
Il offrira des debouches de plusieurs millions de dollars aux entrepreneurs canadiens et chinois.
M. Robert Bertrand (Pontiac-Gatineau-Labelle, Lib.):
Le nombre de visiteurs, cette annee, a augmente de 238 p. 100.
M. Gordon Earle (Halifax-Ouest, NPD):
Des excuses s'imposent.
De jeunes autochtones ont ete arraches a leur famille et a leur collectivite.
Leur langue, leur culture, leurs coutumes et leurs valeurs ont ete supprimees.
Des excuses de la part du gouvernement constitueraient un petit remede qui favoriserait la guerison.
Moins que cela serait deraisonnable, impardonnable et inacceptable.
LA COMMISSION ROYALES SUR LES PEUPLES AUTOCHTONES
Mme Nancy Karetak-Lindell (Nunavut, Lib.):
Le depute s'exprime en inuktitut: ]
Mme Madeleine Dalphond-Guiral (Laval-Centre, BQ):
Monsieur le President, ou sont passes nos epouvantails?
Mais cet automne, nos epouvantails se font etrangement discrets.
Il est vrai que les bonnes nouvelles abondent.
Il y a de quoi refroidir d'ardeur de nos epouvantails economiques.
M. Peter MacKay (Pictou-Antigonish-Guysborough, PC):
Monsieur le President, c'est la semaine de la justice reparatrice.
La veritable justice mise sur les forces au lieu de perpetuer l'echec.
L'INSTITUT DES TECHNOLOGIES DE FABRICATION INTEGREE
Mme Sue Barnes (London-Ouest, Lib.):
Mme Rose-Marie Ur (Lambton-Kent-Middlesex, Lib.):
Ils peuvent etre fiers de leur travail pour les Nations Unies.
M. Ken Epp (Elk Island, Ref.):
Les Canadiens veulent envoyer leurs colis.
L'hon. Lawrence MacAulay (ministre du Travail, Lib.):
M. Ken Epp (Elk Island, Ref.):
L'hon. Lawrence MacAulay (ministre du Travail, Lib.):
Laissons fonctionner le systeme.
M. Ken Epp (Elk Island, Ref.):
Monsieur le President, les organismes de charite sont le plus durement frappes.
Prenez par exemple la Old Brewery Mission, a Montreal, un refuge pour sans-abri.
L'hon. Lawrence MacAulay (ministre du Travail, Lib.):
Les reformistes devraient avoir honte.
M. Jim Hart (Okanagan-Coquihalla, Ref.):
Monsieur le President, les petites et moyennes entreprises sont le moteur de notre economie.
Pendant des semaines, elles ont souffert de la menace d'une greve des postes.
L'hon. Lawrence MacAulay (ministre du Travail, Lib.):
Monsieur le President, le regime de negociation collective donne aux travailleurs le droit de greve.
M. Jim Hart (Okanagan-Coquihalla, Ref.):
Quatre greves en dix ans.
Sept mois de negociations qui se sont soldees par un echec.
L'hon. Lawrence MacAulay (ministre du Travail, Lib.):
Depuis les sept derniers mois, le nombre de jours perdus est inferieur a la moyenne.
Mme Suzanne Tremblay (Rimouski-Mitis, BQ):
Monsieur le President, ma question s'adresse au premier ministre suppleant.
L'affaire d'Option Canada sent de plus en plus mauvais.
L'hon. Andrew Mitchell (secretaire d'Etat (Parcs), Lib.):
Mme Suzanne Tremblay (Rimouski-Mitis, BQ):
L'hon. Andrew Mitchell (secretaire d'Etat (Parcs), Lib.):
Ce parti n'a jamais pu se resigner a accepter ce verdict.
Les bloquistes poursuivent la lutte d'il y a deux ans.
M. Bernard Bigras (Rosemont, BQ):
Monsieur le President, ma question s'adresse au vice-premier ministre.
Il demande au gouvernement federal de ne pas se defiler de ses engagements anterieurs.
L'hon. Allan Rock (ministre de la Sante, Lib.):
M. Bernard Bigras (Rosemont, BQ):
L'hon. Allan Rock (ministre de la Sante, Lib.):
Nous avons l'intention de le faire a la conference, dans quelques semaines.
M. Bill Blaikie (Winnipeg-Transcona, NPD):
L'hon. Pierre S. Pettigrew (ministre du Developpement des ressources humaines, Lib.):
Celui-ci rencontre actuellement les autorites provinciales, les pecheurs, les associations et les representants locaux.
M. Bill Blaikie (Winnipeg-Transcona, NPD):
L'hon. Pierre S. Pettigrew (ministre du Developpement des ressources humaines, Lib.):
M. Bill Matthews (Burin-St.
M. Stan Keyes (secretaire parlementaire du ministre des Transports, Lib.):
Monsieur le President, je remercie le depute pour sa question.
M. Bill Matthews (Burin-St.
Il n'existe pas d'autres operations pour cette compagnie.
M. Stan Keyes (secretaire parlementaire du ministre des Transports, Lib.):
M. Bill Gilmour (Nanaimo-Alberni, Ref.):
Le gouvernement veut retablir les niveaux de 1990 pour l'an 2007.
La Chambre n'arrive pas a savoir quelle est la position du gouvernement.
Qui payera cette somme?
L'hon. Allan Rock (ministre de la Sante, Lib.):
M. Bill Gilmour (Nanaimo-Alberni, Ref.):
C'est le seul but des participants a cette rencontre.
Je pose la question a nouveau: qui payera le prix de cette position?
L'hon. Allan Rock (ministre de la Sante, Lib.):
Nous consultons tous les intervenants.
Les provinces participent activement au processus.
Voila quel est notre objectif.
M. Richard Marceau (Charlesbourg, BQ):
Monsieur le President, ma question s'adresse au solliciteur general.
L'hon. Andy Scott (solliciteur general du Canada, Lib.):
Nous allons tres bientot donner suite a ces recommandations.
M. Richard Marceau (Charlesbourg, BQ):
L'hon. Andy Scott (solliciteur general du Canada, Lib.):
M. Maurice Vellacott (Wanuskewin, Ref.):
L'hon. Allan Rock (ministre de la Sante, Lib.):
Monsieur le President, le depute a tort sur toute la ligne.
On a juge que nous avions raison.
Le depute a tort.
M. Maurice Vellacott (Wanuskewin, Ref.):
Monsieur le President, permettez-moi de demeurer un peu sceptique a cet egard.
Randy Conners et bien d'autres sont morts.
Ces etres chers sont morts pendant que le ministre entravait le travail du juge Krever.
L'hon. Allan Rock (ministre de la Sante, Lib.):
Nous n'avons en rien gene la commission Krever.
Mme Jocelyne Girard-Bujold (Jonquiere, BQ):
L'hon. Alfonso Gagliano (ministre des Travaux publics et des Services gouvernementaux, Lib.):
Mme Jocelyne Girard-Bujold (Jonquiere, BQ):
L'hon. Alfonso Gagliano (ministre des Travaux publics et des Services gouvernementaux, Lib.):
M. Darrel Stinson (Okanagan-Shuswap, Ref.):
Monsieur le President, le premier ministre a montre qu'il en a enfin eu assez.
L'hon. Herb Gray (vice-premier ministre, Lib.):
Nous lui savons gre de son appui.
M. Darrel Stinson (Okanagan-Shuswap, Ref.):
Le premier ministre revendique le droit de nommer des senateurs.
L'hon. Herb Gray (vice-premier ministre, Lib.):
Monsieur le President, le depute demande une modification de la Constitution.
Il s'agit la d'une proposition tres interessante.
Telle est la position du Parti reformiste sur tout.
M. Paul Crete (Kamouraska-Riviere-du-Loup-Temiscouata-Les Basques, BQ):
Monsieur le President, ma question s'adresse au ministre du Developpement des ressources humaines.
L'hon. Pierre S. Pettigrew (ministre du Developpement des ressources humaines, Lib.):
Monsieur le President, j'ai rencontre le ministre des Finances a plusieurs reprises cette semaine.
M. John McKay (Scarborough-Est, Lib.):
L'economie cree des emplois deux fois plus vite que la normale.
Et pourtant, tous les Canadiens ne profitent pas de cette relance economique.
Le taux de chomage reste coince a 9,1 p. 100.
L'hon. Paul Martin (ministre des Finances, Lib.):
Cela represente une reduction d'impots de 1,4 milliards de dollars.
Le depute d'Esquimalt-Juan de Fuca.
M. Keith Martin (Esquimalt-Juan de Fuca, Ref.):
On a relegue cette question aux oubliettes.
L'hon. David Kilgour (secretaire d'Etat (Amerique latine et Afrique), Lib.):
M. Keith Martin (Esquimalt-Juan de Fuca, Ref.):
L'hon. David Kilgour (secretaire d'Etat (Amerique latine et Afrique), Lib.):
Il peut dialoguer directement avec lui de facon ferme, mais sans chercher la confrontation.
Je pretends que cela reflete ce que nous essayons de faire.
M. Gordon Earle (Halifax-Ouest, NPD):
Ma question est claire.
M. Bernard Patry (secretaire parlementaire du ministre des Affaires indiennes et du Nord canadien, Lib.):
Monsieur le President, je remercie le depute pour sa question.
Le rapport n'est pas sur la tablette.
Au contraire, c'est un des principaux documents que nous avons.
Mme Bev Desjarlais (Churchill, NPD):
Le gouvernement liberal declare-t-il officiellement qu'il approuve cette recommandation?
M. Bernard Patry (secretaire parlementaire du ministre des Affaires indiennes et du Nord canadien, Lib.):
Nous tenons a deposer la bonne reponse.
Nous travaillons en collaboration avec les gouvernements provinciaux et les peuples autochtones.
M. Scott Brison (Kings-Hants, PC):
Monsieur le President, ma question s'adresse au secretaire parlementaire du ministre du Commerce international.
Qu'est-ce que le gouvernement entend faire pour remedier au probleme?
M. Julian Reed (secretaire parlementaire du ministre du Commerce international, Lib.):
M. Gerald Keddy (South Shore, PC):
Monsieur le President, ma question s'adresse au ministre du Commerce international.
Les produits ne sont pas toujours tous recharges.
M. Julian Reed (secretaire parlementaire du ministre du Commerce international, Lib.):
Monsieur le President, je le souhaiterais enormement.
M. Joe Jordan (Leeds-Grenville, Lib.):
Monsieur le President, c'etait hier la Journee nationale de l'enfant.
A bien y penser, cela explique peut-etre certains des comportements de la Chambre hier.
Quand allez-vous cesser de les taxer?
Le ministre de la Sante.
L'hon. Allan Rock (ministre de la Sante, Lib.):
Monsieur le President, il est certain qu'une intervention precoce est cruciale pour les enfants.
A titre de president du programme national pour les enfants...
La presidence a du mal a entendre la reponse.
J'invite le ministre de la Sante a poursuivre.
Il s'agit du Programme d'action communautaire pour les enfants.
M. Mike Scott (Skeena, Ref.):
M. Bernard Patry (secretaire parlementaire du ministre des Affaires indiennes et du Nord canadien, Lib.):
Monsieur le President, nous ne faisons rien de tel.
Mme Helene Alarie (Louis-Hebert, BQ):
Monsieur le President, ma question s'adresse au ministre des Finances.
L'hon. Paul Martin (ministre des Finances, Lib.):
M. Peter Stoffer (Sackville-Eastern Shore, NPD):
L'hon. Pierre S. Pettigrew (ministre du Developpement des ressources humaines, Lib.):
Cela ne concerne pas seulement le programme LSPA, mais se passe dans tout le Canada.
C'est comme cela que les choses fonctionnent.
M. Scott Brison (Kings-Hants, PC):
Il veut et devrait etre le verificateur de cet office.
L'hon. Paul Martin (ministre des Finances, Lib.):
Monsieur le President, ce n'est pas la position du verificateur general.
De toute maniere, la position du verificateur general sera clarifiee.
M. Ian Murray (Lanark-Carleton, Lib.):
Monsieur le President, ma question s'adresse au secretaire parlementaire du ministre de l'Industrie.
Les Canadiens s'inquietent du nombre de personnes victimes des manoeuvres frauduleuses des tele-vendeurs.
M. Walt Lastewka (secretaire parlementaire du ministre de l'Industrie, Lib.):
M. Gurmant Grewal (Surrey-Centre, Ref.):
Monsieur le President, la ministre responsable de l'ACDI dirige mal son ministere.
Le depute nous demande d'etablir un autre palier administratif.
L'EXPORTATION DE MARCHANDISES MILITAIRES
M. Daniel Turp (Beauharnois-Salaberry, BQ):
L'hon. Don Boudria (leader du gouvernement a la Chambre des communes, Lib.):
Monsieur le President, cela se rapporte au depot d'un rapport.
Mme Michelle Dockrill (Bras d'Or, NPD):
Monsieur le President, ma question s'adresse au ministre des Ressources naturelles.
Le president de Devco, M. Joe Shannon, pretend travailler tres fort pour moderniser Devco.
M. Gerry Byrne (secretaire parlementaire du ministre des Ressources naturelles, Lib.):
Je suis tres satisfait des efforts qu'ils deploient dans de tres difficiles circonstances.
Il est tres difficile en effet de soutenir la concurrence internationale dans cette industrie.
Ils font de l'excellent travail.
M. Gerry Ritz (Battlefords-Lloydminster, Ref.):
L'hon. Paul Martin (ministre des Finances, Lib.):
Je vais donc la repeter.
Cela met un terme a la periode des questions.
Je tiens a feliciter tous mes collegues a l'occasion de cet heureux anniversaire.
M. Ken Epp (Elk Island, Ref.):
Monsieur le President, j'invoque le Reglement.
Je m'interroge au sujet de votre declaration.
Je suis convaincu que le depute ne tient pas a embrouiller les choses.
Une deputee m'a donne avis de son intention de soulever la question de privilege.
Mme Bev Desjarlais (Churchill, NPD):
Les reformistes n'ont pas propose d'amendements visant a proteger ces employes.
C'est la le critere selon lequel une question de privilege doit etre evaluee.
Je considere donc l'affaire close.
M. Joe Jordan (Leeds-Grenville, Lib.):
Monsieur le President, j'invoque le Reglement.
M. Rob Anders (Calgary-Ouest, Ref.):
Comment pouvons-nous faire cela?
M. Keith Martin (Esquimalt-Juan de Fuca, Ref.):
Monsieur le President, je prends la parole au sujet du meme rappel au Reglement.
Je ne sais pas comment nous pouvons nous lever a la Chambre...
Je ne crois pas qu'il s'agisse ici d'un rappel au Reglement.
Les deputes sont libres de se lever ou de rester assis s'ils le desirent.
Aucun depute n'est tenu de se lever ou de faire quoi que ce soit.
Je crois que l'affaire est close.
REPONSE DU GOUVERNEMENT A DES PETITIONS
LES COMITES DE LA CHAMBRE
M. Maurizio Bevilacqua (Vaughan-King-Aurora, Lib.):
Le depute de Calgary-Ouest souleve la question de privilege.
M. Rob Anders (Calgary-Ouest, Ref.):
Il n'y a pas matiere a question de privilege.
Malheureusement, trop de mariages au Canada se terminent par un divorce.
COMITES DE LA CHAMBRE
M. Bob Kilger (Stormont-Dundas, Lib.) propose:
PROCEDURE ET AFFAIRES DE LA CHAMBRE
Mme Marlene Catterall (Ottawa-Ouest-Nepean, Lib.):
Monsieur le President, je propose:
M. Jack Ramsay (Crowfoot, Ref.):
Les trois premieres petitions demandent au premier ministre de declarer le Canada indivisible.
M. Jack Ramsay (Crowfoot, Ref.):
M. Hec Clouthier (Renfrew-Nippissing-Pembroke, Lib.):
Mme Michelle Dockrill (Bras d'Or, NPD):
Mme Rose-Marie Ur (Lambton-Kent-Middlesex, Lib.):
Monsieur le President, on repondra aujourd'hui a la question no 12. .[Texte]
M. Ted White (North Vancouver):
Monsieur le President, je suggere sur les autres questions soient reservees.
M. Deepak Obhrai (Calgary-Est, Ref.):
Merci, monsieur le President...
M. Scott Brison (Kings-Hants, PC):
Monsieur le President, j'invoque le Reglement.
Je pense qu'on s'engage dans un debat.
Ce n'est pas le moment de revenir la-dessus.
J'invite le depute a discuter de l'affaire avec le ministre.
Le depute de Calgary-Est a la parole.
Le whip en chef du gouvernement a la parole pour un rappel au Reglement.
M. Bob Kilger (Stormont-Dundas, Lib.):
La Chambre est-elle prete a se prononcer?
Plait-il a la Chambre d'adopter l'amendement?
Je donne la parole au whip en chef du gouvernement pour un rappel au Reglement.
M. Bob Kilger (Stormont-Dundas, Lib.):
Le vote porte sur la motion d'approbation du projet de loi C-7.
Plait-il a la Chambre d'adopter la motion?
La Chambre a entendu la proposition du whip en chef du gouvernement.
Quand lirons-nous ce projet de loi pour la troisieme fois?
Ils ont etudie le projet de loi C-7, article par article.
Tous les Canadiens peuvent etre de fiers gardiens de leur milieu marin.
Dans ce cas precis, il s'agit plus precisement du beluga.
Comme je l'ai deja mentionne, c'est prevu dans le projet de loi.
Les deputes savent que la gerance represente une partie importante du mandat de Parcs Canada.
Nous mettons sur pied un reseau de zones marines de conservation dans tout le Canada.
Ensemble, nous nous efforcerons d'assurer la perennite de ces ecosystemes marins critiques.
Je suis heureux d'un tel soutien.
Nous avons tous l'impression que nous venons d'un endroit special.
M. Deepak Obhrai (Calgary-Est, Ref.):
Le processus legislatif que nous amorcons aujourd'hui a la Chambre est important.
Toutes ses etapes sont eminemment importantes.
Les membres du comite ont ete unanimes a refuser la comparution de temoins.
Le comite n'a pas bien fait son travail.
Il n'a pas fait preuve de la diligence voulue.
Tous les autres partis ont approuve la decision sans rien dire.
Et voila maintenant que c'est au tour du comite.
Quelles sont les repercussions pour d'autres usagers du fleuve?
Quelles sont les repercussions en raison du chevauchement des reglements sur les parcs?
On ne peut pas creer des parcs isolement ou en vase clos.
Une consultation s'impose pour s'assurer qu'ils fonctionnent bien.
Cela ne s'est pas produit.
Les parcs constituent une question importante pour le Parti reformiste.
Mme Suzanne Tremblay (Rimouski-Mitis, BQ):
Le 5 juin 1997, l'Assemblee nationale du Quebec adoptait sa Loi 86.
Je suis heureuse de voir que nous assisterons bientot a la fin du processus.
Le president suppleant (M. McClelland):
La Chambre est-elle prete a se prononcer?
Le president suppleant (M. McClelland):
Le depute de Burin-St.
George's voudrait traiter du projet de loi.
Y a-t-il unanimite?
M. Bill Matthews (Burin-St.
Monsieur le President, je remercie les deputes de leur consentement.
Ce doit etre vendredi.
Felicitations pour voir fait de ce parc un succes.
Notre parti croit que ce projet fait beaucoup en ce sens.
Je suis heureux de declarer que nous continuerons de le faire.
Mme Jocelyne Girard-Bujold (Jonquiere, BQ):
Nous l'avons longtemps attendu, ce projet de loi.
Douze ans auparavant, l'idee d'un parc marin a ete lancee.
Mieux vaut tard que jamais, pourrions-nous dire a l'unisson.
Aussi, il nous faut parvenir a un equilibre.
Ce n'est pas peu dire, quatre milliards d'annees, ca donne de la perspective.
On est bien petits, finalement, dans ce bas monde.
Dans cette flore sous-marine, plusieurs especes de poissons en ont fait leur habitat.
M. Peter Stoffer (Sackville-Eastern Shore, NPD):
La preservation de cet ecosysteme marin avantagera tous les Canadiens.
On a garanti la participation des autochtones et accepte l'apport de tous les interesses.
Je ne reclame pas non plus un amendement.
Je remercie la Chambre de m'avoir donne la possibilite de parler de ces questions.
M. Guy Saint-Julien (Abitibi, Lib.):
Le 12 decembre 1996, les deux gouvernements ont depose leurs projets de loi respectifs.
Toutefois, l'Assemblee nationale du Quebec a sanctionne la Loi 86 le 5 juin 1997.
La loi provinciale entrera en vigueur lorsque la loi federale sera proclamee.
Le president suppleant (M. McClelland):
La Chambre est-elle prete a se prononcer?
Le president suppleant (M. McClelland):
Plait-il a la Chambre d'adopter la motion?
M. Gerry Byrne (secretaire parlementaire du ministre des Ressources naturelles, Lib.):
Monsieur le President, j'invoque le Reglement.
Le president suppleant (M. McClelland):
La Chambre a entendu la motion.
Le president suppleant (M. McClelland):
Mme Michelle Dockrill (Bras d'Or, NPD) propose:
Le charbon leur permet de toucher un salaire, qui, a son tour, stimule l'economie.
Plus de 73 000 hommes et femmes travaillent aujourd'hui dans l'industrie du charbon.
Le charbonnage n'est plus ce qu'il etait.
Aujourd'hui, l'exploitation miniere est tres mecanisee et extremement reglementee.
Une machine abat le travail que faisaient autrefois de nombreux hommes.
La technologie ne peut resoudre tous nos problemes, mais elle peut certainement aider.
Dans les annees 30, il y avait 400 petites mines de charbon.
Je puis dire avec fierte que la Devco emploie les meilleurs mineurs du monde.
Les mineurs du Cap-Breton comptent aujourd'hui sur leur societe d'Etat.
Ils veulent faire leur juste part.
Des milliers et des milliers d'emplois seraient perdus.
Je vais brosser un historique de Devco pour les deputes.
Se rendant compte du chaos que cela provoquerait, le gouvernement est intervenu.
En depit de cette situation, Devco est devenue plus efficiente et plus efficace.
Ce plan d'urgence, c'est la mine Donkin.
Donkin est l'avenir de Devco et l'entreprise demeure l'avenir du Cap-Breton.
La mine Donkin est un enorme gisement de charbon sous l'Atlantique.
La mine a deja ete construite aux frais des contribuables.
Les conditions geologiques semblent propices et la qualite du gisement semble meilleure encore.
Les activites pourraient etre accrues des que des profits seraient realises.
La mine Donkin est la mine d'or des mines de charbon.
Les contribuables canadiens ont deja investi 88 millions de dollars dans la mine Donkin.
On peut lire dans le rapport intitule The Donkin-Morien Mine:
Je repete que c'est la mine d'or des mines de charbon.
En temps de guerre comme en tant de crise...
Le president suppleant (M. McClelland):
M. Gerry Byrne (secretaire parlementaire du ministre des Ressources naturelles, Lib.):
Le gouvernement du Canada a fixe un but clair et raisonnable a Devco.
On ne peut insister trop sur l'obligation de reussir pour Devco.
Permettez-moi de decrire brievement certains des progres realises.
Les progres qu'ils ont accomplis jusqu'ici sont dignes d'eloges.
Il ne faut pas oublier pour autant ce qu'il reste a faire.
La Devco prevoit que sa tresorerie s'ameliorera d'ici l'exercice 1999-2000.
Ce n'est pas une mince affaire.
Cette mise a jour sera soumise au gouvernement au debut de 1998.
Le gouvernement s'engage a prendre en consideration toutes les options qui s'offriront.
M. Ghislain Fournier (Manicouagan, BQ):
Les richesses naturelles comme les mines, les forets et l'energie appartiennent aux provinces.
Parlons maintenant plus particulierement du dossier de Donkin.
Le federal n'y a jamais rien fait.
Prenez l'exemple de la Societe quebecoise d'exploration miniere.
Mais ce n'est pas le cas.
M. Peter MacKay (Pictou-Antigonish-Guysborough, PC):
Cette motion porte sur l'avenir economique de bon nombre de ces collectivites.
Il est consternant de penser que cela pourrait se produire.
Le gouvernement est tres imprevoyant s'il songe a prendre une telle mesure.
La situation semble evoluer plutot rapidement.
Que reste-t-il pour Devco?
Si ces deux mines ne peuvent etre exploitees, Devco est foutue.
Il ne lui restera rien.
Voila la realite des affaires economiques, mais le gouvernement a une responsabilite a assumer.
Il doit s'assurer que la situation de la Donkin soit traitee avec grand soin.
Cette responsabilite revient au gouvernement du Canada.
M. Darrel Stinson (Okanagan-Shuswap, Ref.):
Monsieur le President, nous etudions la motion:
Voyons ce qui se passe ici, dans le cas de Devco.
Le Senat est en train d'examiner la question.
Ce sont des questions qu'il faut poser.
La question est tres grave.
Les contribuables canadiens ont investi des centaines de millions de dollars dans ces projets.
Examinons la situation de facon raisonnable.
Pensez-vous qu'une societe privee se priverait de tous ses marches d'exportation?
Pourtant, que fait l'organisme gouvernemental?
Il se retire de son marche d'exportation et se cree un probleme.
Il n'a pas le sens des affaires.
Il n'a pas le sens des affaires et ne l'a jamais eu.
Je trouve plutot amusant que les conservateurs partagent mon avis.
On revient au meme point.
Nous savons que l'industrie miniere est un secteur d'activite extremement important au Canada.
Toutes les entreprises doivent etre rentables si elles veulent survivre.
C'est aussi simple que cela.
Ceux qui doivent soutenir la concurrence sur le marche mondial doivent bien comprendre la situation.
Je ne pense pas que je suis le seul a me poser cette question.
Les contribuables se la posent aussi.
M. Peter Stoffer (Sackville-Eastern Shore, NPD):
C'est vrai que les societes d'Etat comme la Devco ne sont pas efficaces.
Ce n'est pas la faute des travailleurs.
Le probleme, ce ne sont pas les societes d'Etat.
Le probleme, ce sont les gens qui veulent l'assiette au beurre.
Il merite les louanges du Nouveau Parti democratique et de tous les Canadiens.
L'autre jour, au Cap-Breton, j'ai rencontre trois personnes absolument formidables.
Les reformistes auraient considere qu'elle avait trois tort, et l'auraient condamnee sans appel.
J'ai peine a y croire.
Le gouvernement ne peut rien, nous ne pouvons vous aider.
J'enrage quand on aborde ce sujet.
Je voudrais bien le savoir avant l'ajournement de la seance aujourd'hui.
Le Parti liberal est-il favorable a cette motion?
Le president suppleant (M. McClelland):
M. Hec Clouthier (Renfrew-Nippissing-Pembroke, Lib.):
Pour ce faire, Devco doit devenir rentable.
Comme les deputes en ont ete informes, la societe continue de perdre de l'argent.
Cela se prononce Phalen.
Il y a desaccord sur la prononciation de Phalen.
Je viens de la vallee de l'Outaouais, et j'estime que c'est Phalen.
Vous etes en train de vous creuser un trou.
Nous sommes ici du cote du gouvernement, n'est-ce pas?
Oui, nous sommes du cote du leadership.
Nous nous attendons a ce que les projets entrepris soient commercialement viables.
Le mot vrai a une certaine qualite dans une institution aussi venerable que celle-ci.
C'est un mot tres polyvalent.
Les Canadiens reconnaissent que des progres ont ete realises a cet egard.
Le president suppleant (M. McClelland):
Y a-t-il consentement unanime?
Mme Michelle Dockrill (Bras d'Or, NPD):
Il a parle des organismes qui ont approuve ce financement, l'APECA et la SECB.
On ne le leur a pas demande.
Ou pourrait-on trouver ailleurs une telle situation?
Ou cela se rencontre-t-il, je le demande au depute.
C'est sous le gouvernement liberal actuel que cela se passe.
Il l'a fait sur le dos des habitants de l'ile du Cap-Breton.
Le president suppleant (M. McClelland):
(La seance est levee a 14 h 06.)

ext: f


text:
S.O.S. is the Ski or Snowboard Club at UC Davis.  SOS is the largest social Student Organizations club at UC Davis UCD and is devoted to Snow Sports snow sports.

>http://www.skiorsnowboardclub.com

http://www.facebook.com/ucdsos

http://www.theallcal.com<

>Image(alpinesos2.jpg)<

The UC Davis Ski or Snowboard Club, more commonly referred to as SOS, is the largest social club on campus. The purpose behind SOS is to get skiers and snowboarders up to the snow for as cheap as possible, and of course to have a great time in the process.  The Student Organizations club has ranged anywhere from 200800 members in past decades, though currently it hosts around 450 members along with 6 officers and 11 staffers.  Besides scoring its members huge savings, SOS is known for their wild parties and unique Drinking games drinking games.

Membership is $30. This fee scores you an SOS tshirt and a membership that will earn you tons of Living Cheaply great discounts around Davis town (Woodstocks Pizza Woodstocks, Kens Bike & Ski, Ground Zero, etc.).  Members are also able to attend all SOS events (All Cal, Cabin trips, Broomball, Float trip, etc...) for free or at a largely discounted price compared to nonmembers.  We host an annual Members Only sale night at Ground Zero, in which only members gain access to the store with items up to 80% off! (As well as many other useful discounts around Davis!) If you have a yearning to go to the snow but have no car to get yourself there, SOS also organizes carpools for all of our trips!  Being an SOS member also gives you full access to our Members only parties also known as Ragers.

REMEMBER!! We table in The Quad the quad Monday through Thursday from 11AM3PM throughout the year. Stop by, sign up for a trip, steal some stickers, or strike up a conversation with a friendly SOS officer or staff member. We can answer any questions you may have about the club, the social scene at Davis, where babies come from, the meaning of life, etc.

For more information or signups (we accept Pay Pal), check out the website http://www.skiorsnowboardclub.com or email our president at MailTo(sospres AT gmail DOT com)

 Upcoming Events

First Big Meeting, Freeborn Hall. (Date and Time TBA) Want to learn more about SOS? Come to First Big Meeting and well tell you everything you need to know. Weve also been know to throw out some free stuff at FBM.

ALL CAL @ Copper Mountain, Colorado. All Cal is the largest trip SOS runs throughout the entire year. Its a week long trip held during the first week of winter break along with UCSC, UCSB, UCI, UCSD, UCR and UCLA.  During this week about 500800 California college students will take over the town of Copper Mountain to ski, snowboard and party at one of the best ski resorts in the world. The price for the All Cal trip includes luxury accommodation, worryfree transportation, 4 days of riding as well as a plethora of parties and social events. Its a great opportunity to ride amazing snow as well as meet friends from all over California. Come to First Big Meeting for more information. Also, visit http://www.theallcal.com.

Cabin trip season! Every other weekend during Winter Quarter we rent out huge cabins in Tahoe and youre invited to join! Trip costs are super cheap, generally ranging from 4065 dollars for members and a bit more for nonmembers depending on the length of the trip. Trip costs include lodging and beverages. Cheap lift tickets are sold to particular mountains, usually around half price. Cabin trips are an excellent way to meet lots of cool people and have a great time with them. Trips fill up fast so when Cabin Trip season comes around, make sure to sign up ASAP at our table on the MU!

 Contacts

  President: Jonangelo JP, MailTo(sospres AT gmail DOT com)
  Vice President: Chelsea Truckee
  Treasurer: Mia
  Secretary: Alec
  Publicity: Jo

>Image(allcalbigsky.jpg)<

>Image(alpinesos.jpg)<

SOS Staff and Officers

20070221 23:59:39 nbsp I won a Forum Youngblood at the Big Meeting this year... SOS for the win. DanielSong

20090112 20:18:14 nbsp Dude, its time to update this page. Users/condemned2bfree

20100303 15:21:20 nbsp Im a Davis Enterprise reporter looking to head up to Northstar or Sierra this weekend (March 67) to capture local ski/snowboard feature story. If you think youve got something, ring me at 7478052 or email me at jedwards@davisenterprise.net Users/jon.edwards

max_stars_repo_name: voflo/Search

max_issues_repo_name: voflo/Search

max_forks_repo_name: voflo/Search


text:
C  *********************************************************
C  *                                                       *
C  *    TEST NUMBER: 09.01.04.01/03                        *
C  *    TEST TITLE : Error 152 handled by user             *
C  *                                                       *
C  *    PHIGS Validation Tests, produced by NIST           *
C  *                                                       *
C  *********************************************************

      CALL INITGL ('09.01.04.01/03')
      CALL E04152 (1)
      CALL ENDIT
      END


text:
C
C
*+ do_stretch

       subroutine do_stretch(map_array,status)
C      ---------------------------------------
C
C reproject a map to sky-plane
C
C Given:
C   map data
       real*4        map_array(*)
C Returned:
C   error status word
       integer       status
C
C The input map is reprojected to sky-plane.  This routine provides
C a subset of the facilities of reproject-map.
*-

C Local variables
       integer    minirt(8), ip_map, imap, ip_mapo, imapo
       integer    nx, ny

C check status on entry
       if (status.ne.0) return

C read input map
       call map_getmap('Map : ','Default_Map','read',imap,status)
       if (status.ne.0) goto 999
C read new range for map
       call enminirt(minirt,status)
       call plot_getuv('UV-range on input : ','*',minirt,status)
       if (status.ne.0) then
         call cmd_wrerr('STRETCH','Invalid range on input map')
         goto 999
       end if
C stetch the specified region of the map
       nx=minirt(2)-minirt(1)+1
       ny=minirt(3)-minirt(4)+1
C pointer to output map
       call map_alloc_out(0,0,'DIRECT',imapo,ip_mapo,status)
C force input map into core -- return pointer
       call map_alloc_in(imap,'DIRECT',map_array,ip_map,status)
       if (status.ne.0) goto 999

C extract the map data to the output map
       call map_stretch(map_array(ip_map),minirt,map_array(ip_mapo),
     *                  status)

C update the text of the redtape
       call adredt('CREATED','STRETCH',status)

C Tidy up, make output map new current map etc.
999    continue
       call map_end_alloc(imapo,map_array,status)
       call map_end_alloc(imap,map_array,status)

C check status value
       call cmd_err(status,'STRETCH','Failed ')
       end

max_stars_repo_name: CavendishAstrophysics/anmap

max_issues_repo_name: CavendishAstrophysics/anmap

max_forks_repo_name: CavendishAstrophysics/anmap


module ncfile_mod

use ncd_define_mod
use     ncaxis_mod, only: ncaxis_type, read_ncaxis, write_ncaxis, &
                                      close_ncaxis, copy_ncaxis
use      ncvar_mod, only:  ncvar_type, read_ncvar , write_ncvar, &
                           copy_ncvar
use      ncatt_mod, only:  ncatt_type, read_ncatt , write_ncatt, &
                            copy_ncatt
use   ncvarmap_mod, only: ncvarmap_type, copy_varmap

implicit none

!public :: read_ncfile, write_ncfile, close_ncfile
 public :: copy_ncfile

      type ncfile_type
           integer             :: id, ndim, nvar, natt, recdim, nvarmap
           character(len=maxn) :: name
           logical             :: define_mode
           type(  ncaxis_type), pointer ::   Axis(:)
           type(   ncvar_type), pointer ::    Var(:)
           type(   ncatt_type), pointer ::    Att(:)
           type(ncvarmap_type), pointer :: Varmap(:)
      end type ncfile_type

CONTAINS

!-----------------------------------------------------------------------

      function read_ncfile (name) result (File)

      character(len=*), intent(in) :: name
      type(ncfile_type)            :: File
      integer   ::  i, ierr, nv, nvar_total
      integer  ncopn

         File%name = name
         File%id = ncopn (name, NCNOWRIT, ierr)

         call ncinq (File%id, File%ndim, nvar_total, File%natt,  &
                     File%recdim, ierr)

	 if (File%ndim > 0) then
!del         allocate (File%Axis(max_ndim))
             allocate (File%Axis(File%ndim))
             do i=1,File%ndim
                if (i == File%recdim) then
                   File%Axis(i) = read_ncaxis (File%id, i, .true.)
                else
                   File%Axis(i) = read_ncaxis (File%id, i, .false.)
                endif
             enddo
         endif

         File%nvar = nvar_total - File%ndim
         if (File%nvar > 0) then
!del         allocate (File%Var(max_nvar))
             allocate (File%Var(File%nvar))
             do i=1,File%nvar
                File%Var(i) = read_ncvar (File%id, i+File%ndim)
             enddo
         endif

         if (File%natt > 0) then
             !allocate (File%Att(max_natt))
             allocate (File%Att(File%natt))
             do i=1,File%natt
                File%Att(i) = read_ncatt (File%id, NCGLOBAL, i)
             enddo
         endif

      end function read_ncfile

!-----------------------------------------------------------------------

      subroutine write_ncfile (File)

      type(ncfile_type) , intent(inout) :: File
      integer  i, nc, id, ierr, idum
      integer  nccre, ncsfil

         nc = len_trim(File%name)
!!!      File%id = nccre (File%name(1:nc), NCCLOB, ierr)
         id = nccre (File%name(1:nc), NCCLOB, ierr)
         File%id = id
         File%define_mode = .true.

         idum = ncsfil (File%id, NCNOFILL, ierr)

         do i=1,File%ndim
             call write_ncaxis (File%id, File%Axis(i))
         enddo

         do i=1,File%nvar
             call write_ncvar (File%id, File%Var(i))
         enddo

         do i=1,File%natt
             call write_ncatt (File%id, NCGLOBAL, File%Att(i))
         enddo

      end subroutine write_ncfile

!-----------------------------------------------------------------------

      function create_file (name) result (File)

      character(len=*), intent(in) :: name
      type(ncfile_type)            :: File
      integer  ierr, idum
      integer  nccre, ncsfil

           File%id = nccre (name, NCCLOB, ierr)
           idum = ncsfil (File%id, NCNOFILL, ierr)

           File%name    = name
           File%ndim    = 0
           File%nvar    = 0
           File%natt    = 0
           File%recdim  = 0
           File%nvarmap = 0

           File%define_mode = .true.

      end function create_file

!-----------------------------------------------------------------------

  subroutine close_ncfile (File)

   type(ncfile_type), intent(inout) :: File
   integer  ierr

!  ----- update length of record dimension  ------
      if (File%recdim > 0) then
        call close_ncaxis (File%id, File%recdim, File%Axis(File%recdim))
      endif

      call ncclos (File%id, ierr)

  end subroutine close_ncfile

!-----------------------------------------------------------------------

  subroutine empty_ncfile (File)

   type(ncfile_type), intent(inout) :: File

       if (associated(File%Axis  )) deallocate (File%Axis  )
       if (associated(File%Var   )) deallocate (File%Var   )
       if (associated(File%Att   )) deallocate (File%Att   )
       if (associated(File%Varmap)) deallocate (File%Varmap)

  end subroutine empty_ncfile

!-----------------------------------------------------------------------

      subroutine switch_define_mode (File)

      type(ncfile_type), intent(inout) :: File
      integer  ierr

        if (File%define_mode) then
            call ncendf (File%id, ierr)
            File%define_mode = .false.
        else
            call ncredf (File%id, ierr)
            File%define_mode = .true.
        endif

      end subroutine switch_define_mode

!-----------------------------------------------------------------------

   subroutine copy_ncfile (File_in, File_out)

     type(ncfile_type), intent(in)  :: File_in
     type(ncfile_type), intent(out) :: File_out

     integer :: i

!     ---- copy axes ----
      File_out%ndim = File_in%ndim
      allocate (File_out%Axis(File_out%ndim))
      do i = 1, File_out%ndim
         call copy_ncaxis ( File_in%Axis(i), File_out%Axis(i) )
      enddo

!     ---- copy variables ----
      File_out%nvar = File_in%nvar
      allocate (File_out%Var(File_out%nvar))
      do i = 1, File_out%nvar
         call copy_ncvar ( File_in%Var(i), File_out%Var(i) )
      enddo

!     ---- copy global attributes ----
      File_out%natt = File_in%natt
      allocate (File_out%Att(File_out%natt))
      do i = 1, File_out%natt
         call copy_ncatt ( File_in%Att(i), File_out%Att(i) )
      enddo

!     ---- copy variable mapping ----
      File_out%nvarmap = File_in%nvarmap
      allocate (File_out%Varmap(File_out%nvarmap))
      do i = 1, File_out%nvarmap
         call copy_varmap ( File_in%Varmap(i), File_out%Varmap(i) )
      enddo

!     ---- miscellaneous ----
      File_out%recdim = File_in%recdim
      File_out%define_mode = .true.

!     ---- un-copied data ----
!     File_out%id   =
!     File_out%name =

   end subroutine copy_ncfile

!-----------------------------------------------------------------------

end module ncfile_mod


text:
C Copyright (c) 2007-2017 National Technology & Engineering Solutions of
C Sandia, LLC (NTESS).  Under the terms of Contract DE-NA0003525 with
C NTESS, the U.S. Government retains certain rights in this software.
C
C Redistribution and use in source and binary forms, with or without
C modification, are permitted provided that the following conditions are
C met:
C
C     * Redistributions of source code must retain the above copyright
C       notice, this list of conditions and the following disclaimer.
C
C     * Redistributions in binary form must reproduce the above
C       copyright notice, this list of conditions and the following
C       disclaimer in the documentation and/or other materials provided
C       with the distribution.
C
C     * Neither the name of NTESS nor the names of its
C       contributors may be used to endorse or promote products derived
C       from this software without specific prior written permission.
C
C THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
C "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
C LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
C A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
C OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
C SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
C LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
C DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
C THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
C (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
C OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
C

C=======================================================================
*DECK,SINTPE
      SUBROUTINE SINTPE (ICONA,SOLENA,ISRCHR,NISR,RSRCHR,NRSR,
     &                   SOLEB,IDBLK,XB,YB,ZB,
     &                   ICONB,ITT,IBLK,TIMES,CENTER,
     &                   ISTP,IST,INSUB,ICOMPL,DUME)
C
C     ******************************************************************
C
C     SUBROUTINE TO CONTROL INTERPOLATION OF ELEMENT TRANSFORMED
C     INTO NODAL RESULTS FROM MESH-A TO MESH-B
C     INTERPOLATED SOLUTION IS PASSED OUT TO BE RETRANSFORMED
C     INTO ELEMENT RESULTS AND THEN WRITTEN TO MESH-C EXODUS FILE
C
C     Calls subroutines SHAPEF
C
C     Called by MAPVAR
C
C     ******************************************************************
C
C ICONA   INT   Connectivity of donor mesh (1:nelnda,1:numeba)
C SOLENA  REAL  Element variables at nodes for donor mesh
C ISRCHR  INT   Donor mesh element in which point is found
C NISR    INT   Dimension for ISRCHR
C RSRCHR  REAL  Isoparametric coordinates of point in ISRCHR
C NRSR    INT   Dimension for RSRCHR
C SOLEB   REAL  Element variables for RECIPIENT Mesh
C IDBLK   INT   Recipient mesh element block I.D.
C XB,etc  REAL  Recipient mesh coordinates
C ICONB   INT   Recipient mesh connectivity
C ITT     INT   Truth table
C IM      INT   Mesh block being processed (not block I.D.)
C TIMES   REAL  Array of times associated with time step
C CENTER  REAL  Recipient mesh element centroid coordinates
C ISTP    INT   Time step
C INSUB   INT   Number of times into subroutine for this element block
C               1-first time; >1-second,etc time;
C               used to control mapping for many element blocks to one
C ICOMPL  INT   Flag for completion of mapping; used in many to one
C               0-incomplete; 1-complete
C
C     ******************************************************************
C
      include 'aexds1.blk'
      include 'aexds2.blk'
      include 'amesh.blk'
      include 'bmesh.blk'
      include 'ebbyeb.blk'
      include 'ex2tp.blk'
      include 'tapes.blk'
C
      DIMENSION TIMES(*), CENTER(NUMEBB,*)
      DIMENSION ICONA(NELNDA,*), SOLENA(NODESA,NVAREL)
      DIMENSION ITT(NVAREL,*)
      DIMENSION SOLEB(NUMEBB,*), XB(*), YB(*), ZB(*)
      DIMENSION ISRCHR(NISR,*), RSRCHR(NRSR,*), ICONB(nelndb,*)
      DIMENSION SOLN(27), XX(27), YY(27), ZZ(27), DUME(*)
C
C     ******************************************************************
C
      DO 10 IVAR = 1,NVAREL
        IF (ITT(IVAR,IBLK) .EQ. 0)GO TO 10
C
C Initialize SOLEB if first time in subroutine for this element block
C After first time into subroutine
C retrieve SOLEB from storage in EXODUS
C
        IF (INSUB .EQ. 1)THEN
          CALL INIELT(SOLEB,IVAR,TIMES,ISTP,IDBLK,CENTER,DUME)
        ELSE
          CALL EXGEV(NTP4EX,IST,IVAR,IDBLK,NUMEBB,SOLEB(1,IVAR),IERR)
        END IF
C
C Loop on centroids in recipient mesh
C
        DO 30 I = 1,NUMEBB
          IF (ISRCHR(1,I) .NE. 0)THEN
C
C Set parameters for element in donor mesh
C
            S = RSRCHR(5,I)
            T = RSRCHR(6,I)
            R = 0.
            NNODES = 4
            DO 20 J = 1,NNODES
              INODE = ICONA(J,ISRCHR(1,I))
              SOLN(J) = SOLENA(INODE,IVAR)
 20         CONTINUE
C
C Shape function to evaluate interpolation
C
            CALL SHAPEF (3,S,T,R,SOLN,BVALUE)
            SOLEB(I,IVAR) = BVALUE
          ELSE
            if (ivar .eq. 1 .and. icompl .eq. 1) then
              write (nout, 9000) i, idblk
              write (ntpout, 9000) i, idblk
 9000         FORMAT('SINTPE: No donor element found for ',
     *          'recipient Element', i9, ' in block ', i9)
            end if
          END IF
 30     CONTINUE
C
C  If there is more searching to do (i.e. many blocks to one)
C  use EXODUS as temporary storage
C  don't bother to make adjustments to element variables yet
C
        IF (ICOMPL .NE. 1)THEN
          CALL EXPEV(NTP4EX,IST,IVAR,IDBLK,NUMEBB,SOLEB(1,IVAR),IERR)
        ELSE
C
C Do needed modification to element variable values
C write element vars out to EXODUS data base
C
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C ELMASS is special
C
          IF (NAMVAR(nvargp+IVAR)(1:6) .EQ. 'ELMASS')THEN
C
C ELMASS was changed to nodal density prior to processing.
C need to go back from density to element mass now
C
            DO 100 IEL = 1, NUMEBB
              DO 105 I = 1, NNODES
                XX(I) = XB(ICONB(I,IEL))
                YY(I) = YB(ICONB(I,IEL))
                ZZ(I) = ZB(ICONB(I,IEL))
 105          CONTINUE
              CALL VOL(ITYPE,XX,YY,ZZ,VOLUME)
              SOLEB(IEL,IVAR) = SOLEB(IEL,IVAR) * VOLUME
 100        CONTINUE
          END IF
C
C****************************************************************
C apply constraints to variables here as applicable
C
C Plastic strain (EQPS) and tearing
C must greater than or equal to 0.
C
          IF (NAMVAR(nvargp+IVAR)(1:4) .EQ. 'EQPS')THEN
            DO 110 IEL = 1, NUMEBB
              IF (SOLEB(IEL,IVAR) .LT. 0.)THEN
                SOLEB(IEL,IVAR) = 0.
              END IF
  110       CONTINUE
          ELSE IF (NAMVAR(nvargp+IVAR)(1:7) .EQ. 'TEARING')THEN
            DO 120 IEL = 1, NUMEBB
              IF (SOLEB(IEL,IVAR) .LT. 0.)THEN
                SOLEB(IEL,IVAR) = 0.
              END IF
  120       CONTINUE
C
C Hourglass forces and bulk viscosity have no meaning other than on
C the mesh from which they originated, just set them to zero.
C
          ELSE IF (NAMVAR(nvargp+IVAR)(1:2) .EQ. 'HG' .OR.
     &            NAMVAR(nvargp+IVAR)(1:5) .EQ. 'BULKQ')THEN
            DO 130 IEL = 1, NUMEBB
              SOLEB(IEL,IVAR) = 0.
  130       CONTINUE
C
C Shell element basis vectors are computed from the mesh
C It makes no sense to transfer them - set them to zero.
C
          ELSE IF (NAMVAR(nvargp+IVAR)(1:5) .EQ. 'BASEL')THEN
            DO 140 IEL = 1, NUMEBB
              SOLEB(IEL,IVAR) = 0.
 140        CONTINUE
C
C         ELSE IF (NAMVAR(nvargp+IVAR)(1:?) .EQ. ?????)THEN
C           DO ??? IEL = 1, NUMEBB
C             IF (SOLEB(IEL,IVAR) .??. ?.)THEN
C               SOLEB(IEL,IVAR) = ?.
C             END IF
C  ???      CONTINUE
C
          END IF
c************************************************************************
C
C write element variables
C
          CALL EXPEV(NTP4EX,IST,IVAR,IDBLK,NUMEBB,SOLEB(1,IVAR),IERR)
        END IF
   10 CONTINUE
      RETURN
      END

max_stars_repo_name: mathstuf/seacas

max_stars_repo_licenses: BSD-3-Clause NetCDF Zlib MIT

max_issues_repo_name: mathstuf/seacas

max_issues_repo_licenses: BSD-3-Clause NetCDF Zlib MIT

max_forks_repo_name: mathstuf/seacas

max_forks_repo_licenses: BSD-3-Clause NetCDF Zlib MIT


text: Jason Dayne is a local artists local artist. who produces handcrafted silver art jewelry available via his website.  His workshop is located downtown.


program arrayops
    implicit none
    integer :: i, j
    integer, parameter :: M = 5
    real, dimension(M) :: x
    real :: s, t

    ! initialize the array
    do i=1, M
        x(i) = i
    end do

    ! sum
    s = 0.
    do i=1,M
        s = s + x(i)
    end do
    write(*,*)"sum = ",s

    ! average (mean)
    write(*,*)"average = ",s/M

    ! product
    s = 1.
    do i=1,M
        s = s * x(i)
    end do
    write(*,*)"product = ",s

    ! maximum value
    s = x(1)
    do i=2,M
        if (x(i) > s) then
            s = x(i)
        end if
    end do
    write(*,*)"max = ",s

    ! minimum value
    s = x(1)
    do i=2,M
        if (x(i) < s) then
            s = x(i)
        end if
    end do
    write(*,*)"min = ",s

    ! norm
    s = 0.
    do i=1,M
        s = s + x(i) * x(i)
    end do
    write(*,*)"norm = ",sqrt(s)

    ! second largest value
    if (x(1) > x(2)) then
        s = x(2)
        t = x(1)
    else
        s = x(2)
        t = x(1)
    end if
    do i=3,M
        if (x(i) > s) then
            if (x(i) > t) then
               s = t
               t = x(i)
            else
                s =x(i)
            end if
        end if
    end do
    write(*,*)"second largest value = ", s

end program


text:
      SUBROUTINE CALHEAT(ISTL_)
C
C   Subroutine CALHEAT takes the information from the atmospheric boundary
C   file and the wind forcing file and calculates the net heat flux across
C   the water surface boundary.  The heat flux is then used to update the
C   water temperature either in the surface cells, or distributed across
C   the cells in the vertical and into the bottom.  The subroutine has
C   three options these are:
C
C    ISOPT(2)=1: Full surface and internal heat transfer calculation
C                using meteorologic data from input stream.
C                IASWRAD=0:  ADSORB SW SOLR RAD TO ALL LAYERS AND BED
C                IASWRAD=1:  ADSORB SW SOLR RAD TO TO SURFACE LAYER
C    ISOPT(2)=2: Transient equilibrium surface heat transfer calculation
C                using external equilibrium temperature and heat transfer
C                coefficient data from the meteorologic input data.
C    ISOPT(2)=3: Equilibrium surface heat transfer calculation using constant
C                equilibrium temperature and heat transfer coefficients
C                HEQT and TEMO read in through input stream.
C    ISOPT(2)=4: Equilibrium surface heat transfer calculation using algorithm
C                from CE-QUAL-W2.
C
C   The heat flux terms are derived from a paper by Rosati
C   and Miyakoda (1988) entitled "A General Circulation Model for Upper Ocean
C   Simulation".  The heat flux is prescribed by term for the following
C   influxes and outfluxes:
C
C     - Short Wave Incoming Radiation (+)
C     - Net Long Wave Radiation (+/-)
C     - Sensible Heat Flux (convection -)
C     - Latent Heat Flux (evaporation +/-)
C
C   Two formulations of the Latent Heat Flux are provided.  The first is from
C   the Rosati and Miyakoda paper, the second is an alternate formulation by
C   Brady, Graves, and Geyer (1969).  The second formulation was taken from
C   "Hydrodynamics and Transport for Water Quality Modeling" (Martin and
C   McCutcheon, 1999).  The Rosati and Miyakoda formulation will have zero
C   evaporative cooling or heating if wind speed goes to zero.  The Brady,
C   Graves, and Geyer formulation provides for a minimum evaporative cooling
C   under zero wind speed.
C
C
C VARIABLE LIST:
C
C   CLOUDT  = Cloud cover (0 to 10)
C   HCON    = Sensible heat flux (W/m2)
C   HLAT    = Latent heat flux (W/m2)
C   HLWBR   = Net longwave radiation (atmospheric long wave plus back
C             radiation, W/m2)
C   SOLSWRT = Short wave incoming radiation (W/m2)
C   SVPW    = Saturation vapor pressure in mb based upon the water surface
C             temperature
C   TATMT   = Temperature of air above water surface (deg C)
C   TEM     = Water temperature in cell (deg C)
C   VPA     = Vapor pressure of air at near surface air temperature (mb)
C   WINDST  = Wind speed at 10 meters over cell surface (m/s)
C
C MODIFICATION HISTORY:
C
C   Date       Author             Comments
C   ---------- ------------------ ---------------------------------------------
C   06/01/1992 John M. Hamrick    Orignial author
C   11/07/2000 Steven Peene       Cleaned code, provided more detailed
C                                 descriptions, added alternate formulation
C                                 for latent heat flux, separated out
C                                 individual heat flux terms
C   11/01/2005 Paul M. Craig      Added Option 4, to use the Equilibrium Temperature
C                                 algorithym from CE-QUAL-W2.  Also added the sub-option
C                                 under this option to couple or decouple the bottom temperature
C                                 to the water column temperatures.
C                                 Added the ability to input spatially variable bed temps and
C                                 thermally active bed thicknesses.
C                                 Also cleaned up the code and added structure.
C
C CHANGE RECORD
C **  SUBROUTINE CALHEAT CALCULATES SURFACE AND INTERNAL HEAT SOURCES
C **  AND SINKS IN THE HEAT (TEM) TRANSPORT EQUATION
C
      USE GLOBAL
      USE OMP_LIB
      REAL,SAVE,ALLOCATABLE,DIMENSION(:,:)::NETRAD
      REAL,SAVE,ALLOCATABLE,DIMENSION(:)::TBEDTHK
      REAL,SAVE,ALLOCATABLE,DIMENSION(:)::HDEP
      REAL,SAVE,ALLOCATABLE,DIMENSION(:)::RADBOT
      REAL,SAVE,ALLOCATABLE,DIMENSION(:)::FLUXTB
      !REAL,SAVE ::    PTIME
      !REAL,SAVE ::    PMCTOL
      REAL            K_ABOVE
C
      IF(.NOT.ALLOCATED(NETRAD))THEN
        ALLOCATE(NETRAD(LCM,KCM))
        ALLOCATE(TBEDTHK(LCM))
        ALLOCATE(HDEP(LCM))
        ALLOCATE(RADBOT(LCM))
        ALLOCATE(FLUXTB(LCM))
	  RADBOT=0.0 !SCJ
	  FLUXTB=0.0 !SCJ

        ! *** Ininitialze Heat Exchange Terms
        IF (PARTID ==0) PRINT *,'CALHEAT: INITIALIZING'
        CALL HEAT_EXCHANGE

        !PMCTOL=0.1
        NETRAD=0.
        HDEP=0.
        TBEDTHK=0.
        IF((ISTOPT(2).EQ.1.AND.IASWRAD.EQ.0).OR.ISTOPT(2).EQ.4)THEN
          IF(DABEDT.GT.0.)THEN
            IF(PARTID==MASTER_TASK) PRINT *,'CALHEAT: SET CONSTANT THICKNESS TO:',DABEDT
            DO L=2,LA
              TBEDTHK(L)=DABEDT
            ENDDO
          ELSE
            ! *** READ IN THE SPATIALLY VARYING INIT T AND BED THICKNESS (DABEDT)
            IF (PARTID==MASTER_TASK) PRINT *,'CALHEAT: READ IN SPATIALLY VARYING INIT T AND
     & BED THICKNESS: TEMB.INP'
            DO L=2,LA
              TBEDTHK(L)=ABS(DABEDT)
              IF(ISCI(2).EQ.0)TEMB(L)=ABS(TBEDIT)
            ENDDO
            OPEN(1001,FILE='TEMB.INP',ERR=1000,STATUS='OLD')
            DO IS=1,4
              READ(1001,*)
            ENDDO
            DO L1 = 2, LA
              READ(1001,*,END=1000) I,J,T1,T2
              L=LIJ(I,J)
              IF(ISCI(2).EQ.0)TEMB(L)=T1
              TBEDTHK(L)=T2
            ENDDO
 1000       CLOSE(1001)
          ENDIF
        ENDIF

        IF(ISTOPT(2).EQ.1.OR.ISTOPT(2).EQ.4)THEN
          DO L=2,LA
            PSHADE(L)=1.0
          ENDDO
          IF(USESHADE)THEN
            ! *** READ IN THE SPATIALLY VARYING INIT T AND BED THICKNESS (DABEDT)
            IF (PARTID == 0) PRINT *,
     *        'CALHEAT: READ IN SPATIALLY VARYING SHADE: PSHADE.INP'
            OPEN(1001,FILE='PSHADE.INP',ERR=1010,STATUS='OLD')
            DO IS=1,4
              READ(1001,*)
            ENDDO
            DO L1 = 2, LA
              READ(1001,*,END=1010) I,J,T1
              L=LIJ(I,J)
              PSHADE(L)=T1
            ENDDO
 1010       CLOSE(1001)
          ELSE
            IF (PARTID ==0) PRINT *,
     &      'CALHEAT: SET CONST SHADE TO: 1.0 (NO SHADE)'
          ENDIF
        ENDIF

        IF(DEBUG)THEN
          PRINT *,'CALHEAT: Bed Temp(L=2):', TEMB(2)
          OPEN(77,file='calheat.dia',status='unknown')
          CLOSE(77,status='DELETE')
          OPEN(77,file='calheat.dia',status='NEW')
          WRITE(77,998)'TIMEDAY','SRON','ET','TD_C','TA_C','TDEW_F',
     &                  'TAIR_F','FW'
  998     FORMAT(A11,8A9)
        ENDIF
      ENDIF
C
      ! *** DSLLC BEGIN CHANGE
CPMC      DELT=DT2
CPMC      IF(ISTL_.EQ.2)THEN
CPMC        DELT=DT
CPMC      ENDIF
CPMC      DELT=DT2
      IF(ISTL_.EQ.2)THEN
        IF(ISDYNSTP.EQ.0)THEN
          DELT=DT
        ELSE
          DELT=DTDYN
        END IF
      ELSE
        DELT=DT2
      ENDIF
      ! *** DSLLC END CHANGE

      ! *** OVERWRITE THE INPUT SOLAR RAD WITH A COMPUTED ONE
      IF(COMPUTESOLRAD)THEN
        CALL SHORT_WAVE_RADIATION(WINDST(2),RHA(2),TATMT(2),CLOUDT(2),
     &                            PATMT(2),SRO,SRON)
        DO L=2,LA
          ! *** USE COMPUTED SRO
          SOLSWRT(L)=SRON
        ENDDO
      ENDIF

      IF(USESHADE)THEN
        DO L=2,LA
          ! *** APPLY PSHADE FACTORS
          SOLSWRT(L)=SOLSWRT(L)*PSHADE(L)
        ENDDO
      ENDIF

      IF(ISTOPT(2).EQ.1)THEN
        ! *** FULL HEAT BALANCE WITH ATMOSPHERIC LINKAGE

        ! *** SET UP MIN DEPTH
        DO L=2,LA
          HDEP(L)=MAX(HP(L),0.)
        ENDDO

        ! *** NET SHORTWAVE SOLAR RADIATION
        IF(IASWRAD.EQ.0.)THEN
          ! *** ADSORB SW SOLR RAD TO ALL LAYERS AND BED

        ! NET HEAT FLUX = RSN+RAN-RB-RE-RC  (WATT/M2)
!$OMP PARALLEL PRIVATE(SVPW,CLDFAC,RAN,FW,RE,RC,RB,RSN,
!$OMP& UBED,VBED,USPD,TMPVAL,TFAST,TFAST1,C1,TSLOW,TSLOW1)

          ! *** SURFACE LAYER
          TFAST=SWRATNF*(Z(KC)-1.)
          TFAST1=SWRATNF*(Z(KC-1)-1.)
          TSLOW=SWRATNS*(Z(KC)-1.)
          TSLOW1=SWRATNS*(Z(KC-1)-1.)
        IF (NASER == -1 ) THEN ! USE DeepThunder provided values of longwave radiation
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
          DO L=2,LA
            SVPW=(10.**((0.7859+0.03477*TEM(L,KC))/
     &          (1.+0.00412*TEM(L,KC))))

           ! *** Evaporation
            FW=9.2+0.46*WINDST(L)**2;
           !Saturated water vapor pres - actual water vapor press (eqn 2.3.10 in Ji 2008)
            RE=FW*(SVPW-VPA(L));
            ! *** Conduction
            RC=0.47*FW*(TEM(L,KC)-TATMT(L))
            NETRAD(L,KC)=LWDNB(L) - LWUPB(L) -RE-RC
          ENDDO
!$OMP END DO NOWAIT
        ELSE   ! use existing implementation based on Swinbank (1963) formula
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
          DO L=2,LA
            SVPW=(10.**((0.7859+0.03477*TEM(L,KC))/
     &          (1.+0.00412*TEM(L,KC))))

            ! *** Net atmospheric radiation (Diffusive)
            CLDFAC=1.0+0.0017*CLOUDT(L)**2
            RAN=0.51531538e-12*(TATMT(L)+273.15)**6*CLDFAC
            ! *** Evaporation
            FW=9.2+0.46*WINDST(L)**2;
           !Saturated water vapor pres - actual water vapor press (eqn 2.3.10 in Ji 2008)
            RE=FW*(SVPW-VPA(L));
            ! *** Conduction
            RC=0.47*FW*(TEM(L,KC)-TATMT(L))
            ! *** Longwave back radiation
            ! ***  5.443E-8 = 5.67E-8 * 0.97
            RB=5.443E-8*(TEM(L,KC)+273.15)**4 ! upward radiation flux
            NETRAD(L,KC)=RAN-RB-RE-RC
          ENDDO
!$OMP END DO NOWAIT
        END IF
          IF(FSWRATF.LT.1.)THEN
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
            DO L=2,LA
              RSN=SOLSWRT(L)*
     &           (    FSWRATF*(EXP(TFAST*HDEP(L))-EXP(TFAST1*HDEP(L)))
     &          +(1.-FSWRATF)*(EXP(TSLOW*HDEP(L))-EXP(TSLOW1*HDEP(L))))
              NETRAD(L,KC)=NETRAD(L,KC)+RSN
            ENDDO
!$OMP END DO NOWAIT
          ELSE
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
            DO L=2,LA
              RSN=SOLSWRT(L)*(1.-EXP(TFAST1*HDEP(L)))
              NETRAD(L,KC)=NETRAD(L,KC)+RSN
            ENDDO
!$OMP END DO NOWAIT
          ENDIF

          ! *** ALL REMAINING LAYERS
          IF(KC.GT.1)THEN
            DO K=1,KS
              TFAST=SWRATNF*(Z(K)-1.)
              TFAST1=SWRATNF*(Z(K-1)-1.)
              C2=DELT*DZIC(K)*0.2393E-6
              IF(FSWRATF.LT.1.)THEN
                TSLOW=SWRATNS*(Z(K)-1.)
                TSLOW1=SWRATNS*(Z(K-1)-1.)
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
                DO L=2,LA
                  RSN=SOLSWRT(L)*
     &             (    FSWRATF*(EXP(TFAST*HDEP(L))-EXP(TFAST1*HDEP(L)))
     &           +(1.-FSWRATF)*(EXP(TSLOW*HDEP(L))-EXP(TSLOW1*HDEP(L))))
                  NETRAD(L,K)=RSN
                ENDDO
!$OMP END DO NOWAIT
              ELSE
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
                DO L=2,LA
                  RSN=SOLSWRT(L)*
     &               (EXP(TFAST*HDEP(L))-EXP(TFAST1*HDEP(L)))
                  NETRAD(L,K)=RSN
                ENDDO
!$OMP END DO NOWAIT
              ENDIF
            ENDDO
          ENDIF

          ! *** Distribute heat flux to the bed for each grid cell.
          TFAST=SWRATNF*(Z(0)-1.)
          IF(FSWRATF.LT.1.)THEN
            TSLOW=SWRATNS*(Z(0)-1.)
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
            DO L=2,LA
              LE=LEAST(L)
              LN=LNC(L)
              UBED=0.5*( U(L,1)+U(LE,1) )
              VBED=0.5*( V(L,1)+V(LN,1) )
              USPD=SQRT( UBED*UBED+VBED*VBED )
              TMPVAL=(HTBED1*USPD+HTBED2)*(TEM(L,1)-TEMB(L))
              NETRAD(L,1)=NETRAD(L,1)-TMPVAL/0.2393E-6

              ! *** UPDATE BOTTOM
              IF(TBEDIT.GT.0.)THEN
                TEMB(L)=TEMB(L) + (TMPVAL + 0.2393E-6*SOLSWRT(L)*
     &               (FSWRATF *EXP(TFAST*HDEP(L))
     &           +(1.-FSWRATF)*EXP(TSLOW*HDEP(L))))*DELT/TBEDTHK(L)
              ENDIF
            ENDDO
!$OMP END DO NOWAIT
          ELSE
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
            DO L=2,LA
              LE=LEAST(L)
              LN=LNC(L)
              UBED=0.5*( U(L,1)+U(LE,1) )
              VBED=0.5*( V(L,1)+V(LN,1) )
              USPD=SQRT( UBED*UBED+VBED*VBED )
              TMPVAL=(HTBED1*USPD+HTBED2)*(TEM(L,1)-TEMB(L))
              NETRAD(L,1)=NETRAD(L,1)-TMPVAL/0.2393E-6

              ! *** UPDATE BOTTOM
              IF(TBEDIT.GT.0.)THEN
                TEMB(L)=TEMB(L) + (TMPVAL + 0.2393E-6*SOLSWRT(L)*
     &                  FSWRATF*EXP(TFAST*HDEP(L)))*DELT/TBEDTHK(L)
              ENDIF
            ENDDO
!$OMP END DO NOWAIT
          ENDIF

          ! *** NOW FINALIZE THE TEMPERATURE
          DO K=1,KC
            ! *** RHO = 1000.0  Density (kg / m^3)
            ! *** CP  = 4179.0  Specific Heat (J / kg / degC)
            ! *** 0.2393E-6 = 1/RHO/CP
            C1=DELT*DZIC(K)*0.2393E-6
!$OMP DO SCHEDULE(STATIC,CHUNKSIZE)
            DO L=2,LA
              TEM(L,K)=TEM(L,K)+HPI(L)*C1*NETRAD(L,K)
              IF (TEM(L,K) < 0.85) TEM(L,K) = 0.85   ! include this hard code until development of full ice-model
            ENDDO
!$OMP END DO NOWAIT
          ENDDO
!$OMP END PARALLEL

          IF(ISDRY.GT.0.AND.ISTOPT(2).EQ.1)THEN
            DO L=2,LA
              IF(IMASKDRY(L).EQ.1.) TEMB(L)=TATMT(L)
            ENDDO
          ENDIF

        ELSE   ! IF(IASWRAD.EQ.1)THEN

          ! NET HEAT FLUX = RSN+RAN-RB-RE-RC  (WATT/M2)
          C1=DELT*DZIC(KC)*0.2393E-6
          IF (NASER == -1) THEN  ! Deep Thunder computed longwave radiation
            DO L=2,LA
              SVPW=(10.**((0.7859+0.03477*TEM(L,KC))/
     &            (1.+0.00412*TEM(L,KC))))

              ! *** Evaporation
               FW=9.2+0.46*WINDST(L)**2;
               RE=FW*(SVPW-VPA(L));     !VPA fcn(TEMP,AirTemp,RelativeHumidity)

               ! *** Conduction
               RC=0.47*FW*(TEM(L,KC)-TATMT(L))

               NETRAD(L,KC)=LWDNB(L)-LWUPB(L)-RE-RC + SOLSWRT(L)    ! *** ADSORB SW SOLR RAD TO TO SURFACE LAYER
               ! *** NOW FINALIZE THE TEMPERATURE
               TEM(L,KC)=TEM(L,KC)+HPI(L)*C1*NETRAD(L,KC)
             ENDDO
           ELSE   ! Use existing implementationfor LW rad  based on Swinbank (1963)
             DO L=2,LA
               SVPW=(10.**((0.7859+0.03477*TEM(L,KC))/
     &            (1.+0.00412*TEM(L,KC))))

               ! *** Net atmospheric radiation (Diffusive)
               CLDFAC=1.0+0.0017*CLOUDT(L)**2
               ! ** .5153153831e-12 = 1000.0/3600.0*9.37E-6*2.0411E-7*0.97
               RAN=0.51531538e-12*(TATMT(L)+273.15)**6*CLDFAC

               ! *** Evaporation
               FW=9.2+0.46*WINDST(L)**2;
               RE=FW*(SVPW-VPA(L));     !VPA fcn(TEMP,AirTemp,RelativeHumidity)

               ! *** Conduction
               RC=0.47*FW*(TEM(L,KC)-TATMT(L))

               ! *** Longwave back radiation
               ! ***  5.443E-8 = 5.67E-8 * 0.97
               RB=5.443E-8*(TEM(L,KC)+273.15)**4

               NETRAD(L,KC)=RAN-RB-RE-RC + SOLSWRT(L)    ! *** ADSORB SW SOLR RAD TO TO SURFACE LAYER
               ! *** NOW FINALIZE THE TEMPERATURE
               TEM(L,KC)=TEM(L,KC)+HPI(L)*C1*NETRAD(L,KC)
             ENDDO
           END IF
       ENDIF

      ELSEIF(ISTOPT(2).EQ.2)THEN
        ! *** IMPLEMENT EXTERNALLY SPECIFIED EQUILIBRIUM TEMPERATURE FORMULATION
        TMPKC=DELT/DZC(KC)
        DO ND=1,NDM
          LF=2+(ND-1)*LDM
          LL=LF+LDM-1
          DO L=LF,LL
            TEM(L,KC)=TEM(L,KC)-TMPKC*CLOUDT(L)*HPI(L)*(TEM(L,KC)
     &          -TATMT(L))
!            TEM(L,KC)=TEM(L,KC)-TMPKC*HEQT*HPI(L)*(TEM(L,KC)
!     &          -TATMT(L))  !SCJ changed
          ENDDO
        ENDDO

      ELSEIF(ISTOPT(2).EQ.3)THEN
        ! *** IMPLEMENT CONSTANT COEFFICIENT EQUILIBRIUM TEMPERATURE FORMULATION
        DTHEQT=DELT*HEQT*FLOAT(KC)
!        DO ND=1,NDM
!          LF=2+(ND-1)*LDM
!          LL=LF+LDM-1
!          DO L=LF,LL
!            TEM(L,KC)=TEM(L,KC)-DTHEQT*HPI(L)*(TEM(L,KC)-TEMO)
!          ENDDO
!        ENDDO
        DO L=2,LA
          TEM(L,KC)=TEM(L,KC)-DTHEQT*HPI(L)*(TEM(L,KC)-TEMO)
        ENDDO

      ELSEIF(ISTOPT(2).EQ.4)THEN
        ! *** IMPLEMENT W2 EQUILIBRIUM TEMPERATURE FORMULATION
        IF(.NOT.COMPUTESOLRAD)THEN
          ! *** MUST MAKE AT LEAST ONE CALL TO THIS TO INITIALIZE VARIABLES
          CALL SHORT_WAVE_RADIATION(WINDST(2),RHA(2),TATMT(2),CLOUDT(2),
     &                            PATMT(2),SRO,SRON)
        ENDIF
        ! *** SWRATNF - Background/Clear Water Extinction Coefficient
        ! *** SWRATNS - Light Extinction for TSS (1/m per g/m3)
        ! *** FSWRATF - Fraction of Solar Rad Absobed in the Surface Layer
        ! *** HTBED2  - Bottom Heat Exchange Coefficient (W/m2/s)
        PSHADE_OLD=-1.0
        DO L=2,LA
          IF(PSHADE_OLD.NE.PSHADE(L))THEN
            IF(SOLSWRT(L).gt.0.01.OR.PSHADE_OLD.LT.-.99)then
              CALL EQUILIBRIUM_TEMPERATURE(SOLSWRT(L),ET,CSHE)
            ENDIF
            PSHADE_OLD=PSHADE(L)
          ENDIF

          ! *** SURFACE HEAT FLUX
          THICK  =HP(L)*DZC(KC)
          TFLUX = CSHE*(ET-TEM(L,KC))/THICK*DELT
          TEM(L,KC) = TEM(L,KC)+TFLUX
          ! *** BEGIN PMC
          ! *** TEMPORARY FIX UNTIL BUILD IN ICE SUB-MODEL INTO THE HEAT SUB-MODEL (COOK INLET)
          IF(ISTRAN(1)>0)THEN
            IF( TEM(L,KC)<-1.3 )THEN
              TEM(L,KC) = -1.3*(SAL(L,KC)/35.)
            ENDIF
          ELSE
            IF( TEM(L,KC)<0.1 )THEN
              TEM(L,KC) = 0.1
            ENDIF
          ENDIF
          ! *** END PMC

          ! *** BOTTOM HEAT FLUX
          THICK = HP(L)*DZC(1)
          TFLUX = HTBED2*(TEMB(L)-TEM(L,1))*DELT
          TEM(L,1) = TEM(L,1)+TFLUX/THICK
          FLUXTB(L)=TFLUX
        ENDDO

        ! *** Distribute Solar Radiation Across Water Column
        IF(SOLSWRT(2).GT.0.1)THEN
          DO L=2,LA
            K=KC
            IF(ISTRAN(6).GT.0.OR.ISTRAN(7).GT.0)THEN
              TSSS_ABOVE=SNDT(L,K)+SEDT(L,K)
              TSS_ABOVE=TSSS_ABOVE
            ELSE
              TSS_ABOVE=0.0
            ENDIF
            IF(ISTRAN(8).GT.0)THEN
              ! *** Water Quality is Active so account for Chlorophyll and POM
              ! *** If using WQ then use the WQ Coefficients
              WQCHLS_ABOVE=WQCHL(L,K)
              WQCHL_ABOVE=WQCHLS_ABOVE
              POMS_ABOVE=WQV(L,K,4)+WQV(L,K,5)
              POM_ABOVE =POMS_ABOVE
              GAMMA = WQKEB(1) + WQKETSS*TSS_ABOVE   +
     &                           WQKECHL*WQCHL_ABOVE +
     &                           WQKEPOM*POM_ABOVE
            ELSE
              GAMMA = SWRATNF + SWRATNS*TSS_ABOVE
            ENDIF

            TOP=GAMMA*HP(L)*(Z(K-1)-1.)
            EXPTOP=EXP(TOP)
            K_ABOVE=1.

            ! *** ENSURE AT LEAST THE FSWRATF FRACTION OF SRO IS ATTENUATED
            IF((1.-EXPTOP).GT.FSWRATF)THEN
              SRON=SOLSWRT(L)*EXPTOP
            ELSE
              SRON=SOLSWRT(L)*(1.0-FSWRATF)
            ENDIF
            DO K = KS,1,-1
              ! *** Net Extinction Coefficient
              K_ABOVE=K_ABOVE+1.
              IF(ISTRAN(6).GT.0.OR.ISTRAN(7).GT.0)THEN
                TSSS_ABOVE=TSSS_ABOVE+SNDT(L,K)+SEDT(L,K)
                TSS_ABOVE=TSSS_ABOVE/K_ABOVE
              ENDIF
              IF(ISTRAN(8).GT.0)THEN
                ! *** Water Quality is Active so account for Chlorophyll
                ! *** If using WQ then use the WQ Coefficients
                POMS_ABOVE=POMS_ABOVE+WQV(L,K,4)+WQV(L,K,5)
                POM_ABOVE=POMS_ABOVE/K_ABOVE

                WQCHLS_ABOVE=WQCHLS_ABOVE+WQCHL(L,K)
                WQCHL_ABOVE=WQCHLS_ABOVE/K_ABOVE
                GAMMA = WQKEB(1) + WQKETSS*TSS_ABOVE   +
     &                             WQKECHL*WQCHL_ABOVE +
     &                             WQKEPOM*POM_ABOVE
              ELSE
                GAMMA = SWRATNF + SWRATNS*TSS_ABOVE
              ENDIF

              BOT=GAMMA*HP(L)*(Z(K-1)-1.)

              ! *** Compute Net Energy
              EXPBOT=EXP(BOT)
              NETRAD(L,K)=SRON*(EXPTOP-EXPBOT)
              TOP=BOT
              EXPTOP=EXPBOT
            ENDDO
            RADBOT(L)=EXPBOT*SRON
          ENDDO

          ! *** NOW FINALIZE THE TEMPERATURE
          DO K=1,KS
            ! *** RHO = 1000.0  Density (kg / m^3)
            ! *** CP  = 4179.0  Specific Heat (J / kg / degC)
            ! *** 0.2393E-6 = 1/RHO/CP --> Conversion from Watts
            C1=DELT*DZIC(K)*0.2393E-6
            DO L=2,LA
              TEMO=TEM(L,K)+HPI(L)*C1*NETRAD(L,K)
              !IF(ABS(TEM(L,K)-TEMO).GT.PMCTOL)THEN
              !  IPMC=0
              !ENDIF
              TEM(L,K)=TEMO
            ENDDO
          ENDDO

          !DO L=2,LA   ! *** PMC PMC PMC
          !  DO K=1,KS
          !    IF(HP(L).LT.1.0.AND.(TEM(L,K)-PMCTOL).GT.TEM(L,K+1)) THEN
          !      IPMC = 0
          !    ENDIF
          !  ENDDO
          !ENDDO

        ENDIF
        ! *** UPDATE BOTTOM
        IF(TBEDIT.GT.0.)THEN
          IF(SOLSWRT(2).GT.0.01.AND.BETAF.GT.0.)THEN
            DO L=2,LA
              TEMB(L)=TEMB(L) +
     &                (0.2393E-6*RADBOT(L)*DELT*BETAF - FLUXTB(L))
     &                      /TBEDTHK(L)
            ENDDO
          ELSE
            DO L=2,LA
              TEMB(L)=TEMB(L) - FLUXTB(L)/TBEDTHK(L)
            ENDDO
          ENDIF
        ENDIF
      ENDIF

      ! *** APPLY DRY CELL CORRECTIONS
      IF(ISDRY.GT.0)THEN
        DO K=1,KC
          DO L=2,LA
            IF(.NOT.LMASKDRY(L))THEN
              TEM(L,K)=TATMT(L)
              ! *** BEGIN PMC
              ! *** TEMPORARY FIX UNTIL BUILD IN ICE SUB-MODEL INTO THE HEAT SUB-MODEL (COOK INLET)
              IF(ISTRAN(1)>0)THEN
                IF( TEM(L,K)<-1.3 )THEN
                  TEM(L,K) = -1.3*(SAL(L,K)/35.)
                ENDIF
              ELSE
                IF( TEM(L,K)<0.1 )THEN
                  TEM(L,K) = 0.1
                ENDIF
              ENDIF
              ! *** END PMC
            ENDIF
          ENDDO
        ENDDO
      ENDIF

  600 FORMAT(4I5,2E12.4)

      RETURN
      END

************************************************************************
**           S U B R O U T I N E   H E A T  E X C H A N G E           **
**                                                                    **
**                      FROM CE-QUAL-W2 (VER 3.1)                     **
**                                                                    **
************************************************************************

      SUBROUTINE HEAT_EXCHANGE

      USE GLOBAL

******* Tupe declaration

        REAL      JDAY
        INTEGER*4 IDAY

******* Allocate/Dimension declaration

        Real ::  MPS_TO_MPH, W_M2_TO_BTU_FT2_DAY
        Real ::  FLUX_BR_TO_FLUX_SI, BTU_FT2_DAY_TO_W_M2
        Real ::  REFL
        Real ::  AFW,BFW,CFW, BCONV
        Real ::  TDEW_F, TAIR_F, WIND_2M
        Real ::  TIMENEXT

******* Data declaration

        DATA MPS_TO_MPH          /2.23714/,
     .       W_M2_TO_BTU_FT2_DAY /7.60796/,
     .       FLUX_BR_TO_FLUX_SI  /0.23659/
        DATA BTU_FT2_DAY_TO_W_M2 /0.1314/
        DATA AFW /9.2/, BFW /0.46/, CFW /2.0/, BCONV /1.520411/
        DATA REFL/0.06/  ! *** SHORTWAVE RADIATION RELFECTED OFF SURFACE

******* Function declaration

        DEG_F(X) = X*1.8+32.0
        DEG_C(X) = (X-32.0)*5.0/9.0

        TIMENEXT=TIMEDAY+1./24.
      RETURN

************************************************************************
**                 S H O R T  W A V E  R A D I A T I O N              **
************************************************************************

      ENTRY SHORT_WAVE_RADIATION(WSPD,TD,TAIR,CLD,ATMPR,SRO,SRON)

******* Input Conversions
        IF(TD.LT.1.1.AND.IRELH(NASER).EQ.1)THEN
          ! *** TD IS RELATIVE HUMIDITY.  CONVERT TO DEW POINT
          ! *** Jensen et al. (1990) ASCE Manual No. 70 (see pages 176 & 177)

          ! *** Ambient vapor pressure in kPa
          VaporP=TD* 0.611*EXP(17.27*TAIR/(TAIR+237.3))

          ! *** Compute dewpoint temperature (Tdew) in C
          Tdew=(116.9+237.3*LOG(VaporP))/(16.78-LOG(VaporP))

          ! *** Compute wet bulb temperature (Tw)
          !PSR=ATMPR/10.  ! ***  Convert millibars to kPa
          !Gamma=0.00066*PSR
          !Delta=4098.*VaporP/(Tdew+237.3)**2
          !Tw=((Gamma*TAIR)+(Delta*Tdew))/(Gamma+Delta)
        ELSE
          TDEW=TD
        ENDIF

        CLD10 = CLD*10.  ! *** CONVERT FROM 0-1 (EFDC) TO 0-10 (W2)
        IF(WINDH<=0.0)WINDH = 10.0     ! *** FIXING THE WIND SPEED MEASUREMENT HEIGHT

******* British units

        TDEW_F   = DEG_F(TDEW)
        TAIR_F   = DEG_F(TAIR)
        WIND_MPH = WSPD*MPS_TO_MPH
        WIND_2M  = WIND_MPH*(LOG(2.0/0.003)/LOG(WINDH/0.003))

******* Shortwave Radiation

        STANDARD = 15.0*INT(DS_LONG/15.0)

        ! *** Day of the Year
        THOUR    = (TIMEDAY-INT(TIMEDAY))*24.0
        IDAY     =  TIMEDAY-INT(TIMEDAY/365.)*365.
        IDAY     =  IDAY+INT(INT(TIMEDAY/365.)/4.)
        JDAY     = REAL(IDAY)
        PMC1     = (2.*PI*(JDAY-1.))/365.
        EQTNEW   =  0.170*SIN(4.*PI*(JDAY-80.)/373.)-
     &              0.129*SIN(2.*PI*(JDAY-8.)/355.)
        H     = 0.2618*(THOUR-(DS_LONG-STANDARD)*0.066667+EQTNEW-12.0)
        DECL =  0.006918-0.399912*COS(PMC1)  +0.070257*SIN(PMC1)-
     &                   0.006758*COS(2*PMC1)+0.000907*SIN(2*PMC1)-
     &                   0.002697*COS(3*PMC1)+0.001480*SIN(3*PMC1)
        SINAL = SIN(DS_LAT*.017453)*SIN(DECL)+COS(DS_LAT*.017453)
     .          *COS(DECL)*COS(H)
        A0    = 57.2958*ASIN(SINAL)
        IF (A0.GT.0.0) THEN
          SRO  = 2.044*A0+0.1296*A0**2-1.941E-3*A0**3+7.591E-6*A0**4
          SRO  = (1.0-0.0065*CLD10**2)*SRO*24.0
          SRO  = SRO*BTU_FT2_DAY_TO_W_M2
          SRON = SRO*(1.0-REFL) ! *** Adjust for surface reflection
        ELSE
          SRO  = 0.0
          SRON = 0.0
        END IF
      RETURN

************************************************************************
**             E Q U I L I B R I U M  T E M P E R A T U R E           **
************************************************************************

      ENTRY EQUILIBRIUM_TEMPERATURE(SRON,ET,CSHE)

******* British units

        ! *** SRON Should already be adjusted for Shading & Reflection
        SRO_BR   = SRON*W_M2_TO_BTU_FT2_DAY

******* Equilibrium temperature and heat exchange coefficient

        ET    = TDEW_F
        TSTAR = (ET+TDEW_F)*0.5
        BETA  = 0.255-(8.5E-3*TSTAR)+(2.04E-4*TSTAR*TSTAR)
        FW    = W_M2_TO_BTU_FT2_DAY*AFW+BCONV*BFW*WIND_2M**CFW
        CSHE  = 15.7+(0.26+BETA)*FW
        RA    = 3.1872E-08*(TAIR_F+459.67)**4
        ETP   = (SRO_BR+RA-1801.0)/CSHE+(CSHE-15.7)
     .            *(0.26*TAIR_F+BETA*TDEW_F)/(CSHE*(0.26+BETA))
        J     = 0
        DO WHILE (ABS(ETP-ET).GT.0.05.AND.J.LT.10)
          ET    = ETP
          TSTAR = (ET+TDEW_F)*0.5
          BETA  = 0.255-(8.5E-3*TSTAR)+(2.04E-4*TSTAR*TSTAR)
          CSHE  = 15.7+(0.26+BETA)*FW
          ETP   = (SRO_BR+RA-1801.0)/CSHE+(CSHE-15.7)
     .            *(0.26*TAIR_F+BETA*TDEW_F)/(CSHE*(0.26+BETA))
          J     = J+1
        END DO

******* SI units

        ! *** RHO = 1000.0  Density (kg / m^3)
        ! *** CP  = 4179.0  Specific Heat (J / kg / degC)
        ! *** 0.2393E-6 = 1/RHO/CP
        ET   = DEG_C(ET)
        CSHE = CSHE*FLUX_BR_TO_FLUX_SI*0.2393E-6

        IF(DEBUG)THEN
          IF(TIMEDAY.GT.TIMENEXT)THEN
            TD_C=DEG_C(TDEW_F)
            TA_C=DEG_C(TAIR_F)
            WRITE(77,999)TIMEDAY,SRON,ET,TD_C,TA_C,TDEW_F,TAIR_F,FW,CSHE
  999       FORMAT(F11.4,8F9.2,E12.4)
            TIMENEXT = TIMENEXT+1./24.
          ENDIF
        ENDIF
      RETURN

      END


text:
# 1 "/data/home/scott/test/gfortran/src/gfortran.dg/dev_null.F90"
# 1 "<built-in>"
# 1 "<command line>"
# 1 "/data/home/scott/test/gfortran/src/gfortran.dg/dev_null.F90"
! { dg-do run }
! { dg-options "-std=legacy" }
!
! pr19478 read from /dev/null
! Thomas.Koenig@online.de

      character*20 foo
      open(10,file="/dev/null")
      write(10,'(A)') "Hello"
      rewind(10)
      read(10,'(A)',end=100) foo
      call abort
 100  continue
      end

max_stars_repo_name: maurizioabba/rose

max_issues_repo_name: sujankh/rose-matlab

max_forks_repo_name: sujankh/rose-matlab


text:
great_ne(x1,v1).
great_ne(b1,i1).
great_ne(p1,m1).
great_ne(jj1,dd1).
great_ne(aa1,h1).
great_ne(p1,w1).
great_ne(dd1,x1).
great_ne(x1,c1).
great_ne(aa1,c1).
great_ne(k1,v1).
great_ne(cc1,w1).
great_ne(cc1,bb1).
great_ne(kk1,dd1).
great_ne(cc1,h1).
great_ne(q1,f1).
great_ne(cc1,i1).
great_ne(aa1,dd1).
great_ne(hh1,b1).
great_ne(aa1,m1).
great_ne(p1,v1).
great_ne(w1,ff1).
great_ne(k1,ff1).
great_ne(ee1,hh1).
great_ne(k1,a1).
great_ne(ee1,h1).
great_ne(bb1,ff1).
great_ne(u1,x1).
great_ne(t1,b1).
great_ne(x1,w1).
great_ne(jj1,v1).
great_ne(v1,ff1).
great_ne(z1,b1).
great_ne(q1,bb1).
great_ne(ff1,i1).

max_stars_repo_name: manoelfranca/cilppp

max_issues_repo_name: manoelfranca/cilppp

max_forks_repo_name: manoelfranca/cilppp


text:
      PROGRAM TPWRZT
C
C Define the error file, the Fortran unit number, the workstation type,
C and the workstation ID to be used in calls to GKS routines.
C
C     PARAMETER (IERRF=6, LUNIT=2, IWTYPE=1,  IWKID=1)   ! NCGM
C     PARAMETER (IERRF=6, LUNIT=2, IWTYPE=8,  IWKID=1)   ! X Windows
C     PARAMETER (IERRF=6, LUNIT=2, IWTYPE=11, IWKID=1)   ! PDF
C     PARAMETER (IERRF=6, LUNIT=2, IWTYPE=20, IWKID=1)   ! PostScript
C
      PARAMETER (IERRF=6, LUNIT=2, IWTYPE=1,  IWKID=1)
C
C OPEN GKS, OPEN WORKSTATION OF TYPE 1, ACTIVATE WORKSTATION
C
      CALL GOPKS (IERRF, ISZDM)
      CALL GOPWK (IWKID, LUNIT, IWTYPE)
      CALL GACWK (IWKID)
C
C INVOKE DEMO DRIVER
C
      CALL TPWRZT1(IERR)
C
C DEACTIVATE AND CLOSE WORKSTATION, CLOSE GKS.
C
      CALL GDAWK (IWKID)
      CALL GCLWK (IWKID)
      CALL GCLKS
C
      STOP
      END
C
      SUBROUTINE TPWRZT1 (IERROR)
C
C PURPOSE                To provide a simple demonstration of
C                        entry PWRZT with the THREED utility.
C
C USAGE                  CALL TPWRZT1 (IERROR)
C
C ARGUMENTS
C
C ON OUTPUT              IERROR
C                          An integer variable
C                          = 0, if the test was successful,
C                          = 1, otherwise
C
C I/O                    If the test is successful, the message
C
C               PWRZT TEST EXECUTED--SEE PLOTS TO CERTIFY
C
C                        is printed on unit 6.  In addition, 1
C                        frame is produced on the machine graphics
C                        device.  In order to determine if the test
C                        was successful, it is necessary to examine
C                        the plot.
C
C PRECISION              Single
C
C REQUIRED ROUTINES      PWRZT, THREED
C
C REQUIRED GKS LEVEL     0A
C
C LANGUAGE               FORTRAN 77
C
C ALGORITHM              The THREED package is called to establish a
C                        3-D projection onto 2 space, and to draw the
C                        axis lines.  PWRZT is then called to label
C                        the axes for a 3 space plot.
C
C EYE contains the (U,V,Z) coordinate of the eye position.
C
      REAL            EYE(3)
      DATA EYE(1), EYE(2), EYE(3) /3.5, 3.0, 5.0/
C
C Initialize the error parameter.
C
      IERROR = 1
C
C Select normalization transformation number 0.
C
      CALL GSELNT (0)
C
C A call to SET3 establishes the mapping of 3 space coordinates onto
C the coordinate system of the graphics device.
C
      CALL SET3 (.1,.9,.1,.9,0.,1.,0.,1.,0.,1.,EYE)
C
C Draw the 3 space axes.
C
      CALL LINE3 (0.,0.,0.,0.,0.,1.)
      CALL LINE3 (0.,0.,0.,0.,1.,0.)
      CALL LINE3 (0.,0.,0.,1.,0.,0.)
C
C PWRZT is used to label each of the axes and the plot
C
      ICNT = 0
      ISIZE = 30
      LINE = 2
      ITOP = 3
      CALL PWRZT (0.,.5,.1,'V-AXIS',6,ISIZE,LINE,ITOP,ICNT)
C
      LINE = -1
      ITOP = 3
      CALL PWRZT (.5,0.,.1,'U-AXIS',6,ISIZE,LINE,ITOP,ICNT)
C
      LINE = 3
      ITOP = -2
      CALL PWRZT (0.,.1,.5,'Z-AXIS',6,ISIZE,LINE,ITOP,ICNT)
C
      LINE = 2
      ITOP = -1
      ISIZE = 30
      ICNT = -1
      CALL PWRZT (.5,.2,0.,'DEMONSTRATION OF PWRZT WITH THREED',
     1            34,ISIZE,LINE,ITOP,ICNT)
      CALL FRAME
C
      IERROR = 0
      WRITE (6,1001)
C
      RETURN
 1001 FORMAT (' PWRZT TEST EXECUTED--SEE PLOT TO CERTIFY')
C
      END

max_stars_repo_name: tenomoto/ncl

max_issues_repo_name: tenomoto/ncl

max_forks_repo_name: tenomoto/ncl


! { dg-do compile }
! { dg-options "-fcoarray=single" }
!

subroutine cont1(x) ! { dg-error "has the CONTIGUOUS attribute but is not an array pointer or an assumed-shape or assumed-rank array" }
  type t
  end type t
  class(t), contiguous, allocatable :: x(:)
end

subroutine cont2(x) ! { dg-error "has the CONTIGUOUS attribute but is not an array pointer or an assumed-shape or assumed-rank array" }
  type t
  end type t
  class(t), contiguous, allocatable :: x(:)[:]
end

subroutine cont3(x, y)
  type t
  end type t
  class(t), contiguous, pointer :: x(:)
  class(t), contiguous :: y(:)
end

function func() ! { dg-error "shall not be a coarray or have a coarray component" }
  type t
  end type t
  class(t), allocatable :: func[*]
end

function func2() ! { dg-error "must be dummy, allocatable or pointer" }
  type t
    integer, allocatable :: caf[:]
  end type t
  class(t) :: func2a ! { dg-error "CLASS variable 'func2a' at .1. must be dummy, allocatable or pointer" }
  class(t) :: func2
end

subroutine foo1(x1) ! { dg-error "Coarray variable 'x1' at .1. shall not have codimensions with deferred shape" }
  type t
  end type t
  type(t) :: x1(:)[:]
end

subroutine foo2(x2) ! { dg-error "Coarray variable 'x2' at .1. shall not have codimensions with deferred shape" }
  type t
  end type t
  type(t) :: x2[:]
end

! DITTO FOR CLASS

subroutine foo3(x1) ! { dg-error "Coarray variable 'x1' at .1. shall not have codimensions with deferred shape" }
  type t
  end type t
  class(t) :: x1(:)[:]
end

subroutine foo4(x2) ! { dg-error "Coarray variable 'x2' at .1. shall not have codimensions with deferred shape" }
  type t
  end type t
  class(t) :: x2[:]
end

subroutine bar1(y1) ! { dg-error "Allocatable coarray variable 'y1' at .1. must have deferred shape" }
  type t
  end type t
  type(t), allocatable :: y1(:)[5:*]
end

subroutine bar2(y2) ! { dg-error "Allocatable coarray variable 'y2' at .1. must have deferred shape" }
  type t
  end type t
  type(t), allocatable :: y2[5:*]
end

subroutine bar3(z1) ! { dg-error "Allocatable coarray variable 'z1' at .1. must have deferred shape" }
  type t
  end type t
  type(t), allocatable :: z1(5)[:]
end

subroutine bar4(z2) ! { dg-error "Allocatable array 'z2' at .1. must have a deferred shape" }
  type t
  end type t
  type(t), allocatable :: z2(5)
end subroutine bar4

subroutine bar5(z3) ! { dg-error "Array pointer 'z3' at .1. must have a deferred shape" }
  type t
  end type t
  type(t), pointer :: z3(5)
end subroutine bar5

! DITTO FOR CLASS

subroutine bar1c(y1) ! { dg-error "Allocatable coarray variable 'y1' at .1. must have deferred shape" }
  type t
  end type t
  class(t), allocatable :: y1(:)[5:*]
end

subroutine bar2c(y2) ! { dg-error "Allocatable coarray variable 'y2' at .1. must have deferred shape" }
  type t
  end type t
  class(t), allocatable :: y2[5:*]
end

subroutine bar3c(z1) ! { dg-error "Allocatable coarray variable 'z1' at .1. must have deferred shape" }
  type t
  end type t
  class(t), allocatable :: z1(5)[:]
end

subroutine bar4c(z2) ! { dg-error "Allocatable array 'z2' at .1. must have a deferred shape" }
  type t
  end type t
  class(t), allocatable :: z2(5)
end subroutine bar4c

subroutine bar5c(z3) ! { dg-error "Array pointer 'z3' at .1. must have a deferred shape" }
  type t
  end type t
  class(t), pointer :: z3(5)
end subroutine bar5c

subroutine sub()
  type t
  end type
  type(t) :: a(5)
  class(t), allocatable :: b(:)
  call inter(a)
  call inter(b)
contains
  subroutine inter(x)
    class(t) :: x(5)
  end subroutine inter
end subroutine sub

subroutine sub2()
  type t
  end type
  type(t) :: a(5)
contains
  subroutine inter(x)
    class(t) :: x(5)
  end subroutine inter
end subroutine sub2

subroutine sub3()
  type t
  end type
contains
  subroutine inter2(x) ! { dg-error "must have a deferred shape" }
    class(t), pointer :: x(5)
  end subroutine inter2
end subroutine sub3


text:
locate port P at "port_ok.fpp"
port P

ext: fpp

max_stars_repo_name: kevin-f-ortega/fpp

max_issues_repo_name: kevin-f-ortega/fpp

max_forks_repo_name: kevin-f-ortega/fpp